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	83.7%	96.9 / 2.8	43	intergenic (+1191/+69)	LY180_07765/LY180_07770	transposase/sulfatase
Reads supporting (aligned to +/- strand):  ref base C (7/0);  new base T (13/23);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 2.41e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTT  >  CP006584/1645161‑1645688                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                    cGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGAGTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAACGAAAGTCGTGGAAAATATCGACGTCGGACATCTTGTTGTAGGATTATTGATTTTTCGCGCAACCATTTTCCAATATTCCAAGATATGTATCTaa                                                                                                                                                                                                                                       >  2:282099/1‑298 (MQ=11)           cGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGCTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTGGAGTCGTTTTAAGTGTAATTCGGGGCCGAATTGGTAAAGAGAGTCGTGTAAAATATAGAGTTCGCACATCTTGTTGTCTGAGTATTGATTTCTCGCGAAACCATTTGTTCTTATTACAAGATGTTTTTCTACcttaacttaa                                                                                                                                                                                                                             >  2:148364/1‑299 (MQ=11)             aGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGCTTGATGTTCTGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGGAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTCGTCTGATTATTGCGTTTTAGCGAAACAATTTGATCAAATGACAAGATGTGTATCTACCTTACATTAATg                                                                                                                                                                                                                           >  2:697565/1‑299 (MQ=11)                 gCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTACTGCGTGGCAAGCCAACGGTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGTGGCAGAATTGGTAAAGAAAGTCGGGTAAAAAATGAACTTCTCACATCCTGTTGTCTGTATATTTATTTTTAACGAAACCATTTGATCATATGAAAACAGTTTTAACTCCCTTAACCTAATTAttt                                                                                                                                                                                                                       >  2:186636/1‑299 (MQ=11)                 gCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGCGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCCTATGACACGATGTGTATCTACTTTAAATTAATGAttt                                                                                                                                                                                                                       >  2:124741/1‑299 (MQ=9)                 gCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCTAAACCATTTGATCATATAACAATATGTGTATATCCCTTAACTTAATGtttt                                                                                                                                                                                                                       >  2:578386/1‑299 (MQ=9)                                      tgccTTACAATAGGTTGATTGGTTTGATTGTATGTTTGGTTTACGCTTTGGTTGGCTGGTTACGGTCAGGTGCGGTAAGAGTTGGTGTTTTACGGGAAAGTACGGCGTGGCAAGCCATCGTTAAGGTGTTTAGTCGTATTAAGTGTAATTCGGGGTAGAATTGGTGAAGAGAGTCGTGTAAAAGATCGAGTTCGCAGATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGGGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGag                                                                                                                                                                                                  <  2:305362/297‑2 (MQ=11)                                                       cttCGTTTGATGTTATGTTTTGTTCTCGCTTTGCTTGGCAGTTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTACTTCGTGGCGAGCCAACGTTAAGCTGTTGAGTCGATTTAATTTTAATTCGGGGCAGAATTGGGCAAGAGAGTCGTGTAAAATATCGAGTTCGACCATCCTGTTGCCTGATAATTTATTTTTCCACATACCATGTGATAATATGACAAACTTGCTATCTACCTAAACTTAATGATTTTTAAACAAAACAGTAGGGGAAGCATCAgctgcatgtgt                                                                                                                                                                                 >  2:77254/2‑289 (MQ=11)                                                             tttattttttgttTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAgg                                                                                                                                                                           <  2:277947/290‑1 (MQ=255)                                                                           tGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGttt                                                                                                                                                             <  1:671091/299‑1 (MQ=255)                                                                               ctcGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGc                                                                                                                                                         >  1:618379/1‑299 (MQ=255)                                                                                      tggttggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATca                                                                                                                                                  >  1:564367/1‑299 (MQ=255)                                                                                        gttggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGTCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATcaac                                                                                                                                                <  1:131918/299‑1 (MQ=255)                                                                                           ggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCaacaac                                                                                                                                             <  1:482770/299‑1 (MQ=255)                                                                                           ggCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCaacaac                                                                                                                                             >  1:694269/1‑299 (MQ=255)                                                                                             aaGGTTACGGCAAAGTTTGTGAAAGGTGAGAGTTTTTAAGTTAAGTAAGGGTTGGAAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGc                                                                                                                                           <  2:305750/298‑1 (MQ=17)                                                                                                gTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCaa                                                                                                                                        >  1:649635/1‑299 (MQ=255)                                                                                                         aaGTGCGGTAAGGGTGAGGGTTTTACAGTCAAGTAAGGCGTGGTAAGCCAACGTTAAGCGTGTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCAt                                                                                                                               <  2:578147/299‑1 (MQ=17)                                                                                                                                taatataaaGTAGAGCGAGGCAAGCCATGGTGAGTATGTTGAGGCGTTTTAAGGTTAATTGGGGGAAATATTGTTAAAGGATGTTGTTTAAATTGTCGAGTTCGCTGATTTTGTTGTCGGATTATTAATTTTTGGCTAAACCATTTGATCATATGACAATATGTGTATCTAGCTTAACGTAATGATTTTGATAAAAAGCATTAGGGGATTCATCAGGGTAACGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCGAATTGCCTCATCGGTTAACTGATCGCTGGTATaa                                                                                                        <  2:172291/292‑1 (MQ=255)                                                                                                                                  cAGTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAAcc                                                                                                      <  1:263301/299‑1 (MQ=255)                                                                                                                                    gTCAAGTAAGGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCtt                                                                                                    >  1:628126/1‑299 (MQ=255)                                                                                                                                      aaagaaaagggTGGCAAGCCAGCATTAAGGTGGTGGGGTGTGTTAAGTGTAAATCGGGGCGGAATTGGTAAAGAGAGTCGTGTAACATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCAGTTGATGATATGGCAAGATGTGTAGCTACCTTAACTTAATGATTTGGATAAAAATCATTAGGTGATTCATCAGGGTAAAGCATAAAAGTCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCtgtt                                                                                                  <  2:187036/289‑1 (MQ=18)                                                                                                                                                      aGCGAATGTGAAGCTGTGGAGTCGTTTTGAGTGTAATTCGGGGCAGAATTTGTAAAGAGAGTATGGTAAAATATCGAGTTCGCACATTTTGTTTTCTGATTATTGATTTTTGGCGAAGCCATTTGAGCATATGACAAGATGTGTAACTACCTTGACTTAATGATGTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCGGATATAACCTTTTGCGGGGACACGTTCAc                                                                                  <  2:596762/299‑1 (MQ=255)                                                                                                                                                          aaCGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGAtt                                                                              <  1:234158/299‑1 (MQ=255)                                                                                                                                                           aCGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACACGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCCTAAGGGGATTCCTCAGGGTAAAGCATAAAAGGCTGGTCAAGTGATTTGGCACGATCAAAAACGCCAATCACCTCATCGGTAAACTGCTCGCTGCTATAAACTTTTGCGGGGACACGGTCGCGAgtt                                                                             >  2:205803/1‑299 (MQ=18)                                                                                                                                                                  gCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACa                                                                      >  1:214761/1‑299 (MQ=255)                                                                                                                                                                    tGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTACCTAACTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGt                                                                    <  1:296818/299‑1 (MQ=255)                                                                                                                                                                         aGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAgg                                                               >  1:295417/1‑299 (MQ=255)                                                                                                                                                                             gTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAg                                                           <  1:81694/299‑1 (MQ=255)                                                                                                                                                                                ttAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTg                                                        <  1:4314/299‑1 (MQ=255)                                                                                                                                                                                         aTTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGcc                                               >  1:133868/1‑299 (MQ=255)                                                                                                                                                                                                  aGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACACGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCa                                      >  1:6818/1‑299 (MQ=255)                                                                                                                                                                                                               ggagtgttgtAAAATAGCGAGGTGGCTCATTTTGTTGTGTGATTTTTGATGTTTGGCGAAACAATTTGATCATAGTAAAAGATTTGTATCTTCCTTAAAGTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGGCAAGTGTGTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAAGCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCa                         <  2:694269/292‑1 (MQ=255)                                                                                                                                                                                                                 agTCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATa                       <  1:42937/299‑1 (MQ=255)                                                                                                                                                                                                                   tCGTGTAAAATATCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATaaa                     <  1:564178/299‑1 (MQ=255)                                                                                                                                                                                                                      tgtAAAATATAGATTTCGAACAATTTGTTGTCTTGTTAGTGATTTTTGGCGAAGCCATTTGATCATATGACAAGATGTGTATCTAACTTAACTTAATGATTTTGATAAAAGTCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTAAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACGCGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTa                  <  2:256554/299‑1 (MQ=255)                                                                                                                                                                                                                               tCGAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGcc         <  1:147016/299‑1 (MQ=255)                                                                                                                                                                                                                                cgcgTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATCACCTTTCACGGGAACACGTTCACGAATTTTGAACAGTAAAGGGGAGTTGTAATATGCAGTTCATGCAGCGTGAAACCCCACACATTCATCAATGCCa        >  2:172862/4‑299 (MQ=25)                                                                                                                                                                                                                                 gAGTTCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCAc       <  1:256472/299‑1 (MQ=255)                                                                                                                                                                                                                                   gTTCGCAAATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACgg     <  1:507826/299‑1 (MQ=255)                                                                                                                                                                                                                                     tCGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTTAAGTGTTTTGGCACGATCAACAAAGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGCAACGTTCACTATTTTTGAACGGTGAAGGGGAGTTGTAATATGCCGTTCCTGCATGGTGGAATCCCATAAAGAAATCAAAGCCACGGtt   >  2:163048/1‑299 (MQ=37)                                                                                                                                                                                                                                     tCGCAAATTTTGTTGTCTGATTATTTAATTTGTGCGAAACCATGTGATCAGAAGAAAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGTGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGtt   <  2:284687/299‑1 (MQ=255)                                                                                                                                                                                                                                      cGCACATTTTGTTGTCTGATTATTGATTTTTGGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATCTGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAATCCCCGGttt  >  2:235781/1‑299 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                                                    CGGTATACAGCGAGTCGCTCTTTGTGGAGGTGACGATTACCTAACAATCGGTCGATTCGTTTGATGTTATGTTTTGTTCTCGCTTTGGTTGGCAGGTTACGGCCAAGTTCGGTAAGAGTGAGAGTTTTACAGTCAAGTAATGCGTGGCAAGCCAACGTTAAGCTGTTGAGTCGTTTTAAGTGTAATTCGGGGCAGAATTGGTAAAGAGAGTCGTGTAAAATATCGAGTTCGCACATCTTGTTGTCTGATTATTGATTTTTCGCGAAACCATTTGATCATATGACAAGATGTGTATCTACCTTAACTTAATGATTTTGATAAAAATCATTAGGGGATTCATCAGGGTAAAGCATAAAAGGCTGGTCAAGTGTTTTGGCACGATCAACAACGCCAATTGCCTCATCGGTTAACTGATCGCTGATATAACCTTTTGCGGGGACACGTTCACGATTTTTGAACAGTGAAGGGGAGTTGTAATATGCCGTTCCTGCAGCGTGGAATCCCATAAAGTAATCAAAGCCACGGTTT  >  CP006584/1645161‑1645688

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