Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,811,179 G→A intergenic (‑27/+71) 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,1790GA96.9% 99.0 / ‑6.1 32intergenic (‑27/+71)mutM/rpmGDNA‑formamidopyrimidine glycosylase/50S ribosomal subunit protein L33
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/15);  minor base C (0/1);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

AATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTT  >  NC_000913/3810939‑3811422
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
aatgactgATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCCCTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCACCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATggg                                                                                                                                                                                                                                           >  2:194036/1‑251 (MQ=255)
        aTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGACCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCACCGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGc                                                                                                                                                                                                                                    >  1:82698/1‑250 (MQ=255)
               aGCCCTCAGGCAGCTCCAGCAGAAGATATTGAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACa                                                                                                                                                                                                                            <  2:274822/251‑1 (MQ=255)
                ccccTCAGGCATCCCCATCAGTATATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCGCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTTGTTTCATCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACaa                                                                                                                                                                                                                           <  1:179839/250‑1 (MQ=255)
                          aGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCccgc                                                                                                                                                                                                                 <  1:166682/251‑1 (MQ=255)
                                     cagctatatAGCCCGCCGCTGCACGCTAAGCACTGGTTGCTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCCGCGCAAGCGTCCGTTGCGCCCCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCACCTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGa                                                                                                                                                                                                                                 >  2:181617/9‑224 (MQ=255)
                                     cagATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCGTCTGAAACCGGCCAGCGCAAGCGGCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAGCGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGa                                                                                                                                                                                                                                 <  1:181617/224‑1 (MQ=255)
                                     cagATATTTAGCCCGCAGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGCGCGAGGt                                                                                                                                                                                                      >  1:320817/1‑251 (MQ=255)
                                      agATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtt                                                                                                                                                                                                     <  1:257528/251‑1 (MQ=255)
                                      agATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACCGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGACGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGtt                                                                                                                                                                                                     <  2:108340/251‑1 (MQ=255)
                                                                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)
                                                                                                                                                                                                                                         aGATGCAATCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTCGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGattt  <  1:181021/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
AATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCACCAACGAGATGCGGTTCTATGCCGCGGCGGCTGGTTTCAACTTCGGGTAATTCAGGCATAGCATCTCCAGGAATGAACAGATGCAGTCAATATGGGGGCAAGCAGATAACAAAAAACCCCGCCGGAGCGAGGTTTTTTGTTACATCAAAGCGAGAATTATTTGATTTTCGCTTCTTTGTAGATCACGTGCTGGCGAACAACTGGATCGAATTTTTTCAGTTCCAGTTTTTCCGGCTTAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTT  >  NC_000913/3810939‑3811422

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