Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,563,237	T→C	V173A (GTT→GCT)	LY180_07445 →	transposase IS609

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,563,237	0	T	C	76.5%	13.7 / 0.9	17	V173A (GTT→GCT)	LY180_07445	transposase IS609
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (7/6);  total (9/8)
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

GAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAG  >  CP006584/1563086‑1563529                                                                                                                                                        |                                                                                                                                                                                                                                                                                                     aaaatccgCCATATTTTTCGGCCGAACACGGGCCGAAGTCGCTACCGGAATAGCCGTCAGGTGACGGGTGTTGTAAAAAATGTACCTGTCAGCCAGGCCTGCGGTACGTGGTACATCATGATTCAGACAGAAAGTGAAGTATCCCCTCCGGTTCACCCTCCAGCATCAATGGTCGGGCTGGCTGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACACTCTTTGAGCCTGTAAACAGGTTTCAGAAAAACCAGAAGACGCTGGCTAGACTTCAGCGACAGTTAAGCCGCAAgg                                                                                                                                                    <  2:202809/291‑1 (MQ=11)  aaaaCGGCAGTATTTTGCTGCCGAGACTGGGCTTGATGCTCGACCGGACCAGCCGTCATGTCAGGGGTGTGGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGGGGTACATCATTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGt                                                                                                                                                   <  2:49418/298‑1 (MQ=11)                  tCTGCCGAAACTGGGCTGGATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCaacaac                                                                                                                                   <  1:2716/298‑1 (MQ=14)                                  tGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGaaaa                                                                                                                   >  1:114966/1‑297 (MQ=14)                                   ggATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACg                                                                                                                  >  1:201404/1‑298 (MQ=14)                                    gATGCGCTACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGc                                                                                                                 >  1:165474/1‑298 (MQ=14)                                           tACCGGAACAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCTGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGACCCTGTAAACAGTTTTCCGAAAAACAAGAAGACCCTGGCGAGACTTCAGCGCCAGTCAAGCCGCAAGGACAAATTCACCAACAACGGGCAGAAGCAAAAACGCAAAAAAc                                                                                                          >  2:31922/1‑298 (MQ=11)                                                  aTAGCCGTCAGGTCAAGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTAGCCACGCTGTCAGATGGAACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTGAGAGACAGTTAAGCAGAAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAAAACAGCGaat                                                                                                   >  1:35203/1‑296 (MQ=9)                                                     gCCGTTAGGTGACGGGGGGTGTGAAAAATGGCAGTGTGAGGCAGTCCTGCGGGAAGGGGTACATGAGTATTCAGACAGAAAGTGAAGTGTGCACTCGGGTTCACCCTTCATCATCAATGGGCGGGCTGGATGCTGGCGTGGCTAAACTCGCGACGCGGTCAGGTGGCACAGGCTTTGAGCCTGTAAGCAGTTTTCAGAAAAACCAGCAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGAGTGCa                                                                                                <  2:327356/298‑1 (MQ=11)                                                      ccGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCAt                                                                                               >  1:244737/1‑298 (MQ=14)                                                           aGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTg                                                                                          <  1:136141/298‑1 (MQ=14)                                                            ggTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGt                                                                                         <  1:312781/298‑1 (MQ=17)                                                                        tgtgAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGCCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGACACAGTCTTTGAGCCTGTAAACAGTTGTCATAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGAAAAGTCAAATTCACCAAAAACGGGAAGCAGCCAAAACACACAAAACAGCCACCGCATTCCTAAATCACaaaaaac                                                                             >  2:262829/1‑294 (MQ=11)                                                                                      cTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCTGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCt                                                               <  1:31922/298‑1 (MQ=17)                                                                                      cTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCt                                                               <  1:155505/298‑1 (MQ=17)                                                                                                       ggTAAGTGGTACATCAGTATTCAGACAGAACGTGAAGTATCCACTCCGGCTCACCCTTCAGCATCAATGGTCGGACTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCATGTAACCAGCGTTCAGAAAATCCAGAAGACGATGGCGAGACTTCTGCGACAGTTAGGCCGCAAGATCACATTCAGCAACAACTGGCAGAGGCAGAAACACAAACAAAAGCTGCTGAATTACTGTACCGCAAAGAACCGAAGAGGCAACCTCCAAAAAACAACAACgc                                              >  2:20646/1‑297 (MQ=11)                                                                                                                                                   tCCGGCTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAg  >  1:206599/1‑298 (MQ=17)                                                                                                                                                        |                                                                                                                                                                                                                                                                                                     GAAAACAGCCGTATTTTTCTGCCGAAACTGGGCTGGATGCGCTACCGGAATAGCCGTCAGGTCACGGGTGTTGTGAAAAATGTCACTGTCAGCCAGTCCTGCGGTAAGTGGTACATCAGTATTCAGACAGAAAGTGAAGTATCCACTCCGGTTCACCCTTCAGCATCAATGGTCGGGCTGGATGCTGGCGTGGCTAAACTCGCCACGCTGTCAGATGGCACAGTCTTTGAGCCTGTAAACAGTTTTCAGAAAAACCAGAAGACGCTGGCGAGACTTCAGCGACAGTTAAGCCGCAAGGTCAAATTCAGCAACAACTGGCAGAAGCAGAAACGCAAAATACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGACCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAG  >  CP006584/1563086‑1563529

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