Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	1,406	C→G	intergenic (+410/‑173)	SAUSA300_RS14705 → / → SAUSA300_RS14710	replication protein/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	1,406	0	C	G	100.0%	81.6 / NA	29	intergenic (+410/‑173)	SAUSA300_RS14705/SAUSA300_RS14710	replication protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (7/22);  total (7/22)

GCTAATTTTT‑AATTTGCATGTAAATGGGCAGT‑GT‑TAAAAAATTAGTCACTGGTTTTGCTGAAATTTTTGTTTGTGAAATTACCATATATTTATTTGGCTCATATTTGCGTTTTAAGAGCTTGTATATATTTTTAAGTTAAAAACCTATATGATAC  >  NC_007790/1325‑1479                                                                                     |                                                                          gCTAATTTTTAAATTTGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGc                                                           <  1:3617246/100‑1 (MQ=255)  cTAATTTTTAAATTTGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCt                                                          <  1:3639112/100‑1 (MQ=255)   tAATTTTTAAATTTGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCtg                                                         <  1:3861449/100‑2 (MQ=255)   tAATTTTTAAATTTGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTc                                                         <  1:4484414/100‑1 (MQ=255)         ttaaaTTTGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATAtta                                                   <  1:5619554/97‑2 (MQ=255)            aaTTTGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCg                                                <  1:4073232/100‑1 (MQ=255)              tttGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCGCTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGtt                                              <  1:1144220/100‑1 (MQ=255)              tttGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGtt                                              <  1:4644747/100‑1 (MQ=255)               cagCATGTAAATGGGCAGT‑GT‑CTAAAAA‑TAGTTATTGTTTTTGCTCAAATTTTTGTTTTTGAAAATAGCATATATTTATTTGGCTCATATTTGCGTTTTa                                           >  1:2404927/3‑100 (MQ=255)                tGCATGTAAATGGGCAGTATT‑TAAAAAATTAGTCACTGATTTTGTTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGtttt                                            <  1:1420309/100‑1 (MQ=255)                tGCATGTAAATGGGCAGT‑GT‑CTAAAAA‑TAGTTATTGTTTTTGCTCAAATTTTTGTTTTTGAAAATAGCATATATTTATTTGGCTCATATTTGCGTTTTaa                                          >  1:2720352/1‑100 (MQ=255)                tGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGACATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGtttt                                            <  1:5005935/100‑1 (MQ=255)                tGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGtttt                                            >  1:2602478/1‑100 (MQ=255)                tGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGtttt                                            <  1:4580945/100‑1 (MQ=255)                tGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGtttt                                            >  1:1253700/1‑100 (MQ=255)                tGCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCCGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGtttt                                            <  1:270719/100‑1 (MQ=255)                 gCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTa                                           <  1:4020880/100‑1 (MQ=255)                 gCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTa                                           <  1:1104449/100‑1 (MQ=255)                 gCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTa                                           <  1:3108343/100‑1 (MQ=255)                 gCATGTAAATGGGCAGT‑ATCTAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTa                                           <  1:2717283/100‑1 (MQ=255)                        caTGGGCAGT‑GT‑CTAAAAA‑TAGTTATTGTTTTTGCTCAAATTTTTGTTTTTGAAAATAGCATATATTTATTTGGCTCATATTTGCGTTTTAAAAGcatat                                  >  1:5651500/2‑96 (MQ=255)                             ggcagtatctAAAAGATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTAAGAGACTGtat                                <  1:3654012/91‑1 (MQ=255)                             ggcagtatctAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTAAGAGACTGtat                                <  1:4492616/91‑1 (MQ=255)                             ggcagtatctAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTAAGAGACTGtat                                <  1:5244403/91‑1 (MQ=255)                                  tatctAAAAAATTAGTCACTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTAAGAGACTGtatagggt                           <  1:4701949/96‑5 (MQ=255)                                    tgtctaaaaATAGTTATTGTTTTTGCTCAAATTTTTGTTTTTGAAAATAGCATATATTTATTTGGCTCATATTTGCGTTTTAAAAGCATATATAAGtttt                         <  1:3492849/95‑1 (MQ=255)                                                    cTGATTTTGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTAAGAGACTGTATAGATTTTTAAGTACAAAACCtag         <  1:2285133/100‑2 (MQ=255)                                                         cttGCTCAAATTTTTGTTTTTGAAAATAGCATATATTTATTTGGCTCATATTTGCGTTTTAAAAGC‑‑ATATATAAGTTTTAGATATAAAACTATATGaatc  >  1:1939808/2‑97 (MQ=255)                                                          ttGCTGAAATTTTTGTTTGTGAAATTAGCATATATTTATTTGGCTCATATTTGCGTTTTAAGAGACTGTATAGATTTTTAAGTACAAAACCTATATGAtt   >  1:1100434/1‑99 (MQ=255)                                                                                     |                                                                          GCTAATTTTT‑AATTTGCATGTAAATGGGCAGT‑GT‑TAAAAAATTAGTCACTGGTTTTGCTGAAATTTTTGTTTGTGAAATTACCATATATTTATTTGGCTCATATTTGCGTTTTAAGAGCTTGTATATATTTTTAAGTTAAAAACCTATATGATAC  >  NC_007790/1325‑1479

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