Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,555,895	C→A	P277Q (CCA→CAA)	yhhI →	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,555,895	0	C	A	96.8%	83.9 / ‑5.2	31	P277Q (CCA→CAA)	yhhI	predicted transposase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (18/12);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGAT  >  REL606/3555822‑3555961                                                                          |                                                                   gAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAg                                                                  >  2:956628/1‑76 (MQ=21) gAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAg                                                                  <  2:272191/76‑1 (MQ=21)  aaTGGAAAGGGCTGAAGAATTTATGCGTGGCAGTCTCCTTTCGGTCCAT‑ATAGCAGAACAAAAGAAAGAGCAAGaa                                                                >  2:127674/1‑76 (MQ=11)  aaTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGa                                                                 >  2:1500282/1‑76 (MQ=21)   aTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGaa                                                                >  1:1318224/1‑76 (MQ=21)   aTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGaa                                                                <  1:980107/76‑1 (MQ=21)   aTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGaa                                                                >  1:422549/1‑76 (MQ=21)    tGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGaaa                                                               >  1:911350/1‑76 (MQ=21)     ggAAAGGGCTGCAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAAt                                                              >  1:395121/1‑76 (MQ=21)      gAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATg                                                             >  1:1496160/1‑76 (MQ=21)      gAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATg                                                             <  1:1066418/76‑1 (MQ=21)         aGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACg                                                          <  2:1519963/76‑1 (MQ=21)            gCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTc                                                       <  2:1318906/76‑1 (MQ=21)             cTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCa                                                      >  2:183397/1‑76 (MQ=21)              tGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAg                                                     <  2:1201997/76‑1 (MQ=17)              tGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAg                                                     <  2:848264/76‑1 (MQ=17)               gaagaaATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGa                                                    <  1:305794/76‑1 (MQ=14)                                  gTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTctgctg                                 >  2:547957/1‑76 (MQ=21)                                        tttCGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTaa                           <  2:1074477/76‑1 (MQ=21)                                           cGGTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCg                        >  2:1444112/1‑76 (MQ=21)                                             gTCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCt                      >  2:287747/1‑76 (MQ=25)                                              tCCATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTg                     <  1:183985/76‑1 (MQ=255)                                               ccATAATAGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTga                    >  2:1532099/1‑76 (MQ=255)                                                   aataGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaaa                >  1:1335480/1‑76 (MQ=255)                                                   aataGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaaa                >  2:317337/1‑76 (MQ=255)                                                   aataGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGaaa                <  2:788444/76‑1 (MQ=255)                                                    ataGCAGAACAAAAGAAAGAGCCAGAAATGACGGGCAGATATTATATCAGCTCTGCTGATTTAACCGCTGAGAAAt               >  2:469062/1‑76 (MQ=255)                                                     taGCAGAACAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAAtt              >  1:897977/1‑76 (MQ=255)                                                            aCAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCcaca       >  2:1446308/1‑76 (MQ=255)                                                            aCAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCcaca       >  2:1305047/1‑76 (MQ=255)                                                             cAAAAGAAAGAGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAg      <  1:280908/76‑1 (MQ=255)                                                                 agaaagaGCAAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGAt  >  2:478058/1‑76 (MQ=255)                                                                          |                                                                   GAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCTGAGAAATTCGCCACAGCGAT  >  REL606/3555822‑3555961

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