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.T93.9% 97.3 / ‑0.1 33intergenic (+55/‑25)yebT/rsmFMCE domain protein/16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base T (17/14);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

CAGGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑T‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920066‑1920330
                                                                                                                                     |                                                                                                                                     
cAGGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑ttttt                                                                                                                                >  2:293971/1‑139 (MQ=255)
 aGGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑tttttt                                                                                                                               >  2:934215/1‑139 (MQ=255)
      ggCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑tttttttgcca                                                                                                                          <  1:919787/139‑6 (MQ=255)
       gCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAg                                                                                                                          <  1:472659/139‑1 (MQ=255)
            cacTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTa                                                                                                                     <  2:15234/139‑1 (MQ=255)
              cTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACa                                                                                                                   <  1:29801/139‑1 (MQ=255)
                 ccTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGc                                                                                                                >  2:326810/1‑139 (MQ=255)
                 ccTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGc                                                                                                                >  2:874089/1‑139 (MQ=255)
                 ccTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑tttttt                                                                                                                               <  1:13419/123‑1 (MQ=255)
                 ccTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑tttttt                                                                                                                               >  2:13419/1‑123 (MQ=255)
                      tGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTg                                                                                                           <  1:934215/139‑1 (MQ=255)
                           gAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCa                                                                                                      >  2:146761/1‑139 (MQ=255)
                              aGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAaca                                                                                                   <  1:146761/139‑1 (MQ=255)
                                  aaCCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAacaca                                                                                                 >  2:296255/1‑137 (MQ=255)
                                  aaCCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAacaca                                                                                                 <  1:296255/137‑1 (MQ=255)
                                                  tGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACGTGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                               >  1:663385/1‑139 (MQ=255)
                                                  tGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACg                                                                               <  2:663385/139‑1 (MQ=255)
                                                   gAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGc                                                                              <  1:25108/139‑1 (MQ=255)
                                                         ggAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑T‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTg                                                                       >  2:733976/1‑139 (MQ=255)
                                                                     cccAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATgcg                                                            >  1:962781/1‑139 (MQ=255)
                                                                                cccACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGc                                                 >  1:60999/1‑139 (MQ=255)
                                                                                   aCTGCTCCGGCGGGCCTTCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCtt                                              <  2:141304/139‑1 (MQ=255)
                                                                                        tCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑ttttttt                                                                                                                              >  2:100460/1‑52 (MQ=255)
                                                                                        tCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑ttttttt                                                                                                                              <  1:100460/53‑2 (MQ=255)
                                                                                          cGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACg                                         <  2:1043608/137‑1 (MQ=255)
                                                                                          cGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACg                                         >  1:1043608/1‑137 (MQ=255)
                                                                                          cGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGct                                       >  2:759066/1‑139 (MQ=255)
                                                                                          cGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGct                                       >  2:972939/1‑139 (MQ=255)
                                                                                           ggCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGctc                                      <  1:972939/139‑1 (MQ=255)
                                                                                                 ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCg                                            <  1:769371/127‑1 (MQ=255)
                                                                                                 ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCg                                            >  2:769371/1‑127 (MQ=255)
                                                                                                 ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAttt                                >  2:1001859/1‑139 (MQ=255)
                                                                                                 ccTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCATTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCAttt                                >  2:753521/1‑139 (MQ=255)
                                                                                                   tGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑ttttttttgc                                                                                                                           >  1:468815/1‑41 (MQ=255)
                                                                                                   tGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑t‑ttttttttgc                                                                                                                           <  2:468815/45‑5 (MQ=255)
                                                                                                       ccGGAGCGTTTATGCTAAACTGCGCGCCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTgatgat                          <  1:874089/139‑1 (MQ=255)
                                                                                                          gAGCGTTTATGCTAAACTGCGCGCCTG‑t‑tttttttgc                                                                                                                            >  2:976306/1‑34 (MQ=255)
                                                                                                          gAGCGTTTATGCTAAACTGCGCGCCTG‑t‑tttttttgc                                                                                                                            <  1:976306/37‑4 (MQ=255)
                                                                                                           aGCGTTTATGCTAAACTGCGCGCCTG‑TCTTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTc                      >  1:561900/1‑139 (MQ=255)
                                                                                                                             cgcgcCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTgccgcc              <  2:132859/129‑1 (MQ=255)
                                                                                                                             cgcgcCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTgccgcc              >  1:132859/1‑129 (MQ=255)
                                                                                                                             cgcgcCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGcc    >  1:232413/1‑139 (MQ=255)
                                                                                                                              gcgcCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGccc   >  1:535172/1‑139 (MQ=255)
                                                                                                                               cgcCTGTT‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCg  >  1:594492/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                     
CAGGAAGGCAAACACTTCCTGTTGCAGGAAAGTGAACCGAAAGAGTGGCGTGAATGGGGAACTGCGCTTCCCAAATAATGCCCACTGCTCCGGCGTGCCTGCGCCGGAGCGTTTATGCTAAACTGCGCGCCTG‑T‑TTTTTTGCCAGTGGTACATGCTCGTGGCCCAACACACCGTTTATTTCCCGGACGCCTTTCTGACACAAATGCGCGAAGCGATGCCTTCGACGCTCTCATTTGATGATTTTCTTGCCGCCTGTCAGCGCCCG  >  NC_000913/1920066‑1920330

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