Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,645,421	C→G	intergenic (+1215/+45)	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,421	0	C	G	81.2%	29.2 / ‑3.9	19	intergenic (+1215/+45)	LY180_07765/LY180_07770	transposase/sulfatase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base G (8/5);  total (13/6)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGC  >  CP006584/1645165‑1645708                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                 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)                                                                                                                                                                                                                                                       attGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCAAAAAGGAATCAAAGCCACGGTTTTTAGGTTGACATTCTTCCGc  >  2:67050/1‑298 (MQ=255)                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                 ATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAA‑CTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTTTGAGGTTGCCATTCTTCCGC  >  CP006584/1645165‑1645708

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