Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	739,038	A→G	K54K (AAA→AAG)	unknown →	Mobile element protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	739,038	0	A	G	80.7%	30.9 / 6.3	31	K54K (AAA→AAG)	unknown	Mobile element protein
Reads supporting (aligned to +/- strand):  ref base A (3/3);  new base G (14/11);  total (17/14)
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

TCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATGGAAAACGAGGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGATTTAGCGAACCG  >  1/738789‑739115                                                                                                                                                                                                                                                          |                                                                              tCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGGATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAAt                                                                              <  1:558898/251‑1 (MQ=14)    ccTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTggg                                                                           <  1:510499/251‑1 (MQ=14)    ccTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGAGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTggg                                                                           <  1:218177/251‑1 (MQ=9)     cTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGt                                                                          <  1:669858/251‑1 (MQ=14)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                                        >  1:901454/1‑250 (MQ=14)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                                        >  1:771553/1‑250 (MQ=14)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                                        >  1:366529/1‑250 (MQ=14)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                                        >  1:365770/1‑250 (MQ=14)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                                        >  1:222304/1‑250 (MQ=14)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTaa                                                                        >  1:851427/1‑250 (MQ=14)          ataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCg                                                                     <  1:663511/251‑1 (MQ=14)               cccTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGt                                                                 >  1:481328/1‑250 (MQ=14)                        cccTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGttt                                                       >  1:62761/1‑251 (MQ=14)                         ccTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATggg                                                                           >  1:37554/1‑230 (MQ=14)                          cTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCg                                                     >  1:359153/1‑251 (MQ=14)                          cTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCAc                                                              <  1:554589/242‑1 (MQ=14)                            aaaaTAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTg                                                            <  1:607384/242‑1 (MQ=14)                             aaaTAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGt                                                  <  1:219501/251‑1 (MQ=14)                                        ccAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGaaa                                       <  1:463951/251‑1 (MQ=14)                                             gATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGgaagaa                                        <  1:464094/245‑1 (MQ=14)                                             gATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGct                                   >  1:382515/1‑250 (MQ=14)                                                  aaaaGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTAttt                             <  1:899815/251‑1 (MQ=14)                                                         agCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCg                      >  1:553444/1‑251 (MQ=14)                                                          gCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGt                     >  1:863768/1‑251 (MQ=14)                                                                 aaaTCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGa              <  1:692983/251‑1 (MQ=14)                                                                  aaTCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGAt             <  1:10336/251‑1 (MQ=14)                                                                    tCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTg                                                            >  1:573677/1‑202 (MQ=9)                                                                    tCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGt                     >  1:285616/1‑241 (MQ=14)                                                                            cAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGATTTAGCGAAc    >  1:234831/1‑250 (MQ=14)                                                                             aaaCTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCCAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGATTTAGCGAACCg  >  1:442241/1‑251 (MQ=9)                                                                              aaCTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTAttt                             <  1:30635/223‑1 (MQ=14)                                                                                                                                                                                                                                                          |                                                                              TCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATGGAAAACGAGGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGATTTAGCGAACCG  >  1/738789‑739115

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: ‑