Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,600,081 C→T W43* (TGG→TGA)  lsrK ← autoinducer‑2 kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,600,0810CT87.5% 13.0 / ‑5.3 8W43* (TGG→TGA) lsrKautoinducer‑2 kinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (2/5);  total (2/6)
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

CGTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTT‑GTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCA‑TCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCC  >  NC_000913/1599857‑1600306
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  
ngTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTT‑GTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTg                                                                                                                                                                                                          >  1:89950/2‑251 (MQ=255)
                                                       tatataAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTT‑GTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCan                                                                                                                                                    <  1:30701/250‑2 (MQ=255)
                                                                           cGCATCGAACAGGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCGAGTTTTT‑GTTGAGATCAAAGTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTg                                                                                                                                                                                                                         <  1:156334/161‑1 (MQ=255)
                                                                                                                                           ggCGCATACCCTCGCCCGCCAGTTGCCAGTTTTTTGTT‑AGATCAAATTCCATAGACCCAGGAACGTCCTGTACTGCCAGATGCCGTCACGCCGCTTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATTCTGCCGGTGCCTGCATCCATCGCCT‑GCAGGTAGTACTTTTATTCTGAAAGGGTAAAGAGTCGATCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCctc                                                              <  1:266945/251‑1 (MQ=255)
                                                                                                                                           ggCGCATACACTCTCACGCCAGTTGCCATTTTTT‑GTTGAGAGCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCA‑TCAGGTAGTACTTTTATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCctc                                                              <  1:311678/251‑1 (MQ=255)
                                                                                                                                           ggCGCATACACTCCCACGCCAGTTGCCAGTTTTT‑GTTAAGATCAAATTCCATAGAACCAGGAACGTCCTGTTCTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACATCCCGAATACTTCCGGTGCCTGCATCCAGCGCCA‑TCAGGTAGTACTTTGAGGCTGAAATGTTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCctc                                                              <  1:266937/251‑1 (MQ=255)
                                                                                                                                              gcATACACTCACACGCCAGTTGCCAGTTTTT‑GTTCAGAAAAAATTCCATAGTAACAGGAACGTCCGGTACTGCCAGATGCCGTCACTCCGCCTATCCCACTGCTAttt                                                                                                                                                                                                           >  1:308508/1‑108 (MQ=255)
                                                                                                                                               cacacacTCACACGCCAGTTTCCAGTTTTT‑GTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGTCACGCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGc                                                                                                                                                              <  1:131993/149‑1 (MQ=255)
                                                                                                                                                                                                                             gccgtcACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATCTCACAGCACGAATACTTCCGGTGCCTGCCTCCAGCGCCA‑TCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCATCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCcatcat                                             >  1:242428/1‑187 (MQ=255)
                                                                                                                                                                                                                             gccgtcACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCNN‑TCAGgtnnt                                                                                                                                     >  2:9964/1‑99 (MQ=255)
                                                                                                                                                                                                                             gccgtcACTCCGCCTGTCCCACTGCTATTTGATTGCCTCCCAGGTCGAATACCACAGCCCGAATACTCCCGGTGCCTGCATCCAG‑GCCAATCAGGTAGTATTTGGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATCCCCGGTAATAAGGTCATGCAAATTTAACTCCGTAAAATcgcc  <  1:9964/230‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                  
CGTGCCGCTCTGGCATCCACATTGGCGCAGGCCCAGATCGGGGCTCCTTCATTATTATATAAAACAATGCCTTCACGCATCGAACATGCCGAAACGGCAGCGATATACTCCGGGGCTATGCCGGCGTTGTGCAGCGCCTGGCGCATACACTCACACGCCAGTTGCCAGTTTTT‑GTTGAGATCAAATTCCATAGAACCAGGAACGTCCGGTACTGCCAGATGCCGCCACTCCGCCTGTCCCACTGCTATTTGATTGCCTTCCAGGTCGAATATCACAGCCCGAATACTTCCGGTGCCTGCATCCAGCGCCA‑TCAGGTAGTACTTTGATTCTGAAAGGGTAAAGAGTCGAGCCATTAGATTATCCTCGGCTTATGCGTTCTGTCTGTTCCTCTATACGTTCTCCATCATTCCCGGTAATAAGGTCATGCAAATTTAACTACGTAAAATCGCC  >  NC_000913/1599857‑1600306

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑