Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,819,318 G→A intergenic (‑124/+79) ves ← / ← spy environmental stress‑induced protein Ves/ATP‑independent periplasmic protein‑refolding chaperone Spy

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,819,3180GA80.0% 14.5 / 3.7 20intergenic (‑124/+79)ves/spyenvironmental stress‑induced protein Ves/ATP‑independent periplasmic protein‑refolding chaperone Spy
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (0/16);  total (4/16)
Fisher's exact test for biased strand distribution p-value = 2.06e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-03

TCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTAT  >  NZ_CP009273/1819230‑1819400
                                                                                         |                                                                                  
tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga                                                                                   >  1:171346/1‑90 (MQ=255)
                    aTCGGCACCTGGGGTGTTTTACTTTAGTGGTGATCGGTTTTTTTTTACAACTTTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAc                                                               <  1:212392/90‑1 (MQ=255)
                     ccGGCAACTGGGGGATTTTAAGTTAGGGGTGATCGGGTTTTTTTTAAACCTGTTGGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAcc                                                              <  1:687086/89‑1 (MQ=255)
                           aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                        >  2:106117/1‑90 (MQ=255)
                                 gttaTTTACGTTAGTGTTGACCAGTTTTTTTTTACAACTTTTTGAAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                  <  2:119699/87‑1 (MQ=255)
                                 gtgtTTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                  <  2:638712/90‑1 (MQ=255)
                                 gtgtTTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                  <  2:313130/90‑1 (MQ=255)
                                 gtattttACGTTGGTGGTGATCGGGTTTTTTTTACAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                  <  2:601411/86‑1 (MQ=255)
                                 gTGATTTCCGTAAGGGGTGACCAGCTTTTTTTTACAACTTTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa                                                  <  1:391707/90‑1 (MQ=255)
                                 gTGATTTACTTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                  <  2:37933/90‑1 (MQ=255)
                                 gTGATTTACTTTAGGGGTGATCGGTTTTTTTTTACAATTTTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                  <  1:66680/90‑1 (MQ=255)
                                 gTGATTTACGTTAGTGGTGATCAGGTTTTTTTTACAACTGTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa                                                  <  2:130156/90‑1 (MQ=255)
                                  tGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg                                                 <  2:272135/90‑1 (MQ=255)
                                      ttATGTTGGTGGTGATCGGTTTTTTTTTAAATTTGTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                             <  1:606297/90‑1 (MQ=255)
                                      ttACTTTGGGGGTGACAGGTTTTTTTTAAAAATTTTTGGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg                                             <  1:516571/90‑1 (MQ=255)
                                      ttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGGAg                                             >  2:604359/1‑90 (MQ=255)
                                         ctttAGTGGGGATCAGGTTTTTTTTACAACTGTTTGACAAGATAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg                                          <  2:300580/88‑1 (MQ=255)
                                             aGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACaa                                      >  2:545141/1‑90 (MQ=255)
                                                    agcAGGTTTTTTTTAAAACTTTTTGAGAAGAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                               <  2:301987/88‑1 (MQ=255)
                                                    aTCAGGTTTTTTTAACACTTTTTTGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                               <  2:457805/90‑1 (MQ=255)
                                                                                  gagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGGAAGattat  >  1:233955/1‑90 (MQ=255)
                                                                                         |                                                                                  
TCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTAT  >  NZ_CP009273/1819230‑1819400

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: