Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,331,890 C→T L463L (CTC→CTT)  proP → osmolyte:H(+) symporter ProP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,331,8900CT93.8% 47.7 / ‑4.2 16coding (1389/1503 nt)proPosmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/8);  minor base . (0/1);  total (7/9)
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

TTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTCGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAA  >  NC_000913/4331686‑4332120
                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
ttCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGAc                                                                                                                                                                                           <  1:34227/250‑1 (MQ=255)
         ggTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGAtt                                                                                                                                                                                  <  1:841/250‑1 (MQ=255)
              ccGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGGTGGGTTAATTACCGGCGTAAACATGAAAAAGAAGGCAAAAACG‑CCGTTGTAAAGTGCGGCACCCGCGGCGGCAGAAATACAGGGAAG‑GGAGGAAAATTTTTTT‑GGGCATTACGATAAAAACGAGCAGAAAATCGATTATATTTTACACGAGATTTCCgg                                                                                                                                                                             >  1:223821/1‑249 (MQ=255)
                     gACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAAGCGATGATATTGACCACGAGATTGCCGATTTGc                                                                                                                                                                        <  2:139961/248‑1 (MQ=255)
                     gACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCGCGCAGAATCTGATTATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCAGTACGATAATATCGAGCAGAAAATCGAGGATATTGACCACGAGATTGCCGATTTGc                                                                                                                                                                        <  2:139949/248‑1 (MQ=255)
                                 ggcggcCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCTCGAGATTGCCGATTTGCAGGCGAAACGTAccc                                                                                                                                                         <  1:192457/251‑1 (MQ=255)
                                                           aGAATCTGATGATGCCTGCCTATTTCCTGATGGTAGTGGCGGTGGTTGGTTTAATCCCCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAAtt                                                                                                                               <  1:60084/251‑1 (MQ=255)
                                                            gAATCTGTGGATGCCTGCCTATTCCCTGATGGTAGTGGCGGTGGTTGGGTGAATCCCCGGCGTAACAATGAAAGAGACGGAAAAG‑CGTCCGTGGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGAGAAGGAAATCCTTGTC‑GAGCATAACGATAATATCGAGCAGAAAACCTATGATATTGACCACGAGATTGCCGATTTGAAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATtg                                                                                                                              <  2:66/251‑1 (MQ=255)
                                                                                       gATGGTAGTGGCGGTGGTTGTTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAAGCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTg                                                                                                    <  1:190846/250‑1 (MQ=255)
                                                                                                                                                                   gACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGCATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTACGGTTCACTTAATCTCTGACGCGCCTACTCTCCTCCAGGt                        >  2:190786/1‑251 (MQ=255)
                                                                                                                                                                         ggcggcGTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACgg                  >  2:82703/1‑251 (MQ=255)
                                                                                                                                                                               gTCAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggcggcgg            >  1:6350/1‑245 (MQ=255)
                                                                                                                                                                                 cAGATATACAGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGt          >  2:232604/1‑251 (MQ=255)
                                                                                                                                                                                       aaCAGG‑AAGCAAAGGAAATTCT‑GTCCGAGCATTACGATAATATCGAGAAGAAAATCCATGAGATTGACCACGAGATGGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCGCCTCCAGGTTAACGGAGGAGAGTGCAATg    <  2:40644/250‑1 (MQ=255)
                                                                                                                                                                                          aGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGaa  >  1:140402/1‑250 (MQ=255)
                                                                                                                                                                                          aGG‑AAGCGAAGGAAATTCTTGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGaa  >  2:155667/1‑250 (MQ=255)
                                                                                                                                                                                                              |                                                                                                                                                                                                                                       
TTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAAT‑CGTCCGTTGAAAGGTGCGACACCGGCGGCGTCAGATATACAGG‑AAGCGAAGGAAATTCTCGTC‑GAGCATTACGATAATATCGAGCAGAAAATCGATGATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAA  >  NC_000913/4331686‑4332120

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