Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	524,987	C→T	T108T (ACC→ACT)	SAUSA300_RS02500 →	rRNA (cytidine‑2'‑O‑)‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	524,987	0	C	T	100.0%	94.0 / NA	35	T108T (ACC→ACT)	SAUSA300_RS02500	rRNA (cytidine‑2'‑O‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/10);  total (25/10)

TGCGCTCGTATCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACCGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGAAAAGAGA  >  NC_007793/524894‑525081                                                                                              |                                                                                               gctcttccgaTCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCGAATATTAAAGTAGAGACTGTGCCt                                                                                          <  1:6171535/91‑1 (MQ=255)  ctcttccgaTCTGATGTTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCTAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTg                                                                                         <  1:5385707/92‑1 (MQ=255)   cgcTCGTATCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTgg                                                                                        >  1:5725770/1‑100 (MQ=255)   cgcTCGTATCTGATGCTAGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTgg                                                                                        <  1:5273863/100‑1 (MQ=255)    cttccgaTCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTggt                                                                                       <  1:4776613/94‑2 (MQ=255)            ctgctgctgGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGc                                                                               >  1:5735964/1‑100 (MQ=255)              gatgatGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTg                                                                             >  1:6625529/1‑100 (MQ=255)              gATGTTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTg                                                                             <  1:1664679/100‑1 (MQ=255)              gATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTg                                                                             >  1:6152180/1‑100 (MQ=255)               gtGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTgg                                                                            <  1:1623857/99‑1 (MQ=255)               gggCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTgg                                                                            <  1:2031360/98‑1 (MQ=255)                   tGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTg                                                                        <  1:6585435/100‑1 (MQ=255)                                     tGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGt                                                      <  1:2889564/100‑1 (MQ=255)                                        taccggATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGa                                                   >  1:611900/5‑100 (MQ=255)                                         ccTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGAt                                                  >  1:2257769/1‑100 (MQ=255)                                           tGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGGGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTa                                                >  1:2087909/1‑100 (MQ=255)                                           tGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTa                                                >  1:2042345/1‑100 (MQ=255)                                                aCGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTc                                           >  1:3783787/1‑100 (MQ=255)                                                     ttAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGAATGTGCCTGGAACTAATGCTGGGCTGACGGCTTTGTTGGCTAGAGGATTACCTTCATATg                                      >  1:5141883/1‑100 (MQ=255)                                                      tagtagTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGGCCTAATGCTGGGCTGACGGCTTTGATGGCTAGGGGATTACCTTCATATGt                                     >  1:5816347/1‑100 (MQ=255)                                                         tagtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGtatc                                  <  1:5960216/100‑2 (MQ=255)                                                         tagtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCAAGTGGATTACCTTCATATGtata                                  >  1:3496571/1‑100 (MQ=255)                                                         cagtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGtata                                  >  1:261032/2‑100 (MQ=255)                                                          agtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCCTATGtatat                                 >  1:5920191/1‑100 (MQ=255)                                                             ggCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACa                              >  1:4571200/1‑100 (MQ=255)                                                             ggCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACa                              >  1:4356408/1‑100 (MQ=255)                                                             cacAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACa                              >  1:1134304/3‑100 (MQ=255)                                                              gCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACAt                             >  1:3130736/1‑100 (MQ=255)                                                               cAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACAtt                            >  1:2664719/1‑100 (MQ=255)                                                                       cgcccTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATACACATTTTTAGGAt                    >  1:3891344/5‑100 (MQ=255)                                                                        aaGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGAtt                   >  1:406224/1‑100 (MQ=255)                                                                             aaTATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGc              >  1:2879906/1‑100 (MQ=255)                                                                                   aaaGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGaa        >  1:4757174/1‑100 (MQ=255)                                                                                       tAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTATGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGAAAtga    >  1:2268439/1‑100 (MQ=255)                                                                                         gagaCTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGAAAAgaga  >  1:6383104/1‑100 (MQ=255)                                                                                              |                                                                                               TGCGCTCGTATCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACCGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGAAAAGAGA  >  NC_007793/524894‑525081

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