Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	919,052	A→G	intergenic (‑41/‑60)	SAUSA300_RS04535 ← / → SAUSA300_RS04540	hypothetical protein/disulfide oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	919,052	0	A	G	100.0%	108.5 / NA	40	intergenic (‑41/‑60)	SAUSA300_RS04535/SAUSA300_RS04540	hypothetical protein/disulfide oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/23);  total (17/23)

AAGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAG  >  NC_007793/918955‑919145                                                                                                  |                                                                                              aaGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTttgtt                                                                                             <  1:4132340/100‑1 (MQ=255) aaGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTttgtt                                                                                             <  1:1371152/100‑1 (MQ=255) aaGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTttgtt                                                                                             <  1:4987609/100‑1 (MQ=255)  aGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTttgttt                                                                                            <  1:721835/100‑2 (MQ=255)    aCGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGt                                                                                          <  1:3252685/100‑1 (MQ=255)      gTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGtat                                                                                        <  1:3825345/100‑1 (MQ=255)      gTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGtat                                                                                        <  1:2703352/100‑1 (MQ=255)              aCTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTataat                                                                                <  1:454861/100‑1 (MQ=255)              aCTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTataat                                                                                <  1:5412422/100‑1 (MQ=255)                ttCTGCTACTTGATCAAACATCGTTGTAACTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAatag                                                                              <  1:5182290/100‑1 (MQ=255)                                         gTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAAc                                                     >  1:3125040/1‑100 (MQ=255)                                              ctCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATg                                                >  1:4674064/2‑100 (MQ=255)                                               tCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATgg                                               >  1:3346655/1‑100 (MQ=255)                                               tCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATgg                                               >  1:4811305/1‑100 (MQ=255)                                               tCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATgg                                               >  1:4187175/1‑100 (MQ=255)                                                cAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATggg                                              >  1:1217296/1‑100 (MQ=255)                                                cAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATggg                                              <  1:4084533/100‑1 (MQ=255)                                                cAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATggg                                              <  1:2371199/100‑1 (MQ=255)                                                cAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAACggg                                              >  1:5132303/1‑100 (MQ=255)                                                  gTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATggggg                                            <  1:2681878/100‑1 (MQ=255)                                                  gTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATggggg                                            <  1:1891590/100‑1 (MQ=255)                                                  gTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATggggg                                            <  1:136532/100‑1 (MQ=255)                                                  gTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATggggc                                            <  1:5154628/100‑2 (MQ=255)                                                   tAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGa                                           <  1:5438718/100‑1 (MQ=255)                                                   tAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGa                                           >  1:2640139/1‑100 (MQ=255)                                                   tAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAACGGGGGa                                           >  1:4405315/1‑100 (MQ=255)                                                    aGGCATCAGGCAAGTCTCCTTTCGAGAGGTGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATggggaag                                          >  1:5603882/1‑97 (MQ=255)                                                    aGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCaaggagagat                                          >  1:2787892/1‑92 (MQ=255)                                                    aGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGAt                                          <  1:1274090/100‑1 (MQ=255)                                                    aGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGAt                                          >  1:4189222/1‑100 (MQ=255)                                                          cAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAatat                                    >  1:3675587/1‑100 (MQ=255)                                                          cAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAatat                                    >  1:1483228/1‑100 (MQ=255)                                                          cAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAatat                                    <  1:5396939/100‑1 (MQ=255)                                                           aGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAatata                                   >  1:2119927/1‑100 (MQ=255)                                                           aGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAatata                                   >  1:1284165/1‑100 (MQ=255)                                                              cAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATgg                                <  1:3468150/100‑1 (MQ=255)                                                                         tCGAGAGATGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAgt                     >  1:5443582/1‑100 (MQ=255)                                                                              agaTGAATTTGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCCATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGt                <  1:623094/100‑1 (MQ=255)                                                                                      ttGATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATg        <  1:4036580/100‑1 (MQ=255)                                                                                        gATAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATggg      <  1:3073231/100‑1 (MQ=255)                                                                                          tAATCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGc    <  1:24282/100‑1 (MQ=255)                                                                                          cgaTCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGc    <  1:1555602/98‑1 (MQ=255)                                                                                            aTCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCaa  <  1:974948/100‑2 (MQ=255)                                                                                            aTCTTTGTTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAg  <  1:5253231/100‑1 (MQ=255)                                                                                                  |                                                                                              AAGACGTTCAATAACTTCTGCTACTTGATCAAACATCGTTGTATCTTCAGTAGGCATCAGGCAAGTCTCCTTTCGAGAGATGAATTTGATAATCTTTATTAGTATATTATAATAGATAGTGTAATAAAAATCTATCAAACCAATGGGGGATAAATATATGGAGCACGTGAGTGTGGTAGTTTATGGGGCAG  >  NC_007793/918955‑919145

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