Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,623,555	G→T	intergenic (‑61/+1361)	SAUSA300_RS13500 ← / ← SAUSA300_RS13505	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,623,555	0	G	T	100.0%	147.5 / NA	52	intergenic (‑61/+1361)	SAUSA300_RS13500/SAUSA300_RS13505	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (4/48);  total (4/48)

GTGTTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTGGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTAATTCATTAACTGGATCTTTTGTGATTTCTTCT  >  NC_007793/2623459‑2623651                                                                                                  |                                                                                                 gtgtTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGt                                                                                               <  1:1523565/100‑1 (MQ=255) gtgtTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGt                                                                                               <  1:1842435/100‑1 (MQ=255) gtgtTTAATCACATCATCCACTGGCTCTT‑CGATCACTTCTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGt                                                                                               <  1:3909285/100‑1 (MQ=255)  tgtTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTa                                                                                              <  1:1689161/100‑1 (MQ=255)  tgtTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTa                                                                                              <  1:4188403/100‑1 (MQ=255)  tgtTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTa                                                                                              <  1:3428537/100‑1 (MQ=255)  tgtTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTa                                                                                              <  1:2327602/100‑1 (MQ=255)   gtTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTaa                                                                                             <  1:3718442/100‑1 (MQ=255)         tCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTgg                                                                                       <  1:4670115/100‑1 (MQ=255)          cacaTCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTAGa                                                                                      >  1:2748253/1‑100 (MQ=255)           acaTCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTggcg                                                                                     <  1:5647119/100‑3 (MQ=255)              tttgtggtgtgCTCTTCCGATCTCTTTACCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCa                                                                                   >  1:3425352/11‑100 (MQ=25)               gatcCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGATCGGaagagc                                                                                 >  1:5538255/2‑95 (MQ=38)                            tcttccgatcTTTTCCGGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATg                                                                     <  1:2235071/91‑1 (MQ=38)                             t‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGAc                                                                   >  1:1337111/1‑100 (MQ=255)                                  atctttttTCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCtta                                                               <  1:3279384/96‑2 (MQ=38)                                   cACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGc                                                              <  1:407118/100‑1 (MQ=255)                                    cttgttCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGca                                                             <  1:4230036/96‑2 (MQ=255)                                     cTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCgg                                                            <  1:2419975/100‑1 (MQ=255)                                     cTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCgg                                                            <  1:393086/100‑1 (MQ=255)                                      ttttCGTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGt                                                           <  1:301652/100‑1 (MQ=255)                                      ttttCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGt                                                           <  1:3622690/100‑1 (MQ=255)                                      ttttCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGt                                                           <  1:2209935/100‑1 (MQ=255)                                       tttCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTa                                                          <  1:2228889/100‑1 (MQ=255)                                       actCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCggag                                                          <  1:2369352/98‑3 (MQ=255)                                        tttcTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTAt                                                         <  1:4362739/97‑1 (MQ=255)                                        atCCTGTGTCTGGATTCTTGATTACTGGTTTACCTGGTACTTTATCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGtgg                                                         <  1:5295500/99‑3 (MQ=255)                                         tCCTGTGTCTGGATTCTTGATTCCTTGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTAtt                                                        <  1:3507605/100‑1 (MQ=255)                                                 cTGGATTCTTGATTCCTGGTTTGCCTGGTACTTTTTCCGTATGATCTGTTGGTAAGTTTGGAACAAAGATATCTTTATGACCTTGCGATATTTTCTCGcc                                                <  1:3704824/100‑1 (MQ=255)                                                  tGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCa                                               <  1:2975351/100‑1 (MQ=255)                                                    gATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCAct                                             <  1:4326288/100‑2 (MQ=255)                                                    gATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCAca                                             <  1:5289187/100‑2 (MQ=255)                                                     aTTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCActg                                            <  1:4115965/100‑3 (MQ=255)                                                     aTTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCAcac                                            <  1:4935846/100‑3 (MQ=255)                                                      ttCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCAcaaa                                           <  1:1657428/100‑4 (MQ=255)                                                            aTTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTg                                     <  1:3198087/100‑1 (MQ=255)                                                             ttCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGt                                    <  1:811392/100‑1 (MQ=255)                                                             ttCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGt                                    <  1:3097995/100‑1 (MQ=255)                                                              tCCTGGTTTGCCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATGCTGtt                                   <  1:244814/100‑1 (MQ=255)                                                              tCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGtt                                   <  1:1904557/100‑1 (MQ=255)                                                              tCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGtt                                   <  1:2396381/100‑1 (MQ=255)                                                              tCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGtt                                   <  1:58700/100‑1 (MQ=255)                                                               ccTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGGTATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTa                                  <  1:2898312/100‑1 (MQ=255)                                                               ccTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTa                                  <  1:2098321/100‑1 (MQ=255)                                                               ccTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTa                                  <  1:3542448/100‑1 (MQ=255)                                                               ccTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTa                                  <  1:4727811/100‑1 (MQ=255)                                                               ccTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTa                                  <  1:1247406/100‑1 (MQ=255)                                                                cTGGTTTACCTGGTACTTTTTCCGTTTGACCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTaa                                 <  1:4109999/100‑1 (MQ=255)                                                                            gTACTTTTTCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTAATTCATTAACTgg                     <  1:999378/100‑1 (MQ=255)                                                                                    tCCGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTAATTCATTAACTGGATCTTTtg             <  1:2929013/100‑1 (MQ=255)                                                                                      cGTTTGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTAATTCATTAACTGGATCTTTtgtg           <  1:5659763/100‑1 (MQ=255)                                                                                        tttGATCTGTTGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCGCCGAATTCTGTTAATTCATTAACTGGATCTTTTGTGAt         <  1:2639434/100‑1 (MQ=255)                                                                                               tgttGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTAATTCATTAACTGGATCTTTTGTGATTtcttct  <  1:1522792/100‑1 (MQ=255)                                                                                                  |                                                                                                 GTGTTTAATCACATCATCCACTGGCTCTT‑CGATCACTTTTCCTGTGTCTGGATTCTTGATTCCTGGTTTACCTGGTACTTTTTCCGTTTGATCTGTGGGTAAGTTTGGATCAAAGATATCTTTATGACCTTGCGGTATTTTCTCGCCACCGAATTCTGTTAATTCATTAACTGGATCTTTTGTGATTTCTTCT  >  NC_007793/2623459‑2623651

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: ‑