Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	1,168,401	C→A	S606* (TCG→TAG)	LY180_05445 →	glucosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	1,168,401	0	C	A	94.7%	181.1 / NA	59	S606* (TCG→TAG)	LY180_05445	glucosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (26/28);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTC‑TCCGCTCCGCTGT‑GGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATCGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCACCGTTTGTTTCGGTGATTTCCAGC  >  CP006584/1168106‑1168690                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                   ccGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGCCTTATCTC‑TCCGCTCCGCTGT‑GGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGt                                                                                                                                                                                                                                                                                                  <  1:140664/298‑1 (MQ=255)         ccGTGCGGGTTTCCTGACGGGCGTGATTTCTTATCTC‑GCCGCTCCGCTGT‑GGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTGCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGctgctg                                                                                                                                                                                                                                                                                          <  2:322923/298‑1 (MQ=255)                    ttcttgcGGGCTTGATTTTTTATCTC‑TCCCCT‑CGCTGTGGGTTTATGTTCCTCGCGCTCTGTTCTTGATTGCAGGTAGTACTTGCGTTGACCGAACCGCAAGACTTCCTGCAACCACGGCAGTTGTTCCGGGTGTGGCCGCAGTGGCGTCCTGAGCTGGGGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCtttt                                                                                                                                                                                                                                                                               <  2:163196/292‑1 (MQ=255)                           gggCGTGATGTCTTATCTC‑TCCGCTCCGCTGT‑GGTTTATGTTCCTCGCGCTCTCTACTGCATTCCAGGTAGTACATGCGTTGACCTAACCGCAATACTTCCTTCACCCACGCAAGTTTTTCCCGGTGTGGCCGCAGTCGCGTCCTGAGCTGCCGATTGCCCCTCTTGCTTCCATCATGGTCCTTTTGTTCCTGCCGACGGTCTTTAGCATTTTGCTTATCTGTTGCAATGCAACCAAAGAATCCGGCGTCCTCTGCCTCGGTAGATTACCGTTGCTGCAGGAGAGGCTTGTTTcctgt                                                                                                                                                                                                                                                                        >  2:603021/1‑295 (MQ=255)                             gCGTGCTGTCTT‑TCTGTTCCGGTGCGCTGT‑GTTTTATGTCCCTCGCTCGTTGTGCTTGATTGCAGTTAGTGACTGGGTGGACGCAACCGCAATACTGCCTGGAACCACGTCAGTTGTTCCCGGTGTGGCTGCAGTGGCGGCCTCATCTGGCGATTGCACTTTTTGCTTCTACTATGGTGCGGTTGTTCCTGCCGAAGTTATTGCGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGtgtt                                                                                                                                                                                                                                                                      <  2:265684/298‑3 (MQ=255)                                   tGTCTTATCTC‑TCCGCTCCGCTGT‑GGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAAACGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTTCTGCCGAAGTTATTGAGCATTTTGATTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGGGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGc                                                                                                                                                                                                                                                                >  1:187452/1‑298 (MQ=255)                                          tctc‑tcCGCTCCGCTGT‑GGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTa                                                                                                                                                                                                                                                         >  1:630859/1‑298 (MQ=255)                                          tctc‑tcCGCTCCGCTGT‑GGTTTATGTTCCTCGCGCTATCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCTCGTTACATTATAGTTGCCGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTa                                                                                                                                                                                                                                                         >  1:359495/1‑298 (MQ=255)                                                cggtgccgcgtt‑gGTTTTTGTTCCTCGCGCGCTCTGCTGCATTGCTGGTGGTGCATGCGGTGATCGAACCGGAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGGGCTGGCGATTGCACTTTTTGCTTCGAGTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCAt                                                                                                                                                                                                                                                    <  2:8730/287‑1 (MQ=255)                                                      cgctGT‑GGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGTGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGtt                                                                                                                                                                                                                                              >  1:171119/1‑298 (MQ=255)                                                         tgt‑gGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCGGTGCTGCTGGCTCCGGTACGCATGCTGTTcaa                                                                                                                                                                                                                                           >  1:560330/1‑296 (MQ=255)                                                         tgt‑gGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCa                                                                                                                                                                                                                                           >  1:502071/1‑298 (MQ=255)                                                          gt‑gGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCAt                                                                                                                                                                                                                                          <  1:492107/298‑1 (MQ=255)                                                                  aTGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCt                                                                                                                                                                                                                                   >  1:507187/1‑298 (MQ=255)                                                                        cTCGCGCTCTCTACTTCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTg                                                                                                                                                                                                                             <  1:498772/298‑1 (MQ=255)                                                                           gcgcTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCa                                                                                                                                                                                                                          <  1:237380/298‑1 (MQ=255)                                                                                 tctACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGt                                                                                                                                                                                                                               >  1:296288/1‑287 (MQ=255)                                                                                 tctACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGt                                                                                                                                                                                                                    >  1:159313/1‑298 (MQ=255)                                                                                 cctACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCCCGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGt                                                                                                                                                                                                                               <  2:296288/286‑1 (MQ=255)                                                                                    aCTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGAGCGttac                                                                                                                                                                                                                 >  1:507878/1‑296 (MQ=255)                                                                                     cTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAATTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCt                                                                                                                                                                                                                <  1:549178/298‑1 (MQ=255)                                                                                        cATTGCGGGTAGTACACTCGTTGTGCGAACCTCCAGACTTCCTGCAAGCACGGGATTTGTTCCCGATGTGGCCGCAGTGTCGTCCTTGGCTGGCGATTGCATTTTTTGCTTCGACTAGGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTCGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTgg                                                                                                                                                                                                             <  2:635786/298‑1 (MQ=255)                                                                                         aTTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGATTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGc                                                                                                                                                                                                            <  1:653676/298‑1 (MQ=255)                                                                                                                   aaCCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCa                                                                                                                                                                                  >  1:350799/1‑298 (MQ=255)                                                                                                                             cTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGtgatgat                                                                                                                                                                        >  1:396936/1‑298 (MQ=255)                                                                                                                                     gcctccgtcAGTTGTCCGCGGTGGTGCCGCAGTGGTGGTCTGGGCGGGGGATTGCATTTTTTGCTTCGAGTATGGTGCTGGTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCGGGTGAAAAGGAACGAAAGGATACGGCGGCTTCTGGCGCGTTACATTATAGGTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGAGCGTTCCTTGGCTGGGAAGTTGTGTGGAATTCACCGCAGCGTGATGATGACTCCAc                                                                                                                                                                <  2:410535/290‑1 (MQ=255)                                                                                                                                                   ttCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGCAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGCATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGACGTGCTTTTTTCCGTGCTTCTGGCTCCGGTACGCATTCTGTTCCATACGGTCTTCGTTGTCAGCGCCGTCCTTGGCCGGGAAGTGGGGTGGAATTCCCCGCAGCCTGGTGCTGACCCCCCCTCCTcggggggag                                                                                                                                                  >  2:629589/1‑293 (MQ=255)                                                                                                                                                           gtgGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCaaa                                                                                                                                         >  1:297002/1‑298 (MQ=255)                                                                                                                                                               ccGCAGTGGCGTCCTGTGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTTTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTTCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGcc                                                                                                                                     <  1:273196/298‑1 (MQ=255)                                                                                                                                                                 gCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATCCGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCCGCGTGATGATGACCCCACTTCCTGGGG‑GGACGCGTTCAAGCGccgc                                                                                                                                   >  2:24797/1‑296 (MQ=255)                                                                                                                                                                    gTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGc                                                                                                                                <  1:340670/298‑1 (MQ=255)                                                                                                                                                                     tGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTTGCTGGGAAGTGGTGGGGAATTCACCGCCGCGTGATGATGACTCCCCTTCCTGGGG‑GGAACCGTCCAAACGCCAAGGCt                                                                                                                               >  2:434904/1‑298 (MQ=255)                                                                                                                                                                                gCTGGGGATTTGCCTTTTTGGTTGGACTTTGGGGGTGTTTTTGCTGCCGGTGTTGTTGAGTATTTTGGTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGTGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTGCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTGCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCCCAGctgctg                                                                                                                    <  2:507187/298‑1 (MQ=255)                                                                                                                                                                                 cgggggtttgtccTTTTGGCTTCGCCTATGGTGCTGTTGTTCCTGCCGAATTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTTTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGt                                                                                                                   <  2:630859/286‑1 (MQ=255)                                                                                                                                                                                   ggCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTa                                                                                                                 >  1:250688/1‑298 (MQ=255)                                                                                                                                                                                        ttGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGtt                                                                                                            <  1:1989/298‑1 (MQ=255)                                                                                                                                                                                               ttttGCTTCGGCGTTGGGGCTTTTGTTCCTGCGTAATTTATTTGGCGTTTTGGTTTGCTGGGGCAAAGGAACGAAAGGATAGGGCGGCTTCTGGGGCGTTACATTATAGTTGCTGGTGGAAGTGCTTTTTTCCGTGCTGGTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTgg                                                                                                     <  2:159313/298‑1 (MQ=255)                                                                                                                                                                                                tttGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTggg                                                                                                    <  1:625805/298‑1 (MQ=255)                                                                                                                                                                                                    cTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGC‑TTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGATTTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCGGTTAGGGTTAGTGTGGGCTGtt                                                                                               >  1:225618/1‑298 (MQ=255)                                                                                                                                                                                                     ttCGACTGTGGGGCTGTTTGTCCTGCCGACGTTTTTGAGCTTTTTTCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTCGCGCTTTACAGTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGGTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGGGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGtt                                                                                               <  2:171119/298‑1 (MQ=255)                                                                                                                                                                                                     ttCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGtt                                                                                               >  1:233666/1‑298 (MQ=255)                                                                                                                                                                                                             tGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATCGTTGCTGCTGGAAGTGCTTCTTTCCGCGCTGCTGGCTCCGGTACGCATGCTGTTCAATACTGTCCTCGTTGCCCGCCCGTTCCTCGGCTGGGAATACGTTTGTACTTCACCGCACCCCGTTGCTG‑CCCCACCTCCCGGGGTTGAAGCGTTAAAACGCCACGACGCACAGCACCTGCTAGGGCTCGGGTGTGCTGTGGTGATgca                                                                                       >  2:43440/1‑296 (MQ=255)                                                                                                                                                                                                                  cTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCGGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATCGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCGGCGCGTTCCTTGGCCGGGGAGGTGTGTGCAATTCACCTGAGCGTGGTTATTACCCCACTTCGGGGGG‑TGTGGCGTTCAACCGCCACGGCGCACAGCTGCTGCTAGGGCTGGTGTGGCGTGGTGGGGCTGCGTggg                                                                                  >  2:140556/1‑297 (MQ=255)                                                                                                                                                                                                                    gttgttCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCGGGGG‑TGAACCTTTCAAACGCCACGTCCACCAGCCGCTGGTAGGGTAAATGGGGGCTGTCGGGaagccgggcctg                                                                                >  2:268166/1‑287 (MQ=255)                                                                                                                                                                                                                          ctggccagaGTTTTTGGGCATTTTTCTTATCTGGTGCAAGGGACCGCAAGAATGTGGGGGGTTCTGGGGCGTTACTTTATGGTGGCTGCTGTAAGTGCTTTTTTCCGTGGTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGGGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTg                                                                          <  2:502071/290‑1 (MQ=255)                                                                                                                                                                                                                                  aGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCt                                                                  >  1:96522/1‑298 (MQ=255)                                                                                                                                                                                                                                       tttgtcATTTTGCTTATCTGGGGAAAATGCACGAATGAATGCGGCGGCTTTTGGGGTGTTAAATTATTGGTGCTGTTGGAAGTGCTTGTTTGCGTGCTGGTTGGTCGTGTACGCAGGCTTTTCCGTACGGTTTTCGTTGTCAGCGCGTTCCTTGGCTTGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCGGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTct                                                             <  2:396936/293‑1 (MQ=255)                                                                                                                                                                                                                                              ttttGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCtggctggc                                                      <  1:277320/298‑1 (MQ=255)                                                                                                                                                                                                                                                   cTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCTCAGCGGGCTGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCCCAGCTGCTGTTAGGGTTAGTGGGGGCTGTTGGGATGGCGTGGCTGGACCTGCTTTTCCTGGTCTGGCCGGCAccaa                                                 >  2:656184/1‑296 (MQ=255)                                                                                                                                                                                                                                                   cTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGa                                                 <  1:140556/298‑1 (MQ=255)                                                                                                                                                                                                                                                   cTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCGGTTCTGGCTGGCACCGa                                                 >  1:337040/1‑298 (MQ=255)                                                                                                                                                                                                                                                      aTCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTg                                              >  1:12654/1‑298 (MQ=255)                                                                                                                                                                                                                                                               aaaGGAACGAAAGAATACGGCGGCTTCAGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGCACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGt                                     <  1:613586/298‑1 (MQ=255)                                                                                                                                                                                                                                                                           gAATACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCAc                         <  1:395385/298‑1 (MQ=255)                                                                                                                                                                                                                                                                            aaTACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCAcc                        <  1:656184/298‑1 (MQ=255)                                                                                                                                                                                                                                                                             aTACGGCGGCTTCTGGCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGGTAGTGTGGGCTGTTGGGATGGCGTGGCTGGTTCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGCTCCTGTCACCg                       >  1:307986/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                      cTTCTGGCGGCTTACATTATAGTTGCTCCTGGAAGTGCGTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCACCGTTTGTTTCg              <  1:3091/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                           ggCGCGTTACATTATAGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCACCGTTTGTTTCGGTGAt         >  1:270576/1‑298 (MQ=255)                                                                                                                                                                                                                                                                                                  tACATTATAGTTGTGGCTGGAAGTGCTTTTTCCCGTTATGCTGGCTCCGTTACTCATGCTGTTCCTTACGGTCTTTGTTGTTAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGTGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCACCGTTTGTTTCGGTGATTTCCAGc  <  1:497756/298‑1 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                   CCGGTTCACCGTGCGGTGTTCCTGACGGGCGTGATGTCTTATCTC‑TCCGCTCCGCTGT‑GGTTTATGTTCCTCGCGCTCTCTACTGCATTGCAGGTAGTACATGCGTTGACCGAACCGCAATACTTCCTGCAACCACGGCAGTTGTTCCCGGTGTGGCCGCAGTGGCGTCCTGAGCTGGCGATTGCACTTTTTGCTTCGACTATGGTGCTGTTGTTCCTGCCGAAGTTATTGAGCATTTTGCTTATCTGGTGCAAAGGAACGAAAGAATACGGCGGCTTCTGGCGCGTTACATTATCGTTGCTGCTGGAAGTGCTTTTTTCCGTGCTGCTGGCTCCGGTACGCATGCTGTTCCATACGGTCTTCGTTGTCAGCGCGTTCCTTGGCTGGGAAGTGGTGTGGAATTCACCGCAGCGTGATGATGACTCCACTTCCTGGGG‑TGAAGCGTTCAAACGCCACGGCTCACAGCTGCTGTTAGGGTTAGTGTGGGCTGTTGGGATGGCGTGGCTGGATCTGCGTTTCCTGTTCTGGCTGGCACCGATTGTCTTCTCGTTGATCCTGTCACCGTTTGTTTCGGTGATTTCCAGC  >  CP006584/1168106‑1168690

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