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%	184.7 / NA	65	T108T (ACC→ACT)	SAUSA300_RS02500	rRNA (cytidine‑2'‑O‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (40/25);  total (40/25)

TGCGCTCGTATCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACCGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGAAAAGAGA  >  NC_007793/524894‑525081                                                                                              |                                                                                               tGCGCTCGTATCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCt                                                                                          >  1:2625003/1‑100 (MQ=255) gctcttccgaTCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCt                                                                                          <  1:4407895/91‑1 (MQ=255)   tcttccgaTCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTgg                                                                                        <  1:995653/93‑1 (MQ=255)   cttccgatcttGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTgg                                                                                        <  1:3600281/90‑1 (MQ=255)     cTCGTATCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGAc                                                                                      >  1:1647346/1‑100 (MQ=255)       cGTATCTGATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCt                                                                                    >  1:3466918/1‑100 (MQ=255)              gATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTg                                                                             >  1:587348/1‑100 (MQ=255)              gATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTg                                                                             >  1:5582208/1‑100 (MQ=255)              gATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTg                                                                             >  1:69123/1‑100 (MQ=255)              gATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTg                                                                             <  1:4514747/100‑1 (MQ=255)              gATGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTg                                                                             <  1:2572909/100‑1 (MQ=255)               aTGGTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTga                                                                            <  1:2063982/100‑2 (MQ=255)               aTGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTgg                                                                            <  1:5693019/100‑1 (MQ=255)               aTGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTgg                                                                            <  1:4233523/100‑1 (MQ=255)                tGCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTggg                                                                           <  1:1835461/100‑1 (MQ=255)                 gCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGc                                                                          >  1:1733525/1‑100 (MQ=255)                 gCTGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGGGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTggat                                                                          <  1:3403000/100‑3 (MQ=255)                   tGGATTGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTg                                                                        <  1:5681436/100‑1 (MQ=255)                        tGCCCTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGc                                                                   >  1:4876613/1‑100 (MQ=255)                            cTTAATTAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTg                                                               <  1:2838263/100‑1 (MQ=255)                                 ccaGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCT‑GGCTGACGGCTTTGATGGCt                                                         >  1:5050682/3‑100 (MQ=255)                                  tAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCt                                                         <  1:3836938/100‑1 (MQ=255)                                  tAGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCt                                                         <  1:1760874/100‑1 (MQ=255)                                   aGTGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTa                                                        >  1:1116464/1‑100 (MQ=255)                                    gcgATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATggagtg                                                       <  1:4565872/98‑5 (MQ=255)                                     tGATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGt                                                      <  1:1486416/100‑1 (MQ=255)                                      gATCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTg                                                     >  1:5528711/1‑100 (MQ=255)                                      aaTCCTGGATACGAATTAGTAGTGGCAGCCAGTGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTg                                                     >  1:2596377/2‑100 (MQ=255)                                       aTCCTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTgg                                                    <  1:612273/100‑1 (MQ=255)                                          cTGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGAtt                                                 >  1:393027/1‑100 (MQ=255)                                           tGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGGCTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTa                                                >  1:3055516/1‑100 (MQ=255)                                           tGGATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTa                                                >  1:801165/1‑100 (MQ=255)                                            ggATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTAc                                               >  1:635899/1‑100 (MQ=255)                                            ggATACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTAc                                               >  1:4627600/1‑100 (MQ=255)                                              ctACGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCt                                             >  1:3113720/2‑100 (MQ=255)                                                aCGAATTAGTAGTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTc                                           >  1:955568/1‑100 (MQ=255)                                                      tagtagTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGt                                     <  1:390230/100‑1 (MQ=255)                                                      tagtagTGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGt                                     >  1:3571725/1‑100 (MQ=255)                                                        gtagtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGtat                                   <  1:203112/100‑1 (MQ=255)                                                         tagtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCacggccaat                                  <  1:4464641/100‑9 (MQ=255)                                                         tagtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGtatt                                  <  1:617211/100‑2 (MQ=255)                                                         tagtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGtata                                  >  1:4852053/1‑100 (MQ=255)                                                          cgtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGtatat                                 >  1:981948/2‑100 (MQ=255)                                                          agtGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGtatat                                 >  1:505119/1‑100 (MQ=255)                                                           ctGGCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGtatata                                >  1:4755825/2‑100 (MQ=255)                                                             ggCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACa                              >  1:5788592/1‑100 (MQ=255)                                                             ggCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACa                              >  1:3439299/1‑100 (MQ=255)                                                              gCAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACAt                             >  1:2479942/1‑100 (MQ=255)                                                               taGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACAtt                            >  1:4727310/2‑100 (MQ=255)                                                               cAGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACAtt                            >  1:2203678/1‑100 (MQ=255)                                                                aGCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCGTTGATGGCTAGTGGATTACCTTCATATGTATATACAttt                           >  1:5655889/1‑100 (MQ=255)                                                                 gCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACAtttt                          >  1:2399216/1‑100 (MQ=255)                                                                 gCCAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACAtttt                          >  1:1838213/1‑100 (MQ=255)                                                                  ccAGAGAAGCTAATATTAAAGTAGAGACTGTGCCTGGACCTAATGCTTGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACAttttt                         <  1:2655823/100‑1 (MQ=255)                                                                             aaTATTAAAGTAGCGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGc              >  1:2575352/1‑100 (MQ=255)                                                                              atatTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGcc             <  1:1911206/100‑1 (MQ=255)                                                                                atTAAAGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGGATATACATTTTTAGGATTTTTGCCAc           >  1:3935569/1‑100 (MQ=255)                                                                                  tAAAGTAGAGATTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACgt         <  1:3944224/100‑2 (MQ=255)                                                                                   aaaGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATTTACATTTTTAGGATTTTTGCCACGaa        >  1:3906108/1‑100 (MQ=255)                                                                                   aaaGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGaa        >  1:101607/1‑100 (MQ=255)                                                                                   aaaGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGaa        >  1:3220977/1‑100 (MQ=255)                                                                                   aaaGTAGAGACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGaa        >  1:3175479/1‑100 (MQ=255)                                                                                   aaaGTAGAGACTGTGCCTGGAACTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGaa        >  1:5419720/1‑100 (MQ=255)                                                                                        agagACTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGAAAagag   >  1:762484/1‑100 (MQ=255)                                                                                         gagaCTGTGCCTGGACCTAATGCTGGGCTGACGGCTTTGATGGCTAGTGGATTACCTTCATATGTATATACATTTTTAGGATTTTTGCCACGAAAAgaga  <  1:1232497/100‑1 (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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑