Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,212,973 Δ1 bp intergenic (+27/‑67) rrlE → / → rrfE 23S ribosomal RNA/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,212,9730T.95.8% 80.2 / ‑3.5 24intergenic (+27/‑67)rrlE/rrfE23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (10/13);  total (10/14)
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.96e-01

AGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATG  >  NC_000913/4212910‑4213222
                                                               |                                                                                                                                                                                                                                                         
ggCTAACCTGTACTAATGAACCGGGAGGCTTAACCTTACAACGCCGAAGATGTTTTGTCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCAGGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGAAGGGCACTGCCAGGCATc                                                                <  2:594966/250‑1 (MQ=11)
      ccGGTACTAATGTACCGTGAGGCTTAACGTTNCAACGCCGAAGGTGTGTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAAGGCAGAAGCGGTCTGATAAAACAGAAGTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAg                                                         <  2:354806/251‑1 (MQ=11)
      ccGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAg                                                         <  1:240296/251‑1 (MQ=11)
           aCTAATGAACCGTGAGGCTTACCCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTGGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGaaaa                                                     <  2:642566/250‑1 (MQ=12)
              aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGCAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGAANGGAACTGCCCGGCATCAACTTAGAAAAAccc                                                 >  2:351331/1‑251 (MQ=255)
                tGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCACACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCg                                               >  1:254028/1‑251 (MQ=21)
                  aaCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGt                                             <  1:387103/251‑1 (MQ=25)
                   aCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTc                                            >  1:565152/1‑251 (MQ=25)
                      gTGAGGCTTAAACTTACAACGCCGACGGCGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGATGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGCAGTGAAACGCCGTAGCACCGCTGGTAGTGCGGGGTCTCCCCATGCGAGAGTAGGGAACTGNCAGGCATCACATTAGAAAAACACCGGTcaat                                         >  2:349131/1‑248 (MQ=11)
                        gAGGCTTAACCTTACAACGCCGAAGGTGATTTGGCGGAT‑GAGAGAAGATTTTCCGCCCGATACCGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGACTGGCGGCAGTAGCGAGGAGGTCCCACCTGACCCCATGCCGAAATAAGAAGTGAAACGCCGTAGCGCCGATGGCAGTGTGGGGTCTCCCCCTGCGAGAGTAGGGAACTGCCAGGCATCAATATAGAAAAACCCCGGTCCATa                                        >  1:101682/1‑250 (MQ=11)
                            cTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGCGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGcc                                   >  2:248679/1‑251 (MQ=255)
                            cTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGcc                                   >  1:78700/1‑251 (MQ=255)
                              tAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgg                                 >  1:236514/1‑251 (MQ=255)
                               atccTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACCGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCGGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgg                                 <  2:298394/248‑1 (MQ=255)
                                aCCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCgggg                               <  1:702859/251‑1 (MQ=255)
                                    tACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGcc                                   <  1:138877/243‑1 (MQ=255)
                                    tACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGcc                                   >  2:138877/1‑243 (MQ=255)
                                    tACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGttt                            <  1:79182/250‑1 (MQ=255)
                                     aCAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAAATGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGttttt                          <  1:402978/251‑1 (MQ=255)
                                     aCAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAAGTGGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGtttt                           <  2:37818/250‑1 (MQ=255)
                                       aaCGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTg                        <  1:254823/251‑1 (MQ=255)
                                                ggTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAat               >  1:35759/1‑251 (MQ=255)
                                                        ggCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGACGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCTGTCCATAAGGCCGGGGTTTTTTGCATATCAATGAATTGCa       <  2:412167/251‑1 (MQ=255)
                                                         gCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACGCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTACTTGCat      <  2:120471/251‑1 (MQ=255)
                                                           cgtgtGAGAGAAGATTTTCAGCCTGATACAGATTAGATCAGAACGCAGAAGCGGTCTGATAAATCAGAATTTGCCTGGCGGCAGTAGCACGGTGGTCCCACCTGACGCCATGCCGAACGCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGGGGGGGCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAGCCCCGGTCCATAAGGCCGGGGTTGTTTGCATATCAATTATTTGCat      <  2:395039/246‑1 (MQ=17)
                                                              atGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCatgatg  >  1:469918/2‑251 (MQ=255)
                                                               |                                                                                                                                                                                                                                                         
AGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGGCCGGGGTTTTTTGCATATCAATTATTTGCATGATG  >  NC_000913/4212910‑4213222

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