Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,380,664	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	1,380,664	0	A	G	82.1%	52.4 / 7.0	28	K54K (AAA→AAG)	unknown	Mobile element protein
Reads supporting (aligned to +/- strand):  ref base A (3/2);  new base G (12/11);  total (15/13)
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

TCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTTCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGAT  >  1/1380415‑1380730                                                                                                                                                                                                                                                          |                                                                   tCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGGATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAAt                                                                   <  1:558898/251‑1 (MQ=255)    ccTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTggg                                                                <  1:510499/251‑1 (MQ=255)    ccTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGAGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTggg                                                                <  1:218177/251‑1 (MQ=255)     cTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGt                                                               <  1:669858/251‑1 (MQ=255)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                             >  1:901454/1‑250 (MQ=255)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                             >  1:771553/1‑250 (MQ=255)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                             >  1:366529/1‑250 (MQ=255)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                             >  1:222304/1‑250 (MQ=255)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTaa                                                             >  1:365770/1‑250 (MQ=255)        atataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTaa                                                             >  1:851427/1‑250 (MQ=255)          ataAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCg                                                          <  1:663511/251‑1 (MQ=255)               cccTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGt                                                      >  1:481328/1‑250 (MQ=255)                        cccTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGttt                                            >  1:62761/1‑251 (MQ=255)                         ccTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATggg                                                                >  1:37554/1‑230 (MQ=255)                          cTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCg                                          >  1:359153/1‑251 (MQ=255)                          cTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCAc                                                   <  1:554589/242‑1 (MQ=255)                            aaaaTAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTg                                                 <  1:607384/242‑1 (MQ=255)                             aaaTAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGt                                       <  1:219501/251‑1 (MQ=255)                                        ccAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGaaa                            <  1:463951/251‑1 (MQ=255)                                             gATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGgaagaa                             <  1:464094/245‑1 (MQ=255)                                             gATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGct                        >  1:382515/1‑250 (MQ=255)                                                  aaaaGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTAttt                  <  1:899815/251‑1 (MQ=255)                                                         agCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCg           >  1:553444/1‑251 (MQ=255)                                                          gCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGt          >  1:863768/1‑251 (MQ=38)                                                                 aaaTCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGa   <  1:692983/251‑1 (MQ=21)                                                                  aaTCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGAt  <  1:10336/251‑1 (MQ=18)                                                                    tCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTg                                                 >  1:573677/1‑202 (MQ=18)                                                                    tCCCGGAACAAACTACGTCAAATGGGTTTCCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAGTGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGt          >  1:285616/1‑241 (MQ=18)                                                                                                                                                                                                                                                          |                                                                   TCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACTACGTCAAATGGGTTTTCATTTGAATCATAAAACGGTGTTAAAACTGATGAATGCGTTAGGTATTCATTCTATTTTACGCAAGAAAAGACATAGAAAACGAAGCAAAACATCGCATATTGCCCCGAATGTGCTAAATCGTGATTTTACAGCGACGGCGCTCAATCAAAAATGGGTAACCGATGTCACTGAGTTTCGAGTTGGGGAAGAAAAGCTCTATTTTTCACCGTTGATGGAT  >  1/1380415‑1380730

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