Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,693,048	(T)7→8	coding (1124/1125 nt)	SAUSA300_RS08410 ←	coproporphyrinogen III oxidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,693,041	1	.	T	100.0%	115.0 / NA	39	intergenic (‑94/+6)	SAUSA300_RS08405/SAUSA300_RS08410	HrcA family transcriptional regulator/coproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (3/36);  total (3/36)

TTAGCACTCACTTATCTCAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAA‑‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCGTTCTTTTCTACAA  >  NC_007793/1692960‑1693135                                                                                   |                                                                                                ataccactcactTATCTCAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATttat                                                                               >  1:465744/5‑100 (MQ=255)       ctcactTATCTCAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGaa                                                                         >  1:5192903/1‑100 (MQ=255)             tATCTCAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTc                                                                   >  1:700100/1‑100 (MQ=255)                             tataATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCt                                                   <  1:2882324/100‑1 (MQ=255)                                aaTTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTagg                                                <  1:3413340/100‑3 (MQ=255)                                    tATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACtt                                            <  1:4321386/100‑1 (MQ=255)                                    tATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACtt                                            <  1:2915546/100‑1 (MQ=255)                                      tCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACtttg                                          <  1:4611436/100‑2 (MQ=255)                                         aaTTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCct                                       <  1:1256918/100‑1 (MQ=255)                                                 aaaGTAAGTCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGt                               <  1:1186307/100‑1 (MQ=255)                                                  aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:5308559/100‑1 (MQ=255)                                                  aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:1617810/100‑1 (MQ=255)                                                  aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:5464210/100‑1 (MQ=255)                                                  aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:927175/100‑1 (MQ=255)                                                  aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:4898958/100‑1 (MQ=255)                                                  aagtaagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:401249/100‑1 (MQ=255)                                                  aaataagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:2356869/97‑1 (MQ=255)                                                   agtaagtCAATGTTAAAGACTCGAAATTTCAATTTTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:2013250/100‑1 (MQ=255)                                                   agtaagtCAATGTTAAAGACTCGAAATTTCAATTTTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:680124/100‑1 (MQ=255)                                                   agtaagtCAATGTTAAAGACTCGAAATTTCAATTTTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTa                              <  1:3811747/100‑1 (MQ=255)                                                     taagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAgg                           <  1:2074154/100‑2 (MQ=255)                                                      aagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTCCCTATGACTTTCCCTCTATTTGTAAGta                          <  1:1052360/100‑2 (MQ=255)                                                      aagtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTg                          <  1:4997891/100‑1 (MQ=255)                                                       agtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                         <  1:820266/100‑1 (MQ=255)                                                       agtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                         <  1:393715/100‑1 (MQ=255)                                                       agtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                         <  1:5063989/100‑1 (MQ=255)                                                       agtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGc                         <  1:1432636/100‑1 (MQ=255)                                                        gtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                        <  1:4729808/100‑1 (MQ=255)                                                        gtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                        <  1:3570751/100‑1 (MQ=255)                                                        gtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                        <  1:2653343/100‑1 (MQ=255)                                                        gtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                        <  1:5375872/100‑1 (MQ=255)                                                        gtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                        <  1:797101/100‑1 (MQ=255)                                                        gtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                        <  1:1292662/100‑1 (MQ=255)                                                        gtCAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCa                        <  1:764983/100‑1 (MQ=255)                                                          cAATGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCatt                      <  1:5130698/100‑3 (MQ=255)                                                           aaTGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTGTGACTTTCCCTCTATTTGTAAGTGCAATc                     <  1:1174481/100‑1 (MQ=255)                                                           aaTGTTAAAGACTCGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATc                     <  1:1376948/100‑1 (MQ=255)                                                                        cGAAATTTCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCGTTCttgc        <  1:4081257/100‑3 (MQ=255)                                                                              ttCAAT‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCGTTCTTTTCcacat  <  1:5516852/100‑2 (MQ=255)                                                                                   |                                                                                                TTAGCACTCACTTATCTCAAGTGCTAATTATAATTTATCAAATTGGTCAAAGTAAGTCAATGTTAAAGACTCGAAATTTCAA‑‑TTTTTTTAATCATTTATTAGGAAAGCTTCAAAAACCTCATTACCTATGACTTTCCCTCTATTTGTAAGTGCAATCACATCGTTCTTTTCTACAA  >  NC_007793/1692960‑1693135

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