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	84.6%	56.5 / 1.2	26	intergenic (+1191/+69)	LY180_07765/LY180_07770	transposase/sulfatase
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base T (10/12);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATG‑TTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGC  >  CP006584/1645161‑1645680                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            cGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATG‑TTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTAc                                                                                                                                                                                                                               >  1:393182/1‑299 (MQ=31)   gTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATG‑TTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTTAGTGTAATTCGGGGCAGAATTGGTAAAGAGGGTCGTGTAAAATATCGAGTTCGCCCATCTTGTTGTCTTATTATTGCTTTTTCTCGAAACCATTTGAACATCTTACAACATGTTTATCTACct                                                                                                                                                                                                                             >  2:505797/1‑299 (MQ=11)           cGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATG‑TTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAGATATCAAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGAAAAGATTTGTAtacccattaaattaa                                                                                                                                                                                                                     >  2:459633/1‑285 (MQ=0)                    ctTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATG‑TTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACCTCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAAATGTGTATCTACCTTAACTTAATTAttttta                                                                                                                                                                                                            >  2:13315/1‑297 (MQ=11)                           agggttctagtttccTATCACTTGGTCGATTC‑TTTTATGTTTTTGTTGTGTTGTAGCTTTTGGTCGAAGTGTCCGTTCAAGGTCGGAAGGGGTAGGAGGTTTACGGTCTAGTAATGCGTGGCAAGCAAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCATAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAAt                                                                                                                                                                                                     <  2:205109/286‑1 (MQ=11)                                                                         ttttGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAgtg                                                                                                                                                        >  1:612133/1‑299 (MQ=255)                                                                               tctcGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTgg                                                                                                                                                  <  1:546361/299‑1 (MQ=255)                                                                                ctcGCCTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGCGCGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGCTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATCGGCAAAGCCAGTCGAGTAAAATATCGATTTCGCACATTTAGTTGTCTGATTATTGACTATTGGCGAAACCATTCGATCATATGACACGATGTGTATCTACCTTACCTCAATCATTTGAATAAAAATCATTAGGGTATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGCTTTGGc                                                                                                                                                 >  2:132458/1‑299 (MQ=14)                                                                                            ggCAGGTTACGGCCAAGTTCGGTAAGATGGAGATTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATTAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCaacaac                                                                                                                                     <  2:357121/299‑1 (MQ=255)                                                                                            ggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCaacaac                                                                                                                                     <  1:634475/299‑1 (MQ=255)                                                                                              cAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGc                                                                                                                                   >  1:195649/1‑299 (MQ=255)                                                                                                      ggCCAAGTGTGGTACAAGTGAGATTTTTACAGTCAAGGAAGGCGTGGCAAGCCAAAGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGTGCAGTATTGGTAAATAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGTAACCATTTGATCATATGACAAGATTTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGcc                                                                                                                           <  2:577538/299‑1 (MQ=17)                                                                                                                      agTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGATAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGAt                                                                                                           >  1:154379/1‑299 (MQ=255)                                                                                                                         gagaGTTTTGCAGTCAAGTAAGGGGTGGCAAGCAAAGGTTTAGCTGTTGGGTCGTTTTAAGTGTAATGCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGAAAATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATAGGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGc                                                                                                        <  2:375585/299‑1 (MQ=18)                                                                                                                                 taaggtataatAAAGCGTGGGACTCCAAGGTTAATCTGGTGAGTCGTTTTAAGGGTAATTCGGGGCCAGATGGTGAAAGAGAGTCGTGTAAAATATCGAGTTCGCAAAGTTTGTTTTCTGATGATTGAGTTGTTGCGAAACCATTGGATCATATGATAAGATGTGTAGCTACCTTAACTTAATGATCTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATaa                                                                                                <  2:550946/289‑1 (MQ=255)                                                                                                                                       cAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCtttt                                                                                          >  1:389250/1‑299 (MQ=255)                                                                                                                                                 gtgtCAAGGCAACGTTAAGCTGTTGAGTCGTTTTAATTGTAATTCGGGGATGAATTGTTAACGGTAGTCGTGTAAAATGTCGAGTTCGTACATTTTGTTGGCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTGAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGacac                                                                                <  2:489639/299‑1 (MQ=21)                                                                                                                                                    gaaagacatgttttaGCTGTTGTGTCGGTTTAAGTGTAAGTCGGTGCTGAATTTGTAAAGCTGTTCGTTTAAAAGATAGAATTCGCACATTTTATGGTCTGATTATTGTTTGTTTGCGATATCATTTGATCATATGACAAGATGTGTATCTACCTTACTGTAATGATTTTTATAAAAATCATTAGGTGATTAATCATGGTAAAGCAAAAAAGGCTGGTCAAGTTTTTTGGCACGATCATCAACGCCAATTGGCTCATTGGTTAACTTATAGCTGATATACCATCTTGCGGGGACACGtt                                                                             <  2:549388/285‑1 (MQ=255)                                                                                                                                                                           gTCGTTTTAAGTGTAATTCGGGGCAGATTTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCCTACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGAGGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATTAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAggg                                                      <  2:40096/299‑1 (MQ=255)                                                                                                                                                                                      tgtAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAACATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGCTTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAgatatagaa                                           >  2:478938/1‑291 (MQ=255)                                                                                                                                                                                        tAATTCGGGGCAGAATTGGGAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATGATTGATTGTTGGCGAAACCATTTGAGCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATg                                         <  2:612973/299‑1 (MQ=255)                                                                                                                                                                                               gggCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTcc                                  >  1:23749/1‑299 (MQ=255)                                                                                                                                                                                                    gAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAg                             >  1:299719/1‑299 (MQ=255)                                                                                                                                                                                                         ggTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTgg                        <  1:22610/299‑1 (MQ=255)                                                                                                                                                                                                                             atatCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCaaa    >  1:75144/1‑299 (MQ=255)                                                                                                                                                                                                                               atCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGc  >  1:516217/1‑299 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            CGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATG‑TTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGC  >  CP006584/1645161‑1645680

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