Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 5,642 A→C intergenic (+112/+41) yaaX → / ← yaaA DUF2502 domain‑containing protein YaaX/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_0009135,6420AC92.9% 52.6 / ‑3.2 29intergenic (+112/+41)yaaX/yaaADUF2502 domain‑containing protein YaaX/DNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (15/11);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

GGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGG  >  NC_000913/5403‑5881
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
ggAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCATGGGCCAGGGAAACATCACCGCTAAATAACAAAAGCCGGGGAAAAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTGATTCTGCAATATATTTAAACTGCATGCTTTTGTAGGCCGGGTAAGGCgg                                                                                                                                                                                                                                      >  1:213523/1‑250 (MQ=255)
                        gCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACACCACCGCTAAAAGACAAATGCCGGGTAACAATCCGGCATTCAACGCCCGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGGATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTg                                                                                                                                                                                                              >  1:254156/1‑251 (MQ=39)
                                                           ccgcGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCCAGGGCCCGGCAAACACCAACACTAAAAGACCAAATCCCGGTAACAATCCCACAGTCAACGCCTGATGCGCCACCGGCGCGTCTTTTCAGGCCTACGTGACATCTGCAATATATTGAATCTTCATTCTTTTATAGGCCGGGTTAGGCGGTCCCGCCGCCACCCGCCGTGACTGCACACTTATCGCCGCTCGttacgtcttaacaccc                                                                                                                                                                            >  2:111620/1‑235 (MQ=255)
                                                                      aaaGAAAGCGCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCTTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGc                                                                                                                                                                <  1:214483/250‑1 (MQ=255)
                                                                                                  ggTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTACAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAAGAAGCTTCATCAAAGAAGTAAcc                                                                                                                                    >  1:233410/1‑251 (MQ=255)
                                                                                                        ggTCCAGGCAAACATCACCGCTAAATGACAAATGACGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCt                                                                                                                              >  1:337338/1‑251 (MQ=255)
                                                                                                          cccagcccaaCATCACCGCTAAATGTCAAATGCCCTGTCACAATCCGGCGTTCAGCGCCTGATGTGTCGCTGGCGCTCCTTATCAGTCCTTCGTGAATTCTGCACTATATCTAATCTGCATGCTTTTGTAGGCACGATAAGGCGTTCACTCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTATAAACACCAGTTCGCCATTGCTGGAGGAAGCTCCATCAAATAAGTAACCTTCTCTAt                                                                                                                            <  2:15459/243‑1 (MQ=255)
                                                                                                                         ccGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGtt                                                                                                             <  1:48656/251‑1 (MQ=255)
                                                                                                                                  gACAAATGCCGGGTAACAATCCTGCATTCAGCGCCTGATGCGCCTCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGt                                                                                                     <  1:43455/250‑1 (MQ=255)
                                                                                                                                      aaaTGCCGGGTAACAAGCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTgg                                                                                                 <  1:237885/251‑1 (MQ=255)
                                                                                                                                         tGCCGGGTAACAATACAGCCTTCAGCGCCTGATGAGACGCTGGCGCCTATTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGACGGATAAGGCGTTCACGCCGCATCCGGAATTGACTGCAAACTTAACACTGCTCGAAGCGATTAAACACTAGTTCGCCATTGCTGGAGGAAACTTCATCAAAAAAGGAAGCTTCGCTATTAAAACCAGTCAGTTGCTATgggtggtgc                                                                                               >  2:213057/1‑243 (MQ=255)
                                                                                                                                             gggTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCg                                                                                          <  1:250317/251‑1 (MQ=255)
                                                                                                                                                    aattctgCATCCAGCGCCTTAGGCGACGCTGGCGCTCTTTTCCAGGCCTACGTGAATTCTGCAATATATTGAACCTGCATGCTTTTGTAGGTCGTATAAGGCGTTCACGCCGCATCCTGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCTGCCGAGTTTCa                                                                                   <  2:347576/245‑1 (MQ=255)
                                                                                                                                                            atccAGCGCCTAATGCGGCGCTGGCCCTGTCTTACGGGCCTACTGGAATTCTGACAAATATTGAAACTGCATGCTTTTGTCTGCCGGATAAGGCGTTCACGCCGCTTCCGGCCTTGACTGCAAACTTAACGCTGCTCGTAGCGTTTACACACCAGGTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTTGCTATTAAAACCAGGCAGTTGCTCTGGTTTTGTCAGCCGATTTTCaataataaa                                                                           <  2:1491/248‑4 (MQ=255)
                                                                                                                                                                       gATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtcatc                                                                 <  1:191856/250‑1 (MQ=255)
                                                                                                                                                                              tgcgggCGCGACTTATCAGGCATACGTGAATTCTGAAATATATTGAATAGGCATGCTTTTGTAGGCCGGATAAGGCGGTCACGCCGCATCCGGCATTGACTGCAAACGTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGAAACCTTCGCTATTAAAACCAGTCAGTTGCGCTGGTTTGGTCa                                                                                              <  2:125261/210‑1 (MQ=255)
                                                                                                                                                                              cGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCa                                                                                              >  1:125261/1‑214 (MQ=255)
                                                                                                                                                                                    cgcgTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCGGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCttt                                                   >  1:162992/1‑251 (MQ=255)
                                                                                                                                                                                               cAGGCCTACGGGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAg                                         <  1:111620/251‑1 (MQ=255)
                                                                                                                                                                                                 aGGCCTACGTGAATTCTGCAATATATTGAATCTGCCTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGt                                           >  1:248633/1‑248 (MQ=255)
                                                                                                                                                                                                  agcatACGTGAAATCTGCACTATATTGCATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGAATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACACGTTCGCCATTGCTGGAGGAAGCTTCATAAAAGAAGAAACCTTAGATATTAAAACAAGTAAGTAGCACTGGATTAGTAAGCCGATTTTCCATAATACACCGACTCATCAAACCGCGGGCTTTCTAAGCATAGa                                        >  2:233712/5‑250 (MQ=255)
                                                                                                                                                                                                          tgaaTTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCtgattat                              >  1:299986/3‑248 (MQ=255)
                                                                                                                                                                                                                    cAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTATCGTAGAAGATGATTATCTTAAAtt                      >  1:296708/1‑251 (MQ=255)
                                                                                                                                                                                                                          atTGAATCTACATGCTTTGGTAGGACGGATACGGCGGTCGCGCCGCATCAGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACCCCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGATGGGTCAGCCGATTTTCAATAATAAAACGCCTCATCAGACCGCGTGCGTTATTAGCGTCGAAgatgattacatataatttgcagt                >  2:272168/1‑234 (MQ=255)
                                                                                                                                                                                                                            tGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGttc              >  1:255503/1‑251 (MQ=255)
                                                                                                                                                                                                                                tCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTaa                                                                                                                         <  1:334357/140‑1 (MQ=255)
                                                                                                                                                                                                                                tCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACATTCGCTATTaa                                                                                                                         >  2:334357/1‑140 (MQ=255)
                                                                                                                                                                                                                                        cTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCGAgg  >  1:204648/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ttttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGCCTGCAAACTTAACGCTGCTCGTAGCGATTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCtt                                                    >  1:328368/1‑200 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
GGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCGAGG  >  NC_000913/5403‑5881

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