Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 339,986 G→A R73Q (CGA→CAA)  yahH → putative uncharacterized protein YahH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913339,9860GA83.9% 64.3 / 6.2 31R73Q (CGA→CAA) yahHputative uncharacterized protein YahH
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (7/19);  total (10/21)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

ACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGA‑CGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCG  >  NC_000913/339743‑340230
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
aCCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAGTGGGGCGAGAAttt                                                                                                                                                                                                                                                >  1:84422/1‑249 (MQ=255)
   tACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGGTAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑ACCCTTATCCGGCCTACGGGTGGCGCGGGAAttgt                                                                                                                                                                                                                                               >  1:262594/1‑246 (MQ=11)
    acagggtCATGCCGGAGGCGGCGTGAACGCCTTATCCGGCCTCCGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAGTGGCGCGAGAAttgta                                                                                                                                                                                                                                              <  2:340872/243‑5 (MQ=11)
            aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTCCGAGTGGCGCGGGAAT‑TGTATGCCTgcta                                                                                                                                                                                                                                    >  1:156402/1‑247 (MQ=11)
                         cgttgaCGCCTTATCCGGCCTACGAAGGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAGTGGCGCGAGAAT‑TGTAGGCCTGATAAGCGTAGCGCATCa                                                                                                                                                                                                                      <  2:188401/246‑1 (MQ=11)
                                                                                                                                           aCGAATGGCACAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCAt                                                                                                     <  2:123716/250‑1 (MQ=255)
                                                                                                                                            cGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATaa                                                                                                   <  2:20226/251‑1 (MQ=255)
                                                                                                                                               aTGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAAccc                                                                                                <  2:145310/251‑1 (MQ=255)
                                                                                                                                                  gcgcAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCAt                                                                                             <  2:112758/251‑1 (MQ=255)
                                                                                                                                                     cAAGAAGCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGTCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCAc                                                                                          <  2:305253/251‑1 (MQ=14)
                                                                                                                                                          aTCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGCGTGGCGCGAGGATTGTCGGCATGCTATGCGTAGCGCAA‑CAGGCAGTCTGGCGAGGGCGATAAACCCATCACCctct                                                                                     >  2:232117/1‑251 (MQ=11)
                                                                                                                                                              gTAGGCCTGATAAGCTTAGCGCATCCGGCATTGGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCctct                                                                                     <  2:228338/247‑1 (MQ=14)
                                                                                                                                                              gTAGGCCTGATAAGCGTAGCGCATCAGGCCTTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGt                                                                                   <  2:138224/249‑1 (MQ=255)
                                                                                                                                                              gTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCTGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATATCCCCATCACCctct                                                                                     <  1:330938/247‑1 (MQ=14)
                                                                                                                                                               tAGGCCTGATAAGCGTAGCGCAGCAGGCATTTGTCACCATTGCCGTATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGTCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAg                                                                                 <  2:357741/250‑1 (MQ=14)
                                                                                                                                                                               aGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATg                                                                 <  1:176242/250‑1 (MQ=255)
                                                                                                                                                                                               gtcaccatTGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGGCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGGCGt                                                >  1:65770/3‑251 (MQ=18)
                                                                                                                                                                                                   ccattGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTCCGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTg                                              <  1:16142/245‑1 (MQ=255)
                                                                                                                                                                                                   ccattGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCTCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCAGCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGcc                                            <  1:375058/247‑1 (MQ=255)
                                                                                                                                                                                                    cattGCCGGATGCTGC‑TGGAACGCCTTTGCCGGCCTACGAAGGGCGCAAGAAGCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGTATGCGGCGTGAACGCCTGATCCGGCCTCCGAGGGGCTCGAGAATTGTAGGCCTGATAAGCTTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATCACCCCTTCTCCCTCTGTAGCGGACATAACAACATGAAAGAGTTTGTGTTCGTTGCCa                                           <  2:268934/248‑1 (MQ=17)
                                                                                                                                                                                                    cattGCCGGATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGGGAGAATTGTAGGCCTGATAAGCGTCGCGCAT‑CCAGAAGTCTGGGGTTGGGCATAACCCCATCACCCCCTGTAGCGGGCATAACAAAATGAAAGAGCGTGTGGTCGtttccc                                           >  2:349973/4‑247 (MQ=17)
                                                                                                                                                                                                            gATGCGGCGTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAAAGCCTTATTCGGCCTACGAGGGGCGCGGGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGGCATAACAACATGAAAAAGGTTGTGGTAGTTGCCATTGGTGGc                                   >  1:174017/1‑251 (MQ=18)
                                                                                                                                                                                                                cggcgTGAA‑CGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAAAAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAAc                                >  1:171050/1‑250 (MQ=255)
                                                                                                                                                                                                                          cGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCCCCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCa                      <  1:47399/251‑1 (MQ=255)
                                                                                                                                                                                                                                 tCCGGCCTACGAATGGCGCAAGAATCTGTCGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTTCGAGTGGCGCGGGAATTGTAGGCCTGCTAAGCGTAGCGCATACAG‑CAGGCTGGCGTTGGTCATAACCCCATCACCCTCTGTCACGGGCCTAACAAAATGAAAGCCATTGTGGTCGTTGCCATTGTGGGGAACCGCATTTTCaaaagta                >  2:277469/1‑246 (MQ=17)
                                                                                                                                                                                                                                  ccGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAAc              <  1:349511/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ccGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAAc              <  1:270812/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ggCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgc            <  1:97862/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGGTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTCGGGTCGTTGCCATTTGGGGCAACAGCATTATCAAAGATAACgcc           >  2:44461/1‑249 (MQ=21)
                                                                                                                                                                                                                                          cgaatGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCGACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCGCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACgccagcca      <  1:239598/247‑1 (MQ=255)
                                                                                                                                                                                                                                               ggCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCg  <  2:156470/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
ACCTACATGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAATGGCGCAAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGCCACCACTGCCGGATGCGGCGTGGA‑CGCCTTATCCGGCCTACGAGTGGCGCGAGAATCTGTAGGCCTGATAAGCGTAGCGCATCAGGCATTTGTCACCATTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGAGTGGCGCGAGAATTGTAGGCCTGATAAGCGTAGCGCAT‑CAGGCAGTCTGGCGTTGGTCATAACCCCATCACCCTCTGTAGCGGACATAACAACATGAAAGAGCTTGTGGTCGTTGCCATTGGTGGCAACAGCATTATCAAAGATAACGCCAGCCAGTCG  >  NC_000913/339743‑340230

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