Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007790	1,871	C→A	intergenic (+104/+222)	SAUSA300_RS14710 → / ← SAUSA300_RS14715	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007790	1,871	0	C	A	100.0%	135.3 / NA	48	intergenic (+104/+222)	SAUSA300_RS14710/SAUSA300_RS14715	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (38/10);  total (38/10)

AGGGGATTTTGAAACGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGCTTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGA  >  NC_007790/1776‑1965                                                                                                |                                                                                                 tgaggATTTTGAAACGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGAtttt                                                                                              >  1:1901976/4‑100 (MQ=255) agaggATTTTGAAACGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGAtttc                                                                                              >  1:5169927/4‑99 (MQ=255) agaggATTTTGAAACGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGAGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGAttct                                                                                              >  1:1107302/4‑98 (MQ=255)  gaggATTTTGAAACGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGAtttgg                                                                                             >  1:1870342/3‑98 (MQ=255)   gggATTTTGAAACGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTAC‑TTTTGTTTTTCTATAAAGATGATGATTTAGAc                                                                                           >  1:4150949/1‑100 (MQ=255)          tGAAACGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGGGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACggagat                                                                                     >  1:3143279/1‑95 (MQ=255)            caaCGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACAAGAtcgg                                                                                   >  1:4833632/2‑97 (MQ=255)            aaaCGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTAGGACGAGAtcgg                                                                                   >  1:83562/1‑97 (MQ=255)                 aGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTGGGACGAGATCGgaagag                                                                              >  1:1622202/1‑99 (MQ=255)                        tcttGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTAAACAAGATCGgaagagcacacgt                                                                       >  1:255243/1‑92 (MQ=255)                        tcttGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTATAACGAGATCGgaagagcacacgt                                                                       >  1:782360/1‑92 (MQ=255)                        tcttGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTATAACGAGATCGgaagagcacacgt                                                                       >  1:5066121/1‑92 (MQ=255)                        tcttGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTAAGACGAGATCGgaagagcacacgt                                                                       >  1:2540359/1‑92 (MQ=255)                        tcttGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACCTTTTGTTTTTCTATAAAGATGATGATTCTGACGAGATCGgaagagcacacgt                                                                       >  1:1233399/1‑92 (MQ=255)                        tcttGTCTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGATCGgaagagcacacgt                                                                       >  1:2170871/1‑92 (MQ=255)                         cttGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTGGGACGAGATCGgaagagcacacgtc                                                                      >  1:4078970/1‑91 (MQ=255)                           tGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTa                                                                   >  1:594830/1‑100 (MQ=255)                                 tacacTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGa                                                             >  1:3314214/2‑100 (MQ=255)                                 cacaCTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGa                                                             >  1:5024863/1‑100 (MQ=255)                                  acacTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGaa                                                            >  1:3643278/1‑100 (MQ=255)                                                    ctctGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAgacg                                          >  1:984520/4‑100 (MQ=255)                                                     cgtGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAgacga                                         >  1:838729/3‑100 (MQ=255)                                                       tGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGacgacg                                       >  1:3275042/1‑100 (MQ=255)                                                            ttgttgCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCAttt                                  <  1:3272494/100‑1 (MQ=255)                                                             tgttgCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTa                                 >  1:5031695/1‑100 (MQ=255)                                                                tgCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGc                              >  1:1832495/1‑100 (MQ=255)                                                                tgCCTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGc                              >  1:265226/1‑100 (MQ=255)                                                                 gCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGcc                             >  1:1420973/1‑100 (MQ=255)                                                                 gCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGcc                             >  1:361555/1‑100 (MQ=255)                                                                 gCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGcc                             >  1:3914242/1‑100 (MQ=255)                                                                 gCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGcc                             >  1:5178879/1‑100 (MQ=255)                                                                 gCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGcc                             >  1:3126470/1‑100 (MQ=255)                                                                 gCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGcc                             >  1:5028522/1‑100 (MQ=255)                                                                 gCTTACTTTTTGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTAGGcc                             >  1:5230204/1‑100 (MQ=255)                                                                           tGTTTTTCTATAAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAAttt                   >  1:4454519/1‑100 (MQ=255)                                                                                   tataAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTAtggcgcg           <  1:4768965/100‑7 (MQ=255)                                                                                    ataAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTTTGCCGAGAAAATTTATTgtgcgt          >  1:4380545/1‑95 (MQ=255)                                                                                    ataAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgt          >  1:4013234/1‑95 (MQ=255)                                                                                    ataAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgt          >  1:2454038/1‑95 (MQ=255)                                                                                    ataAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAAT‑AAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgta         >  1:1348929/1‑94 (MQ=255)                                                                                     taAAGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgtt         <  1:5058021/100‑7 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        <  1:4311978/100‑8 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        <  1:4269823/100‑8 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        <  1:2124015/100‑8 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        <  1:1591496/100‑8 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        <  1:2939228/100‑8 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑agcgttg       <  1:2840281/100‑7 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GGAAAAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        <  1:2676291/100‑8 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GGAAAAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        <  1:3733048/100‑8 (MQ=255)                                                                                      aaaGATGATGATTTTGACGAGAT‑GAAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTgtgcgttg        >  1:2238225/1‑93 (MQ=255)                                                                                           tgatgaTTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTACGTTgagatc   >  1:4470878/1‑98 (MQ=255)                                                                                            gatgatTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGagaaga  >  1:3497991/1‑100 (MQ=255)                                                                                            aatgATTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTGTGCGTTGagaaga  >  1:3410283/2‑100 (MQ=255)                                                                                                |                                                                                                 AGGGGATTTTGAAACGAGTTTCTTCTTGTTTTCACACTGTTTTTTTATTCCTATTGGTGTTGTTGCTTACTTTTTGTTTTTCTATAAAGATGATGCTTTTGACGAGAT‑GGAAGAAAAATATGATTATCATGAAGATAATAAAAAATAGACGACGCATTTATGCCGAGAAAATTTATTG‑ATGTTGAGAAGA  >  NC_007790/1776‑1965

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