Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,913,699	C→T	*434* (TAG→TAA)	rlmD ←	23S rRNA m(5)U1939 methyltransferase, SAM‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,913,699	0	C	T	96.8%	186.5 / ‑2.0	62	*434* (TAG→TAA)	rlmD	23S rRNA m(5)U1939 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (29/31);  total (30/32)
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.83e-01

GTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTT  >  NC_000913/2913572‑2913833                                                                                                                                |                                                                                                                                       gTAATACCCAGA‑TTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAAcgcgc                                                                                                                            >  1:781955/1‑139 (MQ=255)         ccagattGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaaaa                                                                                                                     <  2:771417/134‑1 (MQ=255)             attGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTAc                                                                                                                 >  2:721323/2‑139 (MQ=255)             attGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTAc                                                                                                                 >  1:543271/2‑139 (MQ=255)               tGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCa                                                                                                               <  2:670280/139‑1 (MQ=255)                gCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCAt                                                                                                              <  2:451970/139‑1 (MQ=255)                    tCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGAt                                                                                                          <  2:196905/139‑1 (MQ=255)                                 aTCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                                                                             <  2:417774/139‑1 (MQ=255)                                           ccAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGaa                                                                                   >  1:118878/1‑139 (MQ=255)                                                  tATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATcc                                                                            >  1:110662/1‑139 (MQ=255)                                                         atGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCg                                                                     >  2:568233/1‑139 (MQ=255)                                                            tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGCAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                                  >  1:254507/1‑139 (MQ=255)                                                            tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                                  >  1:352308/1‑139 (MQ=255)                                                            tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                                  >  1:315387/1‑139 (MQ=255)                                                            tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                                  >  2:535506/1‑139 (MQ=255)                                                            tgCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCa                                                                  >  2:4305/1‑139 (MQ=255)                                                              cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaaaa                                                                                                                     <  2:84346/86‑1 (MQ=255)                                                              cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGaaaa                                                                                                                     >  1:84346/1‑86 (MQ=255)                                                              cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt                                                                <  2:799743/139‑1 (MQ=255)                                                              cACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGt                                                                >  2:397694/1‑139 (MQ=255)                                                                 ttcttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgc                                                             >  1:242020/1‑139 (MQ=255)                                                                   cttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                                                              >  1:391358/1‑104 (MQ=255)                                                                   cttACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                                                              <  2:391358/104‑1 (MQ=255)                                                                    ttACCGCAACCATCGTCCTCTCCTTTAGGTACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCa                                                          <  2:543271/139‑1 (MQ=255)                                                                      aCCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                                                              <  1:688891/101‑1 (MQ=255)                                                                      aCCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTccc                                                                                              >  2:688891/1‑101 (MQ=255)                                                                           aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCgag                                                                                                                         >  2:193519/1‑69 (MQ=255)                                                                           aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCgag                                                                                                                         <  1:193519/69‑1 (MQ=255)                                                                           aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                            <  2:67651/130‑1 (MQ=255)                                                                           aaCCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTcgcg                                                            >  1:67651/1‑130 (MQ=255)                                                                              cATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATc                                                <  2:352308/139‑1 (MQ=255)                                                                                       ctcCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTta                                       >  2:618454/1‑139 (MQ=255)                                                                                         ccTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTtaat                                     <  1:543354/139‑1 (MQ=255)                                                                                         ccTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTtaat                                     <  1:678128/139‑1 (MQ=255)                                                                                           tttAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTtaataa                                   >  1:132210/1‑139 (MQ=255)                                                                                             tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTaa                                                                                                                                 <  1:763997/43‑1 (MQ=255)                                                                                             tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTaa                                                                                                                                 >  2:763997/1‑43 (MQ=255)                                                                                             tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                                 <  1:397694/139‑1 (MQ=255)                                                                                             tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                                 <  2:695009/139‑1 (MQ=255)                                                                                             tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                                 <  2:132210/139‑1 (MQ=255)                                                                                             tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACg                                 <  2:242020/139‑1 (MQ=255)                                                                                             tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTggg                                                                                     >  2:589964/1‑87 (MQ=255)                                                                                             tAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTggg                                                                                     <  1:589964/87‑1 (MQ=255)                                                                                                  aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCg                            <  2:118878/139‑1 (MQ=255)                                                                                                  aCCAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCg                            >  1:503667/1‑139 (MQ=255)                                                                                                    cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                                                                             <  1:696106/72‑1 (MQ=255)                                                                                                    cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCg                                                                                             >  2:696106/1‑72 (MQ=255)                                                                                                    cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCt                                           <  1:76869/122‑1 (MQ=255)                                                                                                    cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCt                                           >  2:76869/1‑122 (MQ=255)                                                                                                    cAGACCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCt                          >  1:653794/1‑139 (MQ=255)                                                                                                       aCCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATg                                                       <  1:213699/107‑1 (MQ=255)                                                                                                       aCCTGCCGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATg                                                       >  2:213699/1‑107 (MQ=255)                                                                                                            ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGa                   <  2:244381/138‑1 (MQ=255)                                                                                                            ccGAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGa                   >  1:244381/1‑138 (MQ=255)                                                                                                              gAAATCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGcc                >  1:649202/1‑139 (MQ=255)                                                                                                                aaTCGGCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAg              <  2:254507/139‑1 (MQ=255)                                                                                                                     gCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTg         <  2:649202/139‑1 (MQ=255)                                                                                                                     gCAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTg         <  2:471997/139‑1 (MQ=255)                                                                                                                      cAAATCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGc        >  2:310395/1‑139 (MQ=255)                                                                                                                       aaaTCGCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCa       <  2:506315/139‑1 (MQ=255)                                                                                                                            gCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGtt  <  1:618454/139‑1 (MQ=255)                                                                                                                            gCAATTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGtt  <  1:721323/139‑1 (MQ=255)                                                                                                                                |                                                                                                                                       GTAATACCCAGACTTGCGATCCATTTTTCCGGATCAAATTCACCAGCCTTATTGATATGTGCACTTCTTACCGCAACCATCGTCCTCTCCTTTAGGGACCAGACCTGCCGAAATCGGCAAATCGCAACTATTTAACGCGCGAGAAAAGTACCATCGATTCCAGATGTCCCGTGTGTGGGAACATATCCAGCATCGCCAGTCGCGCAATGGTATATCCTGCTTTTAATAACGCTTCGCTATCCCGAGCCAGCGTTGCAGGGTT  >  NC_000913/2913572‑2913833

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