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.T97.8% 144.8 / ‑4.1 46intergenic (+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 (22/23);  total (22/24)
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

GGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920072‑1920330
                                                                                                                               |                                                                                                                                      
ggCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑tttttt‑ttgcca                                                                                                                          >  2:1241806/1‑134 (MQ=255)
 gCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAg                                                                                                                          <  2:128094/139‑1 (MQ=255)
           ccTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTCTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAg                                                                                                                          <  2:1300775/129‑1 (MQ=255)
           ccTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTCTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAg                                                                                                                          >  1:1300775/1‑129 (MQ=255)
              gTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCg                                                                                                             <  2:1250646/139‑1 (MQ=255)
                                                   ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCg                                                                                                             <  2:763762/102‑1 (MQ=255)
                                                   ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCg                                                                                                             >  1:763762/1‑102 (MQ=255)
                                                        tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTAc                                                                                                                    <  1:992176/90‑1 (MQ=255)
                                                        tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGacac                                                                   >  2:441775/1‑139 (MQ=255)
                                                        tGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACGGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTAc                                                                                                                    >  2:992176/1‑90 (MQ=255)
                                                          cgcTTCCCAAATAATGCCCACTGCTCGGGCGTGCCTGCGCCGGAGTGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGcc                                                                                                        <  2:1487314/100‑1 (MQ=255)
                                                          cgcTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGcc                                                                                                        >  1:1487314/1‑100 (MQ=255)
                                                                         gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑tttttt‑ttg                                                                                                                             >  2:677093/1‑61 (MQ=255)
                                                                         gCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑tttttt‑ttg                                                                                                                             <  1:677093/63‑3 (MQ=255)
                                                                              cTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑TTTTTTGTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCttg                                             <  1:1428428/139‑2 (MQ=255)
                                                                                gCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                >  1:313889/1‑86 (MQ=255)
                                                                                gCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                <  2:313889/86‑1 (MQ=255)
                                                                                  tCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACGGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACg                                         <  1:235757/139‑1 (MQ=255)
                                                                                     ggCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   <  1:1075260/78‑1 (MQ=255)
                                                                                     ggCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   >  2:1075260/1‑78 (MQ=255)
                                                                                         tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAt                                  >  1:397599/1‑139 (MQ=255)
                                                                                         tGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAt                                  >  1:1392985/1‑139 (MQ=255)
                                                                                           ccTGCGCCGGAGTGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGt                                                                                                                         <  1:1028198/50‑1 (MQ=255)
                                                                                           ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGt                                                                                                                         >  2:1028198/1‑50 (MQ=255)
                                                                                           ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑tttttt‑ttgcc                                                                                                                           <  2:693579/47‑5 (MQ=255)
                                                                                           ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑tttttt‑ttgcc                                                                                                                           >  1:693579/1‑43 (MQ=255)
                                                                                           ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑tttttt‑ttgcc                                                                                                                           <  1:359950/47‑5 (MQ=255)
                                                                                           ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑tttttt‑ttgcc                                                                                                                           >  2:359950/1‑43 (MQ=255)
                                                                                              gcgcCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTtgat                             <  1:441775/139‑1 (MQ=255)
                                                                                                gcCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGAc                                                                                <  2:186922/86‑1 (MQ=255)
                                                                                                gcCGGAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGAc                                                                                >  1:186922/1‑86 (MQ=255)
                                                                                                   ggAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                <  1:1374815/67‑1 (MQ=255)
                                                                                                   ggAGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                >  2:1374815/1‑67 (MQ=255)
                                                                                                     aGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCg                                                                                                             <  1:1315738/52‑1 (MQ=255)
                                                                                                     aGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCg                                                                                                             >  2:1315738/1‑52 (MQ=255)
                                                                                                     aGCGTTTATGCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTc                      >  1:817100/1‑139 (MQ=255)
                                                                                                              gCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTAtt                                                                                        >  1:10367/1‑64 (MQ=255)
                                                                                                              gCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTAtt                                                                                        <  2:10367/64‑1 (MQ=255)
                                                                                                              gCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGCGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATg                                                  <  1:479217/102‑1 (MQ=255)
                                                                                                              gCTAAACTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGCGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATg                                                  >  2:479217/1‑102 (MQ=255)
                                                                                                                 aaaCTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCg                                                                                                             >  2:692569/1‑40 (MQ=38)
                                                                                                                 aaaCTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCg                                                                                                             <  1:692569/40‑1 (MQ=255)
                                                                                                                 aaaCTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGAt                                                   >  1:871275/1‑98 (MQ=255)
                                                                                                                 aaaCTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGAt                                                   <  2:871275/98‑1 (MQ=255)
                                                                                                                  aaCTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCg                                            >  1:1186366/1‑104 (MQ=255)
                                                                                                                  aaCTGCGCGCCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCg                                            <  2:1186366/104‑1 (MQ=255)
                                                                                                                       cgcgcCTGTTTTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGAtt                         >  1:805409/1‑119 (MQ=255)
                                                                                                                       cgcgcCTGTTTTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGAtt                         <  2:805409/119‑1 (MQ=255)
                                                                                                                       cgcgcCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                >  1:704660/1‑47 (MQ=39)
                                                                                                                       cgcgcCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAacacac                                                                                                <  2:704660/47‑1 (MQ=39)
                                                                                                                       cgcgcCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGa                                                                                 <  2:1153334/62‑1 (MQ=255)
                                                                                                                       cgcgcCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGa                                                                                 >  1:1153334/1‑62 (MQ=255)
                                                                                                                         cgcCTGT‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCg  <  2:424054/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                      
GGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑‑TTTTTT‑TGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920072‑1920330

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