Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,633,220	T→C	K54K (AAA→AAG)	unknown ←	Mobile element protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,633,220	0	T	C	80.7%	47.5 / 5.7	31	K54K (AAA→AAG)	unknown	Mobile element protein
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (11/14);  total (14/17)
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.79e-01

CGGTTCGCTAAATCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCCTCGTTTTCCATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGA  >  1/1633143‑1633469                                                                              |                                                                                                                                                                                                                                                          cGGTTCGCTAAATCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTGGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGttt                                                                              <  1:442241/251‑1 (MQ=9)   gTTCGCTAAATCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTg                                                                             <  1:234831/250‑1 (MQ=14)            aTCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGAtt                                                                   >  1:10336/1‑251 (MQ=255)             tCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGAttt                                                                  >  1:692983/1‑251 (MQ=255)                    aCGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGc                                                           <  1:863768/251‑1 (MQ=255)                    aCGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGa                                                                     <  1:285616/241‑1 (MQ=25)                     cGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGct                                                          <  1:553444/251‑1 (MQ=255)                            aaaTAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGtt                                                                               >  1:30635/1‑223 (MQ=14)                            aaaTAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCtttt                                                   >  1:899815/1‑251 (MQ=255)                                  agCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATc                                              <  1:382515/250‑1 (MQ=255)                                      tttCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTgg                                         >  1:463951/1‑251 (MQ=255)                                       ttcttcCCCAACTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATc                                              >  1:464094/1‑245 (MQ=255)                                                 aCTCGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTAttt                              >  1:219501/1‑251 (MQ=255)                                                    cGAAACTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAg                           <  1:359153/251‑1 (MQ=255)                                                      aaaCTCAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAggg                         <  1:62761/251‑1 (MQ=255)                                                           cAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGa                                                                     <  1:573677/202‑1 (MQ=255)                                                           cAGTGACATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTAtttt                             >  1:607384/1‑242 (MQ=255)                                                             gTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAg                           >  1:554589/1‑242 (MQ=255)                                                                aCATCGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAggg                <  1:481328/250‑1 (MQ=255)                                                                    cGGTTACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtat           >  1:663511/1‑251 (MQ=255)                                                                       ttACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat         <  1:851427/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat         <  1:771553/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat         <  1:901454/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat         <  1:366529/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat         <  1:365770/250‑1 (MQ=255)                                                                       ttACCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTtatat         <  1:222304/250‑1 (MQ=255)                                                                         aCCCACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAg      >  1:669858/1‑251 (MQ=255)                                                                          cccATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAgg                          <  1:37554/230‑1 (MQ=255)                                                                          cccACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCTCATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAgg     >  1:218177/1‑251 (MQ=255)                                                                          cccACTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAgg     >  1:510499/1‑251 (MQ=255)                                                                             aTTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATCCGGGGCAATATGCGATGTTTTGCTTCGTTTTCTATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGa  >  1:558898/1‑251 (MQ=255)                                                                              |                                                                                                                                                                                                                                                          CGGTTCGCTAAATCCATCAACGGTGAAAAATAGAGCTTTTCTTCCCCAACTCGAAACTCAGTGACATCGGTTACCCATTTTTGATTGAGCGCCGTCGCTGTAAAATCACGATTTAGCACATTCGGGGCAATATGCGATGTTTTGCCTCGTTTTCCATGTCTTTTCTTGCGTAAAATAGAATGAATACCTAACGCATTCATCAGTTTTAACACCGTTTTATGATTCAAATGGAAACCCATTTGACGTAGTTTGTTCCGGGATTTCAATCGCTTTCTTTTTCATCTTTGGCTTCACTATTTTAGGGCGTTGAGGGTTTATATAAGGGGA  >  1/1633143‑1633469

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