Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	76,319	C→T	A17A (GCC→GCT)	SAUSA300_RS00345 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	76,319	0	C	T	95.7%	86.9 / ‑2.9	46	A17A (GCC→GCT)	SAUSA300_RS00345	transposase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (4/40);  total (5/41)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

CAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTA  >  NC_007793/76280‑76411                                        |                                                                                             cAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGaa                                  <  1:5483144/100‑1 (MQ=14)      ataAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtgg                             <  1:4312791/100‑1 (MQ=17)       taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTGCTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                            <  1:2170976/100‑1 (MQ=12)       taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                            <  1:276138/100‑1 (MQ=17)       taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                            <  1:446804/100‑1 (MQ=17)       taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggt                            <  1:4772786/100‑1 (MQ=17)       taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggg                            <  1:1043518/100‑2 (MQ=17)       taAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggc                            <  1:4373098/100‑2 (MQ=17)       taAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggg                            <  1:3982293/100‑2 (MQ=17)        aaaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACgtggtg                           >  1:4161174/1‑100 (MQ=17)         aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGGTATATCTGAAATATTAAGGGAACGTGgtgt                          <  1:336532/100‑1 (MQ=12)         aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                          <  1:1934819/100‑1 (MQ=17)         aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                          <  1:1631172/100‑1 (MQ=17)         aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                          <  1:1017028/100‑1 (MQ=17)         aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                          <  1:5523303/100‑1 (MQ=17)         aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                          <  1:4147623/100‑1 (MQ=17)         aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgt                          <  1:1733533/100‑1 (MQ=17)          aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGgtgtg                         <  1:2046252/100‑2 (MQ=17)          aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                         <  1:2576524/100‑1 (MQ=17)          aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                         <  1:3244414/100‑1 (MQ=17)          aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                         <  1:4205644/100‑1 (MQ=17)          aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                         <  1:4812265/100‑1 (MQ=17)          aCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTa                         <  1:4476784/100‑1 (MQ=17)           cAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTac                        <  1:3301425/100‑2 (MQ=17)           cAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTaa                        <  1:3604105/100‑1 (MQ=17)               gtAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGt                    >  1:3963071/2‑100 (MQ=17)                tAACAAGGATGTTATCACTGTAGCTGTTGGATACTATATAAGATATGAATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGAGGTGTAAACGtt                   >  1:463457/1‑100 (MQ=11)                 aaCAAGAATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTtc                  >  1:3873190/1‑100 (MQ=12)                   cAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTtcat                >  1:1758162/1‑100 (MQ=18)                             aTCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACgg      <  1:2289211/100‑1 (MQ=17)                               cACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTc    <  1:1539199/100‑1 (MQ=18)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACggggc   <  1:2331975/100‑4 (MQ=12)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACgggcc   <  1:3382645/100‑4 (MQ=12)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt   <  1:804930/100‑1 (MQ=17)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt   <  1:5396788/100‑1 (MQ=17)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt   <  1:4865197/100‑1 (MQ=17)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt   <  1:4644517/100‑1 (MQ=17)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt   <  1:4532301/100‑1 (MQ=17)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt   <  1:1789592/100‑1 (MQ=17)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt   <  1:228773/100‑1 (MQ=17)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCt   <  1:3317970/100‑1 (MQ=17)                                aCTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTAATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTAAACGGTCt   <  1:1351343/100‑1 (MQ=255)                                 cTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACggggta  <  1:418502/100‑5 (MQ=12)                                 cTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTa  <  1:4794744/100‑1 (MQ=17)                                 cTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTa  <  1:2462044/100‑1 (MQ=17)                                 cTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTa  <  1:1375069/100‑1 (MQ=17)                                        |                                                                                             CAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTA  >  NC_007793/76280‑76411

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