Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,152 T→C 10.0% intergenic (+57/‑123) rnr → / → rlmB exoribonuclease R, RNase R/23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,1520TC10.0% 184.7 / 9.9 60intergenic (+57/‑123)rnr/rlmBexoribonuclease R, RNase R/23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base C (3/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

AAAAAGCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTA  >  NC_000913/4409015‑4409291
                                                                                                                                         |                                                                                                                                           
aaaaaGCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATgcgg                                                                                                                                  >  4:1065846/1‑149 (MQ=255)
aaaaaGCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATgcgg                                                                                                                                  >  2:8057/1‑149 (MQ=255)
aaaaaGCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATgcgg                                                                                                                                  >  3:435996/1‑149 (MQ=255)
aaaaaGCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATgcgg                                                                                                                                  >  3:1269891/1‑149 (MQ=255)
  aaaGCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATgcggcg                                                                                                                                >  4:334785/1‑149 (MQ=255)
   aaGCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGt                                                                                                                               <  3:334785/149‑1 (MQ=255)
      ccATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGaa                                                                                                                            <  1:864405/149‑1 (MQ=255)
       cATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATAAGGCGttcac                                                                                                                           >  1:120609/1‑145 (MQ=255)
         tCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCAt                                                                                                                                             >  2:125534/1‑129 (MQ=255)
         tCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCAt                                                                                                                                             <  1:125534/129‑1 (MQ=255)
           ggCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCt                                                                                                                       <  1:768971/149‑1 (MQ=255)
             cGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                                                                     <  1:48506/149‑1 (MQ=255)
                 aaCGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                 <  4:1269891/149‑1 (MQ=255)
                 aaCGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                 <  4:435996/149‑1 (MQ=255)
                  aCGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                <  3:1065846/149‑1 (MQ=255)
                                         cAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                         >  1:796869/1‑149 (MQ=255)
                                           aagcgaagcgTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCACGCCGCATCCGGCATGaa                                                                                                                            >  3:1188897/1‑112 (MQ=255)
                                           aagcgaagcgTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCACGCCGCATCCGGCATGaa                                                                                                                            <  4:1188897/112‑1 (MQ=255)
                                             gcgaagcgTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCACGCCGCATCCGGCATGa                                                                                                                             <  4:636379/109‑1 (MQ=255)
                                             gcgaagcgTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCACGCCGCATCCGGCATGa                                                                                                                             >  3:636379/1‑109 (MQ=255)
                                                 agcgTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGc                                                                                 >  1:12772/1‑149 (MQ=255)
                                                      gcggcgAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAgaga                                                                            <  2:12772/149‑1 (MQ=255)
                                                              gaaAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAg                                                                    >  4:1450091/1‑149 (MQ=255)
                                                                 aaaaGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc                                                                 <  3:1450091/149‑1 (MQ=255)
                                                                 aaaaGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc                                                                 <  1:8057/149‑1 (MQ=255)
                                                                 aaaaGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc                                                                 >  3:73225/1‑149 (MQ=255)
                                                                        gCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAt                                                                                             <  3:1057153/114‑1 (MQ=255)
                                                                        gCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAt                                                                                             >  4:1057153/1‑114 (MQ=255)
                                                                         cAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATaa                                                                                                       >  1:1119407/1‑103 (MQ=255)
                                                                         cAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATaa                                                                                                       <  2:1119407/103‑1 (MQ=255)
                                                                          agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGa                                                              >  4:765355/1‑143 (MQ=255)
                                                                          agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGa                                                              <  3:765355/143‑1 (MQ=255)
                                                                          agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCACGCCGCATCCGGCATGaa                                                                                                                            >  2:380477/1‑81 (MQ=38)
                                                                          agagTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCACGCCGCATCCGGCATGaa                                                                                                                            <  1:380477/81‑1 (MQ=38)
                                                                            agTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTa                                                      <  4:73225/149‑1 (MQ=255)
                                                                                aTCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                         >  2:1332368/1‑110 (MQ=255)
                                                                                aTCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                         <  1:1332368/110‑1 (MQ=255)
                                                                                 tCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                                                                                       <  4:201854/111‑1 (MQ=255)
                                                                                 tCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAata                                                                                       >  3:201854/1‑111 (MQ=255)
                                                                                   aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAgaga                                                                            <  4:592213/120‑1 (MQ=255)
                                                                                   aaTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAgaga                                                                            >  3:592213/1‑120 (MQ=255)
                                                                                    aTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                            <  2:392826/87‑1 (MQ=255)
                                                                                    aTACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                            >  1:392826/1‑87 (MQ=255)
                                                                                        ccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGTACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCa                                          >  4:848851/1‑149 (MQ=255)
                                                                                        ccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                         <  3:730074/102‑1 (MQ=255)
                                                                                        ccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                         >  4:730074/1‑102 (MQ=255)
                                                                                         ctctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATAtt                                                                                   >  4:737114/1‑107 (MQ=255)
                                                                                         ctcgTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATAtt                                                                                   <  3:737114/107‑1 (MQ=255)
                                                                                              tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCg                                                                                                                                         >  4:1216832/1‑48 (MQ=255)
                                                                                              tAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCg                                                                                                                                         <  3:1216832/48‑1 (MQ=255)
                                                                                                 aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAgg                                            >  3:599597/1‑138 (MQ=255)
                                                                                                 aagaagAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAgg                                            <  4:599597/138‑1 (MQ=255)
                                                                                                      gagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                                                                <  2:151406/97‑1 (MQ=255)
                                                                                                      gagGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                                                                >  1:151406/1‑97 (MQ=255)
                                                                                                             aGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt                                                                        <  2:1157195/98‑1 (MQ=255)
                                                                                                             aGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCAt                                                                        >  1:1157195/1‑98 (MQ=255)
                                                                                                                         aaaTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGaaa         >  2:582185/1‑149 (MQ=255)
                                                                                                                            tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                         >  1:823711/1‑66 (MQ=255)
                                                                                                                            tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCaa                                                                                         <  2:823711/66‑1 (MQ=255)
                                                                                                                                cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTa  <  3:848851/149‑1 (MQ=255)
                                                                                                                                        aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt                                                       >  4:702887/1‑88 (MQ=255)
                                                                                                                                        aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGt                                                       <  3:702887/88‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                           
AAAAAGCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTA  >  NC_000913/4409015‑4409291

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