Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,920,205 (T)7→8 intergenic (+62/‑18) yebT → / → rsmF MCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,920,1981.T95.5% 64.0 / ‑3.1 22intergenic (+55/‑25)yebT/rsmFMCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (10/11);  total (10/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 = 9.33e-01

GGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920068‑1920330
                                                                                                                                   |                                                                                                                                     
ggAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAGTGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGtttttttt                                                                                                                              >  1:464235/1‑138 (MQ=255)
ggAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGAGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGtttttttt                                                                                                                              <  2:93525/139‑2 (MQ=255)
     gCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  1:269490/139‑1 (MQ=255)
            cTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACa                                                                                                                   <  2:582087/139‑1 (MQ=255)
               ccTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGc                                                                                                                >  1:638669/1‑139 (MQ=255)
                cTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCCTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTATTTTTGCCAGTGGTACATGCt                                                                                                               <  2:464235/139‑1 (MQ=255)
                         gAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCa                                                                                                      >  1:392745/1‑139 (MQ=255)
                             gTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacac                                                                                                  <  1:67272/139‑1 (MQ=255)
                               gAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                <  2:674709/139‑1 (MQ=255)
                                 aCCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCg                                                                                              <  2:566300/139‑1 (MQ=255)
                                                         aaCTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGa                                                                      <  2:392745/139‑1 (MQ=255)
                                                             gcgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGAcaca                                                                  >  1:284419/1‑139 (MQ=255)
                                                              cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgcca                                                                                                                          <  1:65718/81‑6 (MQ=255)
                                                              cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgcca                                                                                                                          >  2:65718/1‑76 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          >  1:216018/1‑67 (MQ=255)
                                                                             gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAg                                                                                                                          <  2:216018/67‑1 (MQ=255)
                                                                                   tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGGGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             >  1:469964/1‑74 (MQ=255)
                                                                                   tGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCg                                                                                                             <  2:469964/74‑1 (MQ=255)
                                                                                        cGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGct                                       >  2:179998/1‑139 (MQ=255)
                                                                                          gCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgcc                                                                                                                           >  2:19325/1‑48 (MQ=255)
                                                                                          gCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGttttttttgcc                                                                                                                           <  1:19325/52‑5 (MQ=255)
                                                                                                   cgcCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAac                                                                                                    >  2:194134/1‑67 (MQ=255)
                                                                                                   cgcCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAac                                                                                                    <  1:194134/67‑1 (MQ=255)
                                                                                                         aGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTc                      >  2:724028/1‑139 (MQ=255)
                                                                                                                      aaCTGCGCGCCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCa         <  1:179998/139‑1 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   >  1:703024/1‑44 (MQ=255)
                                                                                                                           cgcgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   <  2:703024/44‑1 (MQ=255)
                                                                                                                             cgcCTGTTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCg  >  1:328756/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
GGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920068‑1920330

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