Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 6,239 G→T A74E (GCG→GAG)  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,2390GT89.5% 32.1 / ‑2.3 19A74E (GCG→GAG) yaaADNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (11/6);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

GTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTT  >  NC_000913/5992‑6481
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
gTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCg                                                                                                                                                                                                                                                 >  1:382642/1‑251 (MQ=255)
               gatgatATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGATTCGGCCTGCAAGCCGGTGTAAACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCt                                                                                                                                                                                                                                  >  1:170415/1‑251 (MQ=255)
                         ccccAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGt                                                                                                                                                                                                                            <  1:272168/247‑1 (MQ=255)
                             agcaTTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGAGTAGACCTCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCCGTGATGaa                                                                                                                                                                                                                      >  2:55928/4‑249 (MQ=255)
                                       cAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTATCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGt                                                                                                                                                                                                          >  1:323078/1‑251 (MQ=255)
                                             tttCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATAc                                                                                                                                                                                                    >  1:110341/1‑251 (MQ=255)
                                                cccTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATAccc                                                                                                                                                                                                  <  1:57765/250‑1 (MQ=255)
                                                                                        tAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTAGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTAGGCCTGCAAGCCGGTGTAGACATCACCTTTAAAAGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAATCTGGCTGCCAGTCTTGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGAATCAGCGTGCTAATCTTCgg                                                                                                                                                         >  1:255995/1‑251 (MQ=255)
                                                                                                cATTAAATCGAGCGGGCGCAAGACGCCATACAAGCCGGAAAGCATACGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTCGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGAATTCGCCGGCGTGATATCTGGCTGCCAGTCATGCAAGCGAGCGGCGTTGTTACCCGCCAGTTTGTCGCTGCTGGGCATCAGCGTGCTCATCttaggcggcatac                                                                                                                                                 >  2:206061/1‑239 (MQ=255)
                                                                                                                          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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
GTTCAGCTTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAATCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGTCCTTGCAGGAAATTT  >  NC_000913/5992‑6481

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