Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	1,632,353	T→C	E25E (GAA→GAG)	SAUSA300_RS08070 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	1,632,353	0	T	C	95.3%	108.3 / ‑2.6	43	E25E (GAA→GAG)	SAUSA300_RS08070	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (5/36);  total (5/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCTTTTTACG  >  NC_007793/1632280‑1632449                                                                          |                                                                                                 cAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGc                                                                        <  1:4721570/100‑1 (MQ=255)                                 tcactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATa                                        >  1:2946612/7‑100 (MQ=255)                                 tcactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATa                                        >  1:2390683/7‑100 (MQ=255)                                  cactcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATAc                                       >  1:4568314/6‑100 (MQ=255)                                   actcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACa                                      <  1:2683281/96‑1 (MQ=255)                                    ctcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACaa                                     <  1:5736548/97‑1 (MQ=255)                                    ctcaTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACaa                                     >  1:2131754/4‑100 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                                  <  1:2720558/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                                  <  1:2887236/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                                  <  1:2204886/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                                  <  1:484466/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                                  <  1:521497/100‑1 (MQ=255)                                       aTACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAAcc                                  <  1:1147150/100‑1 (MQ=255)                                        tACTTTTTCTTTCAGAATTTAAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCa                                 <  1:4389652/100‑1 (MQ=255)                                                 tttCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCtagt                        <  1:2599159/100‑4 (MQ=255)                                                    cagaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGTCCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGAt                     <  1:5182523/90‑1 (MQ=255)                                                    cagaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACGTGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGAt                     <  1:3133114/90‑1 (MQ=255)                                                    cagaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGAt                     <  1:4723097/90‑1 (MQ=255)                                                     agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa                    <  1:748042/91‑1 (MQ=255)                                                     agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa                    <  1:499003/91‑1 (MQ=255)                                                     agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa                    <  1:5320058/91‑1 (MQ=255)                                                     agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa                    <  1:4556545/91‑1 (MQ=255)                                                     agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa                    <  1:3958322/91‑1 (MQ=255)                                                     agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa                    <  1:5750230/91‑1 (MQ=255)                                                     agaatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATa                    <  1:3093557/91‑1 (MQ=255)                                                      gaatttaaaaCTGCATCAATCTCTTTAGCATATTTCTGACCACATGCAGATAATAATATAACTAAGCCTAAAATTAATACAAGCCATCGCCTCATGATaa                   <  1:5477808/92‑1 (MQ=255)                                                      gaatttGAAACTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaa                   <  1:3280461/97‑1 (MQ=255)                                                       aatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaag                  <  1:3838150/93‑2 (MQ=255)                                                       aatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaag                  <  1:1675201/93‑2 (MQ=255)                                                       aatttaaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATaaa                  <  1:1627176/93‑1 (MQ=255)                                                           taaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCa              <  1:2054829/97‑1 (MQ=255)                                                            aaaaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCAc             <  1:4724733/98‑1 (MQ=255)                                                            aaaaCTGCATCAATCTCTTTACCATATTTCTGCCCACATGCATATAATAATATAACTAATCCTAAAACTAATACAATCCATCGCCTCATGATAAATTCAc             <  1:1246990/98‑1 (MQ=255)                                                              aaCTGCATCAATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACtg           <  1:4634092/100‑2 (MQ=255)                                                                 tgcatcaATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCtt        <  1:984235/98‑1 (MQ=255)                                                                 tgcatcaATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCtg        <  1:5047446/98‑2 (MQ=255)                                                                  gcatcaATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCttt       <  1:474679/99‑1 (MQ=255)                                                                  gcatcaATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCttt       >  1:1610698/2‑100 (MQ=255)                                                                   catcaATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCtttt      <  1:3985602/100‑1 (MQ=255)                                                                   catcaATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCtttt      <  1:3445915/100‑1 (MQ=255)                                                                   catcaATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCttgt      <  1:4625396/100‑3 (MQ=255)                                                                     tcaATCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCtttttt    <  1:3891646/100‑2 (MQ=255)                                                                       aaTCTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCtttttttt  <  1:904919/100‑4 (MQ=255)                                                                          |                                                                                                 CAACATTTTTATACTTCTTTGTATTTTTCACATCTTGAATACTTTCTTTTTCTAACTTTGTAACTTCATCAATTTCTTTATCATATTTCTGACCACATGCAGATAATAATATAACTAATCCTAAAATTAATACAAACCATCGCCTCATGATAAATTCACTCCTTTTTACG  >  NC_007793/1632280‑1632449

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