Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,811,231 T→C intergenic (‑79/+19) mutM ← / ← rpmG DNA‑formamidopyrimidine glycosylase/50S ribosomal subunit protein L33

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,811,2310TC96.8% 104.2 / ‑6.3 31intergenic (‑79/+19)mutM/rpmGDNA‑formamidopyrimidine glycosylase/50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (14/16);  minor base G (0/1);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCA  >  NC_000913/3811003‑3811472
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
aagctcttgtttgtctcTTAAACGGTAGTTCTCTGTTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCCCCTCCGAATGAAGAATGGTTGCACCACCGAGATGCGTCTCTATGCCGCGACGGATGTTGTCAACTTCTGGTAATCCAGGCATATCATCCCCAGGAATGAACAGATTCACTCAATATGGGTGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAattat                                                                                                                                                                                                                             <  2:421395/235‑1 (MQ=255)
   cACTGGTTGGTCGCTTAATAGGTAGCTCTCTTCTGATACCGGCCCGCGCAAGCCTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCCTAGAATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTg                                                                                                                                                                                                                          >  2:357638/1‑251 (MQ=255)
        gttggtCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGAtttt                                                                                                                                                                                                                     >  1:300324/1‑251 (MQ=255)
                  tAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCtttt                                                                                                                                                                                                           >  2:3013/1‑250 (MQ=255)
                        gTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATGTTCGCTTCTTTGTAGATc                                                                                                                                                                                                     <  1:112050/251‑1 (MQ=255)
                          agatcgcgtctGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCTTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCAc                                                                                                                                                                                                   <  1:422956/243‑1 (MQ=255)
                                          cGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACAAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCCACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACt                                                                                                                                                                                   >  1:187623/1‑251 (MQ=255)
                                                        gTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAAtttttt                                                                                                                                                                     >  1:38519/1‑251 (MQ=255)
                                                            gTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCCGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGt                                                                                                                                                                 >  2:142288/1‑251 (MQ=255)
                                                                                          aCCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGttc                                                                                                                                   <  1:184659/251‑1 (MQ=255)
                                                                                                                       ggcTGGTTTCATCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTAcc                                                                                                      <  2:140306/251‑1 (MQ=255)
                                                                                                                                  aCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAgaagaa                                                                                            >  2:111806/1‑250 (MQ=255)
                                                                                                                                               cAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGAtttt                                                                              <  2:46892/251‑1 (MQ=255)
                                                                                                                                                   cATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTTTAGATCACGTGCTGGCGAACAACTGGATc                                                                                                                                                                              <  2:201566/151‑1 (MQ=255)
                                                                                                                                                   cATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATc                                                                                                                                                                              >  1:201566/1‑151 (MQ=255)
                                                                                                                                                    aTAGCATCTCCAGGAATGAACAGATGCCATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACcag                                                                                                    >  2:283496/1‑224 (MQ=255)
                                                                                                                                                    aTAGCATCGCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTGCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACcag                                                                                                    <  1:283496/224‑1 (MQ=255)
                                                                                                                                                        cATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAAt                                                                     <  1:350193/251‑1 (MQ=255)
                                                                                                                                                         aTCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAAt                                                                     >  1:326162/1‑250 (MQ=255)
                                                                                                                                                           ctcCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATAcc                                                                  >  1:177287/1‑251 (MQ=255)
                                                                                                                                                            tcCAGGACTGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAAACCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTATCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCt                                                                 >  1:88431/1‑251 (MQ=255)
                                                                                                                                                                   aagaaaagatgccatCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAACGCGAGAATTATTTGATTTTCGATTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGACGAAACCAGCTTGATTTTCTCACGAATACCTTTAGcc                                                           >  2:201838/15‑250 (MQ=255)
                                                                                                                                                                         aGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGattt                                                    <  1:181021/251‑1 (MQ=255)
                                                                                                                                                                             gcaatCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGATACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGatttattt                                                <  1:205688/247‑1 (MQ=255)
                                                                                                                                                                             gcaatCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGatttatt                                                 <  2:384101/246‑1 (MQ=255)
                                                                                                                                                                              caatCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTc                                               <  1:3013/248‑1 (MQ=255)
                                                                                                                                                                                          gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAagtactt                                   <  1:267725/251‑1 (MQ=255)
                                                                                                                                                                                          gggCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGGTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCATCGTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAagtactt                                   <  1:410146/251‑1 (MQ=255)
                                                                                                                                                                                                         aaaaaaCCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCAc                    <  1:401279/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggAGGTTTTCGGTTACATCAAGGCTGGATTTATTGGATTTTCGCTTCTTTGTGAATCAGGTGCTGGCAAACAACTGGATCGATTTTTTTCAGTTCCGGTTTTTCCGGCTAATTACTTTTTTTCTTCGTAGTGGTAAAGAGTTGACCGGTCCCAGCGGAAGAAACCAGCTTGATTTTCTCACGAGTACCTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGGTCAGCCa  <  2:38519/250‑1 (MQ=255)
                                                                                                                                                                                                                             ggTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGACTACCTTTAGCCATGATTTATTTCCTCTACGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGcc   >  2:350538/1‑248 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
AAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCA  >  NC_000913/3811003‑3811472

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