Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,686,226	C→A	D70Y (GAT→TAT)	SAUSA300_RS08375 ←	tRNA (N(6)‑L‑threonylcarbamoyladenosine(37)‑C(2))‑ methylthiotransferase MtaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,686,226	0	C	A	100.0%	155.2 / NA	51	D70Y (GAT→TAT)	SAUSA300_RS08375	tRNA (N(6)‑L‑threonylcarbamoyladenosine(37)‑C(2))‑ methylthiotransferase MtaB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (6/45);  total (6/45)

TATCTTGTGTACCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATCAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATCAGCATT  >  NC_007793/1686130‑1686322                                                                                                 |                                                                                                 tATCTTGTGTACCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAgg                                                                                               <  1:721977/100‑1 (MQ=255) tATCTTGTGTACCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAgg                                                                                               <  1:5130916/100‑1 (MQ=255)     ttGTGTACCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGAtta                                                                                           <  1:3669915/100‑2 (MQ=255)      tgtgTACCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGAtttt                                                                                          >  1:5284653/1‑100 (MQ=255)       gtgtACCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTg                                                                                         <  1:613050/100‑1 (MQ=255)       gtgtACCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTg                                                                                         <  1:5974570/100‑1 (MQ=255)           aCCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCtt                                                                                     <  1:3292922/100‑1 (MQ=255)                  actacATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCAc                                                                              >  1:675587/1‑100 (MQ=255)                  actacATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCAc                                                                              >  1:4978274/1‑100 (MQ=255)                  actacATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTACTGcca                                                                              <  1:314247/100‑3 (MQ=255)                      cATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCt                                                                          <  1:4340263/100‑1 (MQ=255)                       aTCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTa                                                                         <  1:4203773/100‑1 (MQ=255)                       aTCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTa                                                                         <  1:804157/100‑1 (MQ=255)                            cacCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTaggact                                                                    <  1:1274763/100‑6 (MQ=255)                             tccAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAattatt                                                                   >  1:1473878/2‑100 (MQ=255)                             acCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAattatt                                                                   >  1:1776375/1‑100 (MQ=255)                               cAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTg                                                                 <  1:1083312/100‑1 (MQ=255)                               cAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTg                                                                 <  1:1041204/100‑1 (MQ=255)                                              tttCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATc                                                  <  1:1421776/100‑1 (MQ=255)                                                  aGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAGGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCAccc                                              <  1:5823612/100‑1 (MQ=255)                                                  aGCTGATGAAGTTTGCGCATAACAACCTATTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCAcca                                              <  1:1568470/100‑2 (MQ=255)                                                   gCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCg                                             <  1:5048560/100‑1 (MQ=255)                                                     tgatgaAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGtg                                           <  1:752512/100‑2 (MQ=255)                                                     tgatgaAGTTTGCGCATAACAACCTCTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTa                                           <  1:4629739/100‑1 (MQ=255)                                                       atgaAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTAtc                                         <  1:3324940/100‑2 (MQ=255)                                                         gaAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCgccagca                                       <  1:868551/100‑7 (MQ=255)                                                          aaGTTTGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTAtcctt                                      <  1:3112129/100‑5 (MQ=255)                                                             tttGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTAc                                   <  1:432224/100‑1 (MQ=255)                                                              ttGCGCATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACt                                  <  1:1215109/100‑1 (MQ=255)                                                                 cgcATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGtc                               <  1:4299949/100‑2 (MQ=255)                                                                  gcATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTAc                              <  1:2461391/100‑1 (MQ=255)                                                                  gcATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTAc                              <  1:209949/100‑1 (MQ=255)                                                                  gcATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTAc                              <  1:1577711/100‑1 (MQ=255)                                                                  gcATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTAc                              <  1:5426930/100‑1 (MQ=255)                                                                   cATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACa                             <  1:256250/100‑1 (MQ=255)                                                                   cATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACa                             <  1:3544874/100‑1 (MQ=255)                                                                   cATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACa                             <  1:1583526/100‑1 (MQ=255)                                                                   cATAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACa                             <  1:5434034/100‑1 (MQ=255)                                                                    cttccgatcTGTTGCGCAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTActg                            <  1:261786/92‑3 (MQ=255)                                                                    aTAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTgccaag                            <  1:1816781/100‑6 (MQ=255)                                                                     tAACAACCTGTTACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAg                           <  1:3488716/100‑1 (MQ=255)                                                                                tACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACa                >  1:5163032/1‑100 (MQ=255)                                                                                tACACAGATTACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACa                <  1:5950878/100‑1 (MQ=255)                                                                                        ttACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATc        <  1:2956923/100‑1 (MQ=255)                                                                                        ttACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATc        <  1:3872940/100‑1 (MQ=255)                                                                                        ttACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATc        <  1:3944164/100‑1 (MQ=255)                                                                                         tACAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATCa       <  1:2620675/100‑1 (MQ=255)                                                                                           cAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATCAga     <  1:3293620/100‑2 (MQ=255)                                                                                           cAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATCAGc     <  1:779060/100‑1 (MQ=255)                                                                                           cAGCATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATCAGc     <  1:5620130/100‑1 (MQ=255)                                                                                              cATAAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGGCTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATCAGCagg  <  1:277368/100‑3 (MQ=255)                                                                                                 |                                                                                                 TATCTTGTGTACCAACTACTACATCGACACCAGGAATTTCCATAATTTCAGCTGATGAAGTTTGCGCATAACAACCTGTTACACAGATTACAGCATCAGGATTTTGTCTTATTGCACGTCTAATTATTTGACGACTTTTTTTATCACCCGTATTCGTTACTGTACAAGTATTAATAACAAATACATCAGCATT  >  NC_007793/1686130‑1686322

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