Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 897,973 Δ1 bp 12.9% intergenic (+44/‑16) potI → / → ybjO putrescine ABC transporter permease/DUF2593 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913897,9730T.12.9% 98.1 / 9.5 31intergenic (+44/‑16)potI/ybjOputrescine ABC transporter permease/DUF2593 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (13/14);  new base . (2/2);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

TCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTT  >  NC_000913/897829‑898118
                                                                                                                                                 |                                                                                                                                                 
tCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTcgcg                                                                                                                                               <  4:53059/149‑1 (MQ=255)
tCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTcgcg                                                                                                                                               >  4:39879/1‑149 (MQ=255)
tCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTcgcg                                                                                                                                               <  4:22316/149‑1 (MQ=255)
    cGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAt                                                                                                                                           <  1:277041/149‑1 (MQ=255)
         cgGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGttt                                                                                                                                      <  1:176842/149‑1 (MQ=255)
                    gTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGAc                                                                                                                           <  4:304504/149‑1 (MQ=255)
                     tCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAt                                                                                                                                           <  2:186103/132‑1 (MQ=255)
                     tCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATG‑CGCGGCAt                                                                                                                                           >  1:186103/1‑132 (MQ=255)
                     tCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACg                                                                                                                          >  1:132527/1‑149 (MQ=255)
                     tCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACg                                                                                                                          >  4:277789/1‑149 (MQ=255)
                                    ggTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGAtttt                                                                                                           <  2:132527/149‑1 (MQ=255)
                                           gATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaa                                                                                                    <  4:400782/149‑1 (MQ=255)
                                                cTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACAt                                                                                               <  4:249989/149‑1 (MQ=255)
                                                 tCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATc                                                                                              >  4:257672/1‑149 (MQ=255)
                                                            aaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaaa                                                                                                  >  2:421264/1‑134 (MQ=255)
                                                            aaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGaaaa                                                                                                  <  1:421264/134‑1 (MQ=255)
                                                                   gATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCt                                                                            >  2:9030/1‑149 (MQ=255)
                                                                                         aGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAg                                                      <  1:9030/149‑1 (MQ=255)
                                                                                                 cTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggcggc                                              >  3:356223/1‑149 (MQ=255)
                                                                                                 cTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTggcggc                                              <  3:39879/149‑1 (MQ=255)
                                                                                                  tGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGGGCTGGTGCAGGTGgcggcg                                             >  3:354500/1‑149 (MQ=255)
                                                                                                     aGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                       <  2:438146/88‑1 (MQ=255)
                                                                                                     aGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                                       >  1:438146/1‑88 (MQ=255)
                                                                                                           tAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCATGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCatt                                    >  3:299199/1‑149 (MQ=255)
                                                                                                             aaaTTTGCCAACCTGGCTACATAATGCCCCG‑CATGTCGCGGCATTGTTTTCATGGAAGTCGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCattat                                  <  4:408512/149‑1 (MQ=255)
                                                                                                                                aacaaTGCCGCGACATG‑CGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGAc               >  4:394727/4‑149 (MQ=255)
                                                                                                                                     tGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGt                                                          >  4:164367/1‑101 (MQ=255)
                                                                                                                                     tGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGt                                                          <  3:164367/101‑1 (MQ=255)
                                                                                                                                     tGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCt          >  2:170857/1‑149 (MQ=255)
                                                                                                                                          cg‑cATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCtgatga     <  3:257672/149‑1 (MQ=255)
                                                                                                                                              aTGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGAttt  >  3:344466/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCG‑CATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGCTGATGATTT  >  NC_000913/897829‑898118

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑