Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,645,301	T→G	intergenic (+1095/+165)	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,301	0	T	G	79.2%	57.2 / 8.7	24	intergenic (+1095/+165)	LY180_07765/LY180_07770	transposase/sulfatase
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base G (8/11);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 4.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

CGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTT  >  CP006584/1645161‑1645592                                                                                                                                             |                                                                                                                                                                                                                                                                                                     cGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTTTTTATTTTTCGCGAAACCATTTGATCATTTGAccaaaatgttatcca                                                                                                                                         >  2:570310/1‑284 (MQ=0)    tataCAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTAGGTTGGCAGGTTACGGACAAGTTCGGCAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGATGACTCGCTTTAAGTGTAATTCTGGGCAGACTTGGTAAAGAGAGTCGTGTAAAATATCGGATTTGCAAGTCTGGTTGTCTGATGATTGCTTTTTGGCGGAACCATTTGATCATATGACAAGAcgttttaccacct                                                                                                                                      >  2:142881/1‑288 (MQ=1)         aGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTGTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGGTCAAAATATCGAGTTCGCACATCTTGTTGGCTGATTAGTGTTTTTTCGCGAAACCTATTGACCATTTGACAAGCTgtgttacccccttcact                                                                                                                                 >  2:807286/1‑285 (MQ=1)          gCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCTTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTCAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTAGTGTAAAATATAGAGTTCGCACATCTTGTTGTCAGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATTTTTTTCTACcttaactt                                                                                                                                >  2:237133/1‑298 (MQ=9)                               tGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCCGGTTACGTCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTAGTTTTAAGTGTACTTCGGGGCAGAATTGATAACGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACAATTTGTACATATTACCAGAGTTGTATCTAACTTAACTTAATTATTTTTAATAAAAAtat                                                                                                           >  2:425356/1‑296 (MQ=11)                                                        ttCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCAt                                                                                 >  1:433573/1‑298 (MQ=21)                                                                    tatgttTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTc                                                                     <  1:108146/298‑1 (MQ=255)                                                                                tcGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGc                                                         <  1:682812/298‑1 (MQ=255)                                                                                tcGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTG‑TGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCa                                                        >  1:650442/1‑298 (MQ=255)                                                                                  gCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCAc                                                       <  1:474059/298‑1 (MQ=255)                                                                                   cTTTGGCTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAACCATAAAAGGCTGGTCAAGTGTTTTGGCACg                                                      >  1:47614/1‑298 (MQ=255)                                                                                     ttggttggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGAt                                                    >  1:685253/1‑298 (MQ=255)                                                                                         gtggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATcaac                                                <  2:209349/297‑1 (MQ=255)                                                                                               ggTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGcc                                          <  1:77704/298‑1 (MQ=255)                                                                                                                gaaagaggaggagTTTTACAGTAAAGTAAGGGGTGGCAAGACAACTTTAGGCGTTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCTTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTa                         <  2:936/288‑1 (MQ=17)                                                                                                                gTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTGAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTa                         >  1:638036/1‑298 (MQ=255)                                                                                                                     agTGAGAGGTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGGGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGa                    <  2:590400/298‑1 (MQ=255)                                                                                                                          gagTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCAt                                                                                 <  1:90172/232‑1 (MQ=25)                                                                                                                          gagTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCAt                                                                                 >  2:90172/1‑232 (MQ=25)                                                                                                                             ttttACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGat            >  1:26893/1‑298 (MQ=255)                                                                                                                                 aCAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCGATTGCCTCATCGGTTAACTGATCGCTGATATaa        >  1:212996/1‑298 (MQ=255)                                                                                                                                    gTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCt     <  1:699804/298‑1 (MQ=255)                                                                                                                                     gcAAGTAAGGGGTGGCAAGCCAACGTTGAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGAAGAATTGGTAAAGAGAGTCTGGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCtt    <  2:565246/297‑1 (MQ=255)                                                                                                                                       aaGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCtttt  <  1:775190/298‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                                                                                                                                                                     CGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGA‑TTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTT  >  CP006584/1645161‑1645592

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