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,3490CT97.1% 99.4 / ‑6.3 34E37E (GAG→GAAyaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (19/14);  minor base A (0/1);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGAT  >  NC_000913/6104‑6587
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
cgcaGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATg                                                                                                                                                                                                                                           >  1:216321/4‑251 (MQ=255)
    gTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGACTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCCGTTTCCGCGCTTCATGGATc                                                                                                                                                                                                                                       >  1:175522/1‑251 (MQ=255)
             aCAACCCGTAAAGCATCCGCAATTGCTGTGTGGCACAATCGAATTCTTCTTCGCTTAAGTTTTCGCCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGACTCGCCGGGCGGGCATTCTCCGGCGTAAAATCGGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgctg                                                                                                                                                                                                                               <  1:321392/250‑1 (MQ=255)
                  ccGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAAtt                                                                                                                                                                                                                         >  1:446879/1‑251 (MQ=255)
                    ggAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGt                                                                                                                                                                                                                       >  1:323322/1‑251 (MQ=255)
                      aaaGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCt                                                                                                                                                                                                                     >  1:119746/1‑251 (MQ=255)
                              cGCAAATGCTGTGGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCAAGAATCGCCGGGGGGGTATTCTCCGTCGTGAAATCTGTCTCCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCGACTGCTGGGAATTGTCTAACAGCTc                                                                                                                                                                                                             <  1:462138/251‑1 (MQ=255)
                                cAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGCTGCGCATCCGCGTGCTAATCTGCGGGGGCGTCGGTTTCCGCGCTTCATGGGTTAACTGATGGGCATTGTCTGACAGCTCCg                                                                                                                                                                                                           >  2:84389/1‑251 (MQ=255)
                                                                 gCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCACGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGGTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCttt                                                                                                                                                                          >  2:103536/1‑251 (MQ=255)
                                                                  cTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTACCAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTg                                                                                                                                                                         >  2:17798/1‑251 (MQ=255)
                                                                     agggttccgcccTGAAAGCGGGTGTAGACATCACCTTAAACCGCCAGAACCGCCGGGCGGGCATTCCCCGGCGGGAAATCGGGCTGCCAGTCATGAAAGCGAGGGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTa                                                                                                                                                                      <  2:93578/241‑1 (MQ=255)
                                                                             ggccggaaagcctgtggaaaaaTCACCTTTAAACGCCAGAATCGCCTGGGGGGCATCTCCCGGGGTGAAAGGGGGCTGCCATGAATGAAAGCGAGCGGCGGTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAAGTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCg                                                                                                                                                              <  2:175522/230‑1 (MQ=255)
                                                                                          gtgtaacaatccccTTTAAACGCCAGACCCGCCGGGGGGGTATTTCCCGGGGGGAAAGCTGGCGGCCAGTAATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGa                                                                                                                                                 <  2:58918/239‑1 (MQ=255)
                                                                                                    cactgtataaccctaGAATCGCTCGGCGGGCATTCCCCGGCGTGAAATCTGCCTCCAAGTCTTGAAAGCGCGGGGCGTTGACCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATc                                                                                                                                       <  2:387777/237‑1 (MQ=255)
                                                                                                       cTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGc                                                                                                                                    >  1:325209/1‑251 (MQ=255)
                                                                                                          tAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgctg                                                                                                                                                                                                                               >  2:31334/1‑157 (MQ=255)
                                                                                                          tAAACGCCAGAACCGCCGGGCGGGCATTCCCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAActgctg                                                                                                                                                                                                                               <  1:31334/157‑1 (MQ=255)
                                                                                                                               gggCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAAttt                                                                                                            >  1:48326/1‑251 (MQ=255)
                                                                                                                                 ggcatcctccGGGGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAGTCAGCATATCCAGTCCTTGCAGGAAATTTa                                                                                                           <  1:326153/243‑1 (MQ=255)
                                                                                                                                                                                 gTTGATACCCGCCGGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGGGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGACAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGAAAGTTGTGAtt                                                           >  2:376052/1‑251 (MQ=255)
                                                                                                                                                                                      gcgccgtcAGTTTGTCGCGTATGCGCATCAGCGTGTTAAACGGCGGAGGTGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAAACAGCAAATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTAc                                                         <  2:266782/241‑1 (MQ=255)
                                                                                                                                                                                      cccccGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTcc                                                      <  2:376440/249‑1 (MQ=255)
                                                                                                                                                                                                 ttGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATAGACTGCTGGGGATTGTCTAACAGCTCAGGCAGCGTATCGCGCGTGGTGGTCAACGGGCTTTGGTAATCACGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCAAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAAcatc                                           >  2:172576/1‑251 (MQ=255)
                                                                                                                                                                                                        tGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGCTAGTTGTGATTACTCCTGCGAAACATCATCCCAc                                    >  2:365076/1‑251 (MQ=255)
                                                                                                                                                                                                         gATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCa                                     <  1:103536/249‑1 (MQ=255)
                                                                                                                                                                                                         gATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCAcg                                   <  1:376052/251‑1 (MQ=255)
                                                                                                                                                                                                               cATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCg                             >  2:103847/1‑251 (MQ=255)
                                                                                                                                                                                                                aTCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGTTTACTCCTGCGAAACATCATCCCACGCGTCCgg                            >  1:135348/1‑251 (MQ=255)
                                                                                                                                                                                                                   agcattcTAATCCGCGGAGGCGTCAGTTTCCGCGCTGCAATGATTAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCTTGGGTGTCAACGGTCTTTGTTAATCAAGCGTTTTCGCATTTGAAATAAGAATCAGCATAGCCAGTCCTGGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGGGGTGATCGATAGTTGTGAGTACTCCTGCGAAACATCATCCCACGCCCCCGgaga                         <  2:300866/245‑1 (MQ=255)
                                                                                                                                                                                                                       tGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGc                     >  1:344521/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAAAGGGCTTTGGTAATCAAGCGTTTTAGCAGGTGAAATAAGAATCAGCATATCAAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCggc              >  1:452265/1‑247 (MQ=255)
                                                                                                                                                                                                                                       cGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCg     >  1:71278/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cGTCAGTTTCCGCGCATCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCg     <  1:196885/251‑1 (MQ=255)
                                                                                                                                                                                                                                          catgtgccGCGCTTAATGGAGCAACGGCTGGGAATTGTCTAACAGCGCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAAAGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGGCCGATATCCGGAt  <  2:28122/245‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              
CGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTTATGCCGACTTTAGCAAAAAATGAGAATGAGTTGATCGATAGTTGTGATTACTCCTGCGAAACATCATCCCACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGAT  >  NC_000913/6104‑6587

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: