Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,645,397	C→T	intergenic (+1191/+69)	LY180_07765 → / ← LY180_07770	transposase/sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,645,397	0	C	T	82.3%	32.6 / ‑4.6	18	intergenic (+1191/+69)	LY180_07765/LY180_07770	transposase/sulfatase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (8/6);  total (12/6)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTT  >  CP006584/1645165‑1645688                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                     ataCAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACCATAGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGACACGTTCGGTACGCGTGAGAGTTTTACAGTCCAGTAATGCGAGGCAAGCAAACGATAAGCTTATGAGACGTTTAATGTGTAATTAGGGGCAGAAATGGTAAAGAGAGTCGTTTAACAGACCAAGTTCGCACATATTGTTGTCTGATTATTTCTTTTTCGCGAACCAATTTTATCATATTACAAGAttgatatacacctt                                                                                                                                                                                                                                     >  2:264077/1‑285 (MQ=1)     aGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAAAGTTAAGCTGTTGAGTCGTTTTAAGGGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAACTATAGAGTTCGAACATCTTGTTGTATGATTATTTATTTTTCCCGAGACCACTTGATCCTATGACAAGATTTGTAttccacttaact                                                                                                                                                                                                                                 >  2:136679/1‑287 (MQ=1)          gTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATCCGTTTGCTGTTATGTTTTGTTCTCGCGTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGTCAAGCCAACGTTAAGCTGTTGAGTAGTTTTAAGTTTAATTCAGGGCAGCAATGGTAAAGAGAGTCGTGTAACATATCGAATTCGCACCACTTGTTGTATGATTATTGATTTTTCGCGAACACATTTGACCAATTGACTAGATGTTTTTCTCCCGTAACATAAtt                                                                                                                                                                                                                            >  2:281756/1‑297 (MQ=11)                            gACGATTACCTAAAAATCGGTCGAGTCGTGTGATGTTATGTTTTGTTCCCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGTGTGAGCGTTTTCCAGTCAAGTAATGCGTGGCAAGCTAACGCTAAGCTGTTGAGTCGTTTTAAGTGTAATTGGGCGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGCATTAGCACATCTAGTTCTCTGATTATTGATTTGTAGCGAAACCATTTGATCATATGACAAGATGTTTATCTACCTTAACTTAATGCtttttacgcaaatcgat                                                                                                                                                                                                          >  2:57676/1‑285 (MQ=1)                                                                    ttttGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAgt                                                                                                                                                                  >  1:188745/1‑298 (MQ=255)                                                                                          aGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGACACGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGTGGCAGCATTGGTAAAAAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTATGATTATTGATGTTTAGCGAAACCATTTGATCATATGTCAAGAGGTTAATCTACCGTAACTTAATCATTTTGCTAAAAATAATTATGTAATTCACAAGGTTAATACCTAAAAGACTGGTAAAATGTTTTTGAAAGATcaacaaaca                                                                                                                                            >  2:193890/1‑295 (MQ=11)                                                                                              tACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTTGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAAt                                                                                                                                        <  1:29080/298‑1 (MQ=255)                                                                                                 ggCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGc                                                                                                                                     >  1:335693/1‑298 (MQ=255)                                                                                                   ccAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCt                                                                                                                                   <  1:317974/298‑1 (MQ=255)                                                                                                                                  cAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCttt                                                                                                   <  1:193890/298‑1 (MQ=255)                                                                                                                                                  aGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGGTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTAGCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGCTAAAGCATAAACGGCTGTACAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTCA‑CTGACGATTAATCTAACATTTTGCGGGGACACGTTCa                                                                                   >  2:154626/1‑298 (MQ=18)                                                                                                                                                      aaggtttatcgtttGAGTCGTGTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAAGTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGag                                                                               <  2:303727/286‑2 (MQ=255)                                                                                                                                                      aaCGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGCGTATCTACCTTAACTTACTGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATACAAGGCTGTACAATTGTTTTGGCCCGATAAACAACGCCAATTGCCTCATCGGTTAACCTGATC‑CTGATAAAACATTCTTCgggtaagagctaactca                                                                               >  2:143276/1‑284 (MQ=17)                                                                                                                                                                                     agtCGGGGAAGTATTGGTAAAGAGATTTGTGTTAAATATCGAGTTCGGATATTTTGTTGTCTGGTTATTGATTTTTGGCTGAACCATTGGATCATATGACAATAGGTGTATCTACCTTAACGTAATGATTTTGATAATAATCATTAGGGGATTCATAAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGCGACACGTTCACGATTTGTGAACAGTGAAGGGGAGTTGTAATATGc                                                <  2:224245/296‑1 (MQ=255)                                                                                                                                                                                        ctggacAGAATTGGTAAAGAGTGTCGGGTAAAATTTCGAGTTCGCACATTTTGTTGTCTGATTCTTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAGAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGt                                             <  2:188745/293‑1 (MQ=255)                                                                                                                                                                                                          gagagTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGACTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑ATGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATCTTTGAACGGTGCAGGGGAGTTGTAATTTGCCGTTCCTGCACCGTGGAATcc                           >  1:29317/1‑298 (MQ=255)                                                                                                                                                                                                                           tCGAGTTCACACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTTTATATATCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAATCCTAAAAGTATGGTCAATTGTTTTGTCACGATCAAAGACGCAACTTGCCTCATCGGTTAA‑CTGAGCGCTGATATAACCCTTTGCGGGGCACCGTTCACTACTTTTGAACAGTGAAACGGAGTTGTCATCGTCCATTACTGTAGCGGGGAATAACATAAGATAAtaacaac          >  2:107256/1‑292 (MQ=16)                                                                                                                                                                                                                                   gCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGttt  >  1:118772/1‑298 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                     ATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTT  >  CP006584/1645165‑1645688

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