Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,349 C→T E37E (GAG→GAA yaaA ← DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009136,3490CT92.3% 68.6 / ‑3.6 26E37E (GAG→GAAyaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (13/11);  total (15/11)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGAT  >  NC_000913/6114‑6579
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
cATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCAACATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCTGCGTGAAATCTGGCTGACAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCTTCAGTTTCCGCGCCTCATGCATCACctgc                                                                                                                                                                                                                         >  1:98650/1‑251 (MQ=255)
                      aaatttctggtggacAAAATCGAAATCGTCTTCGATAAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACGTTTAACCGCCAGAATCGCCTGGCTGGCATTCTCCGGGGTGAAATCTGGGTGCCAGTCATGAAAGCGCGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGGCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATGGTCTAACAGCTCCg                                                                                                                                                                                                   <  2:88161/237‑1 (MQ=255)
                                     aaaaTCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGCTACCCGCCAGATTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCCGCGTATAGcgcg                                                                                                                                                                                    >  1:60617/1‑251 (MQ=255)
                                       aatcgaaatcgTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAACGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCCTCAGCGTGCTAAACTGCGGAGGCGTCAGTTTCCGCGCTTCATGGCTCAACTGCTGGGAATTGTCTAACAGCTCCGGCCGCGTATAGagcggg                                                                                                                                                                                  >  2:225063/1‑249 (MQ=255)
                                                          gAAGGTTTCGGCCTGCAAGCCGGTGTAGAAATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGATGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGAGTCAGTCTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCACAGGGCTTTGGt                                                                                                                                                               >  1:268704/1‑251 (MQ=255)
                                                                  cGGCCTGCAAGCCGGTGTAGACATCACCTTTAACCGCCAGAACCGCTGGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCTAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAg                                                                                                                                                        <  1:23691/250‑1 (MQ=255)
                                                                  cGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGGGGTCAACGGGCTTTGGTAATCAAGc                                                                                                                                                       >  1:146304/1‑251 (MQ=255)
                                                                                           aCCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGAGGCGTTGATACCAGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCAGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATCGCCCGTGGTGGTAAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGCATcg                                                                                                                              >  1:234557/1‑250 (MQ=255)
                                                                                                      ccAGAATCGCCTGGCGGTCATTCTCCGGCGTAAACTCTGGTGGCCAGTAATGAAAGCGAGCGGGGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGt                                                                                                                   <  1:349395/251‑1 (MQ=255)
                                                                                                                     gggCATTCTCCGGCGTGAATGCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAAttt                                                                                                    <  1:108975/251‑1 (MQ=255)
                                                                                                                                  gggaaaaTCGGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAg                                                                                       <  1:55928/247‑1 (MQ=255)
                                                                                                                                          cTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACTCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCGGCATATCCAGTCCTTGCAGGAGATTTATGCCGACTTTAGCAAAAaat                                                                               >  1:101835/1‑251 (MQ=255)
                                                                                                                                                gCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAaaagagaat                                                                         >  1:150608/1‑251 (MQ=255)
                                                                                                                                                     tCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGtt                                                                    >  1:180184/1‑251 (MQ=255)
                                                                                                                                                                    gcgTTGATACCAGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGCGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAACTTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTtgtgt                                                     >  2:21978/1‑250 (MQ=255)
                                                                                                                                                                     cgTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAGGCGTTTTCGCAGGTGAAATACGAATCAGAATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTCGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGAt                                                    >  1:361464/1‑251 (MQ=255)
                                                                                                                                                                        tGATACCCGCCAGTTTGTCGCTGCTGCGCATCAGCGTGCTCATCTGCGGCGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGGGATCAACGGGCTTTGGTATTCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAACACAAGAGAATGCGTTGATCGATAGTTGTGATTAc                                                 >  2:77666/1‑251 (MQ=255)
                                                                                                                                                                                    gTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAca                                     <  1:311906/251‑1 (MQ=255)
                                                                                                                                                                                    gTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAc                                      >  1:67962/1‑250 (MQ=255)
                                                                                                                                                                                        gTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAcatca                                  <  1:77666/250‑1 (MQ=255)
                                                                                                                                                                                             tGATGCGCATCAGCGTGCTAATCTGCGGAGTCGTCAGTTTCCGCGCTTCATGGAGCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGt                                                                                                                                                     <  2:143593/130‑1 (MQ=255)
                                                                                                                                                                                             tGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGt                                                                                                                                                     >  1:143593/1‑130 (MQ=255)
                                                                                                                                                                                                   gcATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTcc                      <  1:225063/251‑1 (MQ=255)
                                                                                                                                                                                                       cAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGga                   <  1:100022/250‑1 (MQ=255)
                                                                                                                                                                                                              cTAATCTGCGGAGGCGTCCGTTTCCGCGCTTCATGGATCACCTGCTGGGAATTGTCTAACCGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCCGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGACCGGAGAAAGCTg           >  2:144405/1‑251 (MQ=255)
                                                                                                                                                                                                                       ggAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGGCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGat  <  2:110341/251‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      
CATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGAT  >  NC_000913/6114‑6579

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