Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,215,913	T→C	Y306H (TAT→CAT)	SAUSA300_RS06010 →	16S rRNA (cytosine(967)‑C(5))‑methyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,215,913	0	T	C	100.0%	142.2 / NA	53	Y306H (TAT→CAT)	SAUSA300_RS06010	16S rRNA (cytosine(967)‑C(5))‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (50/3);  total (50/3)

AAGGGCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAATATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGGTGTAATGAGACATAAGCCGG  >  NC_007793/1215815‑1216003                                                                                                   |                                                                                           taGGGCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATAT‑AAAAAATTACGATTAACAAATATTAAAGCTTTTCAacat                                                                                          >  1:4957009/2‑100 (MQ=255) caGGGCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATAT‑AAAAAATTACGATTAACAAATATTAAAGCTTTTCAacat                                                                                          >  1:5385307/2‑100 (MQ=255)  aGGGCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAacat                                                                                          >  1:4595288/1‑100 (MQ=255)  aGGGCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAacat                                                                                          >  1:581244/1‑100 (MQ=255)   gggCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTGATTAATTTTAATAT‑AAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatga                                                                                        >  1:3438336/1‑100 (MQ=255)   gggCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatg                                                                                         >  1:5083277/1‑100 (MQ=255)   gggCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatg                                                                                         >  1:3080701/1‑100 (MQ=255)   gggCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatg                                                                                         >  1:77044/1‑100 (MQ=255)   cggCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATAT‑AAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatga                                                                                        >  1:3282774/2‑100 (MQ=255)    ggCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatga                                                                                        >  1:1595369/1‑100 (MQ=255)    ggCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatga                                                                                        >  1:4382318/1‑100 (MQ=255)    ggCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATAT‑AAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatgat                                                                                       >  1:1011448/1‑100 (MQ=255)    ctcAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatga                                                                                        >  1:1509492/3‑100 (MQ=255)     gCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatgat                                                                                       >  1:2308641/1‑100 (MQ=255)     gCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatgat                                                                                       >  1:5385240/1‑100 (MQ=255)     gCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatgat                                                                                       >  1:5177561/1‑100 (MQ=255)     acAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatgat                                                                                       >  1:3896962/2‑100 (MQ=255)      cAAGTTGACGCTTCAGATATACATGATCACAGAATAGACTTAATTAATTTTAATATAAAAAAATTACGATCAACAAATATTAAAGCTTTTCAACatgatg                                                                                      >  1:852413/1‑100 (MQ=255)      cAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACatgatg                                                                                      >  1:1266505/1‑100 (MQ=255)          ctGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGAc                                                                                  >  1:3981709/2‑100 (MQ=255)          ctGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGAc                                                                                  >  1:888574/2‑100 (MQ=255)           tGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACa                                                                                 >  1:3450280/1‑100 (MQ=255)            gACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACaa                                                                                >  1:5503188/1‑100 (MQ=255)             gcGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACaaa                                                                               >  1:1907037/2‑100 (MQ=255)                  tCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACagatcgga                                                                          >  1:4379579/1‑93 (MQ=255)                  tCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCtt                                                                          >  1:262052/1‑100 (MQ=255)                   cAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACagatcagaa                                                                         >  1:155509/1‑92 (MQ=255)                   cAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTa                                                                         >  1:337196/1‑100 (MQ=255)                   cAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTa                                                                         >  1:1405685/1‑100 (MQ=255)                        ataCATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATa                                                                    >  1:3670074/1‑100 (MQ=255)                        ataCATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATa                                                                    >  1:5012060/1‑100 (MQ=255)                          aCATGATCACAAAATATACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATaaa                                                                  >  1:1032363/1‑100 (MQ=255)                          aCATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGGCGCGACAAAACCTTATGATaaa                                                                  >  1:4680165/1‑100 (MQ=255)                          aCATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATaaa                                                                  >  1:5088887/1‑100 (MQ=255)                          aCATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATaaa                                                                  >  1:4261008/1‑100 (MQ=255)                             tGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACa                                                               >  1:3997222/1‑100 (MQ=255)                              gATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACAt                                                              >  1:866926/1‑100 (MQ=255)                              gATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACAt                                                              >  1:518497/1‑100 (MQ=255)                              gATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACAt                                                              >  1:3758586/1‑100 (MQ=255)                              gATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACAt                                                              >  1:3863488/1‑100 (MQ=255)                              gATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACAt                                                              >  1:2651310/1‑100 (MQ=255)                                 cacaAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACg                                                           >  1:2461106/1‑100 (MQ=255)                                   caAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGAt                                                         >  1:534651/1‑100 (MQ=255)                                   caAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGAt                                                         >  1:1685603/1‑100 (MQ=255)                                    aaaaTAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATAAGATa                                                        <  1:1475771/100‑1 (MQ=255)                                      aaTAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAg                                                      >  1:2183112/1‑100 (MQ=255)                                         aGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATa                                                   >  1:644790/1‑100 (MQ=255)                                                                     gtACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAg                       >  1:2039385/2‑100 (MQ=255)                                                                      tACGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAgg                      >  1:5253038/1‑100 (MQ=255)                                                                       aCGATTAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGgt                     >  1:2074736/1‑100 (MQ=255)                                                                            tAACAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAgatcggaa                >  1:2402278/1‑93 (MQ=255)                                                                             aaCAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGGTGTAATg               >  1:1025907/1‑100 (MQ=255)                                                                             aaCAAATATTAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGGTGTAATg               <  1:2309609/100‑1 (MQ=255)                                                                                      tAAAGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGGTGTAATGAGACATAAg      >  1:5198189/1‑100 (MQ=255)                                                                                        aaGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGGTGTAATGAGACATAGGcc    >  1:2973522/1‑100 (MQ=255)                                                                                        aaGCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGGTGTAATGAGACATACGcc    >  1:154681/1‑100 (MQ=255)                                                                                          gCTTTTCAACATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGGTGTAATGAGACATAAGccag  <  1:570422/100‑3 (MQ=255)                                                                                                   |                                                                                           AAGGGCAAGTTGACGCTTCAGATATACATGATCACAAAATAGACTTAATTAATTTTAATATAAAAAAATTACGATTAACAAATATTAAAGCTTTTCAATATGATGCGACAAAACCTTATGATAAAACATACGATAAGATACTTGTTGATGCACCATGTAGCGGATTAGGTGTAATGAGACATAAGCCGG  >  NC_007793/1215815‑1216003

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