Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 897,968:1 +A 19.4% intergenic (+39/‑21) 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,9681.A19.4% 86.7 / 26.3 31intergenic (+39/‑21)potI/ybjOputrescine ABC transporter permease/DUF2593 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (12/13);  new base A (3/3);  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 = 8.12e-01

AACATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGC  >  NC_000913/897821‑898109
                                                                                                                                                    |                                                                                                                                             
aaCATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCgcgc                                                                                                                                              <  3:423686/149‑1 (MQ=255)
  cATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCAt                                                                                                                                            <  4:1041132/149‑1 (MQ=255)
   aTTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCcgcgacat                                                                                                                                            <  1:1117698/149‑5 (MQ=255)
   aTTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATg                                                                                                                                           <  2:460969/149‑1 (MQ=255)
                gcgGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATGCGCGGCAttat                                                                                                                               >  1:79771/1‑146 (MQ=255)
                gcgGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGACATGCGCGGCAttat                                                                                                                               <  2:79771/148‑3 (MQ=255)
                                      tCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCa                                                                                                                          <  1:706390/131‑1 (MQ=255)
                                      tCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCa                                                                                                                          >  2:706390/1‑131 (MQ=255)
                                              tATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTg                                                                                                         >  3:378585/1‑140 (MQ=255)
                                              tATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTg                                                                                                         <  4:378585/140‑1 (MQ=255)
                                                tCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                              <  4:946258/149‑1 (MQ=255)
                                                tCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                              >  3:946258/1‑149 (MQ=255)
                                                           gcgcAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCa                                                                                                                          >  3:934019/1‑110 (MQ=255)
                                                           gcgcAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCa                                                                                                                          <  4:934019/110‑1 (MQ=255)
                                                                gAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                              <  3:430442/133‑1 (MQ=255)
                                                                gAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTaa                                                                                              >  4:430442/1‑133 (MQ=255)
                                                                gAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACAt                                                                                      <  4:314807/141‑1 (MQ=255)
                                                                gAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACAt                                                                                      >  3:314807/1‑141 (MQ=255)
                                                                 aaaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACACAATGCCGCGACATGCGCGGCAttat                                                                                                                               >  4:391080/1‑97 (MQ=38)
                                                                 aaaaaCAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACACAATGCCGCGACATGCGCGGCAttat                                                                                                                               <  3:391080/99‑3 (MQ=38)
                                                                                gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGa                                                                                                                   <  2:138024/96‑1 (MQ=255)
                                                                                gcgACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGa                                                                                                                   >  1:138024/1‑96 (MQ=255)
                                                                                         aGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGaa                                                                <  4:308450/138‑1 (MQ=255)
                                                                                         aGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGaa                                                                >  3:308450/1‑138 (MQ=255)
                                                                                                       ggCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                          <  3:1238590/82‑1 (MQ=255)
                                                                                                       ggCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGtt                                                                                                          >  4:1238590/1‑82 (MQ=255)
                                                                                                              gACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCt                                                         <  1:616646/124‑1 (MQ=255)
                                                                                                              gACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCt                                                         >  2:616646/1‑124 (MQ=255)
                                                                                                                     aaaTTTGCGAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCattat                         >  3:1348151/1‑149 (MQ=255)
                                                                                                                                          caaTGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                  <  2:557715/102‑1 (MQ=255)
                                                                                                                                          caaTGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTgg                                                  >  1:557715/2‑103 (MQ=255)
                                                                                                                                            aTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAg                                             >  2:35282/1‑106 (MQ=255)
                                                                                                                                            aTGCCGCGACATGCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAg                                             <  1:35282/106‑1 (MQ=255)
                                                                                                                                                cgcgCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc  <  2:123213/145‑1 (MQ=255)
                                                                                                                                                cgcgCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGc  >  1:123213/1‑145 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
AACATTGATCCTCGGCGCGGTCGGAATTGTCGGATTTATCGCCTGGTATCTGATGGCTCGCGCAGAAAAACAGCGGATACGCGACATCCAGCGTGCAAGACGTGGCTGAAGACACTAAAATTTGCCAACCTGGCTACATAATGCCGCGCATGTCGCGGCATTGTTTTCATGGAAGACGAAACGTTGGGATTTTTTAAGAAAACATCTTCATCTCATGCTCGCCTGAATGTGCCTGCGCTGGTGCAGGTGGCGGCGCTCGCCATTATTATGATCCGTGGCCTCGACGTGC  >  NC_000913/897821‑898109

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