Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 915,319 (A)7→8 intergenic (‑462/+33) lysO ← / ← aqpZ L‑lysine exporter/water channel AqpZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,3121.A93.5% 99.4 / ‑0.2 31intergenic (‑455/+40)lysO/aqpZL‑lysine exporter/water channel AqpZ
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (20/9);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATG  >  NC_000913/915196‑915446
                                                                                                                      |                                                                                                                                       
tGCAGGGGGCGGCGGAGATTGGGCGCCCCCGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGc                                                                                                                  <  1:285120/140‑1 (MQ=37)
         cggcggcgATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAAGATCTTAAAAAAGAGCCCTGACGTTACGCCAGGCCTTc                                                                                                         >  1:299507/1‑140 (MQ=37)
             ggcgATTGGCCGCCCCCTGCGCGCTCCTTGCCCCAATGGAAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg                                                                                                     <  2:273150/140‑1 (MQ=38)
                        gccCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa                                                                                          >  1:6840/1‑140 (MQ=255)
                           cccTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTATCTCATGAATGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc                                                                                       <  2:258741/140‑1 (MQ=255)
                                         tGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt                                                                         >  1:127817/1‑140 (MQ=255)
                                         tGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGt                                                                         >  2:126890/1‑140 (MQ=255)
                                                            tGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCg                                                      >  1:559671/1‑140 (MQ=255)
                                                              cTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGAt                                                    >  2:781877/1‑140 (MQ=255)
                                                              cTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAACAGGGTCCGGTAAATCAGACCACCGAt                                                    >  2:583156/1‑140 (MQ=255)
                                                               ttAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATa                                                   >  2:792758/1‑140 (MQ=255)
                                                                 aGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAAt                                                 >  1:845635/1‑140 (MQ=255)
                                                                 aGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCGTTTCCAGAAGGGTCGGGTAAAGCAGACCACCGATAAt                                                 >  1:416544/1‑140 (MQ=255)
                                                                   cTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCATTTGCAGGAGGGTCCGATAAAGCAGACCACCGATAATgc                                               >  2:236932/1‑140 (MQ=255)
                                                                                  caacaaGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGca                                <  1:14743/140‑1 (MQ=255)
                                                                                  caacaaGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGca                                <  1:126890/140‑1 (MQ=255)
                                                                                  caacaaGATGATGAATCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGca                                <  2:6840/140‑1 (MQ=255)
                                                                                           gatgaATCAACATATCACAACATCTCAAAAAAAAGGCGTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCaga                       <  1:236932/140‑1 (MQ=255)
                                                                                               aaTCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAgaagaa                   >  1:70072/1‑140 (MQ=255)
                                                                                                 tCAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAAcc                 >  2:812511/1‑140 (MQ=255)
                                                                                                  cAACATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACca                <  2:559671/140‑1 (MQ=255)
                                                                                                   aaCATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACcac               <  2:845635/140‑1 (MQ=255)
                                                                                                    aCATATAACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACcaca              >  1:847612/1‑140 (MQ=255)
                                                                                                         taACAACATCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgttgt         >  2:487255/1‑140 (MQ=255)
                                                                                                               caTCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgtt           <  2:685154/132‑1 (MQ=255)
                                                                                                               caTCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAgtt           >  1:685154/1‑132 (MQ=255)
                                                                                                               caTCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAAt   >  2:649182/1‑140 (MQ=255)
                                                                                                               caTCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAAt   >  2:45346/1‑140 (MQ=255)
                                                                                                               caTCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAAt   >  1:178017/1‑140 (MQ=255)
                                                                                                                aTCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATg  >  2:567054/1‑140 (MQ=255)
                                                                                                                aTCTTAAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATg  >  1:790038/1‑140 (MQ=255)
                                                                                                                      |                                                                                                                                       
TGCAGGGGGCGGCGGCGATTGGCCGCCCCCTGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATG  >  NC_000913/915196‑915446

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