Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,092,703 A→G F41L (TTC→CTC)  frvA ← putative PTS enzyme IIA component FrvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,092,7030AG95.8% 73.5 / ‑3.8 24F41L (TTC→CTC) frvAputative PTS enzyme IIA component FrvA
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (7/16);  total (7/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

TTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTT‑TGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGG‑AGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAATG  >  NC_000913/4092465‑4092949
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
gtctgcatcaGTGTGCGGATGATGTTGCCCTGATCCTCTTCGCCGCTTGTG‑GGCACGCCGAGTCAGATCCAGCAATTGTCGGCGTTGCCATCGCTGGCTGTCCAGTCAATAGCCTGCGCTTTGCGTGCGAATAATCCGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGATCCAAATCCAGTAGATTGCATTTGTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGGGTGAGTc                                                                                                                                                                                                                                             <  2:406617/243‑1 (MQ=255)
         atggtgCCGATGATTTTGACCTGATCCTCTTCGCCGCTT‑TGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATaaa                                                                                                                                                                                                                                    <  1:66079/248‑1 (MQ=255)
            gtgCCGATGATTTTGACCTGATCCTCTTCGCCGCTT‑TGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAAcc                                                                                                                                                                                                                                 >  1:113128/1‑251 (MQ=255)
                cGATGATTTTGACCTGATCCTCTTCGCCGCTT‑TGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGTTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGttt                                                                                                                                                                                                                             >  1:386732/1‑251 (MQ=255)
                 gatgatTTTGACCTGATCCTCTTCGCCGCTT‑TGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGtttt                                                                                                                                                                                                                            >  1:166867/1‑251 (MQ=255)
                                                     gCACGCCGAGGCCGATCAAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCACTAGCATGGGCTTTGCAGGCGAATAATACGAACGGCTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACAGGCACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGCGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTTGCGCCCACTGCTTCTGAACGGa                                                                                                                                                                                          >  2:413281/1‑250 (MQ=255)
                                                                cAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGtttt                                                                                                                                                                                                                            >  1:238107/1‑205 (MQ=255)
                                                                cAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGtttt                                                                                                                                                                                                                            <  2:238107/205‑1 (MQ=255)
                                                                  gATCCAGCAATTTACGTCTTCGCCAGCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCTGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGACAACTGCTTCAGAACGGAATAAGCGCCATTGc                                                                                                                                                                            <  2:116797/251‑1 (MQ=255)
                                                                      cAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTg                                                                                                                                                                        <  1:197034/251‑1 (MQ=255)
                                                                              gACGTCTTCGCCAACGCTGGCTTTCCAGTCAATAGCCTGGGCGTTGCGGGCGAATAATACGAACGGTTGGTTAACGCAGGCGCTGTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCa                                                                                                                                                                <  2:202770/251‑1 (MQ=255)
                                                                                     tCGCCATCGCTGGCTGTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATGTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGGCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGTCCTGAATATTCAGGTCa                                                                                                                                                           <  2:412835/249‑1 (MQ=255)
                                                                                      cGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACCCCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAAt                                                                                                                                                         <  1:25343/250‑1 (MQ=255)
                                                                                                        aGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCt                                                                                                                                      <  1:183070/251‑1 (MQ=255)
                                                                                                                   tGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGCGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCCTTAATc                                                                                                                            >  2:342172/1‑250 (MQ=255)
                                                                                                                         ttGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGAcgc                                                                                                                      <  2:223005/250‑3 (MQ=255)
                                                                                                                                  aataatACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatat                                                                                                            <  2:162743/251‑1 (MQ=255)
                                                                                                                                     aataCGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCTGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatatat                                                                                                          <  1:113172/250‑1 (MQ=255)
                                                                                                                                                 tgttTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTg                                                                                              <  1:56004/250‑1 (MQ=255)
                                                                                                                                                                    tACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCAGTTTTGCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGG‑AGCGATTAAg                                                                          <  1:92387/250‑1 (MQ=255)
                                                                                                                                                                                                   tCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGG‑AGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGAc                                          >  1:124427/1‑251 (MQ=255)
                                                                                                                                                                                                         agagTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAGCTTGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGG‑AGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTg                                    <  1:59802/251‑1 (MQ=255)
                                                                                                                                                                                                                                gcagGGTTTGCAGTAGCTGGGGTGAGTCGGTGATAAAACCGTTTTGTAACGCTTTTGTCGCCAACTGCTTCAGAACGGAATAAGCGACATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATTGGACTCCATAAAATGTATTTTAATCAATTGTAATG‑TGGTAGCGATTAAGTGATGAGAGGCGATAGGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGAttt             <  1:176552/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aGGAGCTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGG‑AGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAATg  <  1:342172/251‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       
TTGCGACACAGTGTGCCGATGATTTTGACCTGATCCTCTTCGCCGCTT‑TGCGGCACGCCGAGGCAGATCCAGCAATTGACGTCTTCGCCATCGCTGGCTTTCCAGTCAATAGCCTGCGCTTTGCGGGCGAATAATACGAACGGTTGTTTAACGCAGGCGCTTTTACCGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGG‑AGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAGGTCAGATTTGCGGAGTAATG  >  NC_000913/4092465‑4092949

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