Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	4,469,025	+A	coding (748/1350 nt)	LY180_21345 →	permase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	4,469,024	1	.	A	89.3%	169.4 / 8.1	57	coding (747/1350 nt)	LY180_21345	permase
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base A (31/19);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TT‑A‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGAT  >  CP006584/4468748‑4469317                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                         cGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTCGCGCCGGGTGATTTCGCGACCTTCCCGGTGATTACGCCACTGGTAGGTCTAGCGGTGATCACTGGCCCGGAACAACTGACAATCCCGCGTGGCATTCTGCTGACCATTCTCGGTATCTCAATTTGCGGTTTGACCTTACATCCTAACGTCCATTGCTCCGGCCTTTGCGCCCTGCCTTCATGGCGCGATCCAAACGGCATCTCCCGGTTTGCAAGCCCGGCC‑A‑TT‑ACGGGgcgcgcgcacccccgcc                                                                                                                                                                                                                                                                                    >  2:626900/1‑288 (MQ=255)                  tGGTTATTAAAAACCCGCTTTATGGGCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTTACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCt                                                                                                                                                                                                                                                                   <  1:565773/299‑1 (MQ=255)                    gTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGTATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTgg                                                                                                                                                                                                                                                                 >  1:595863/1‑299 (MQ=255)                              aCCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGG‑AAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGtgatg                                                                                                                                                                                                                                                      >  1:128434/1‑299 (MQ=255)                                   cTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACgg                                                                                                                                                                                                                                                  >  1:398136/1‑299 (MQ=255)                                        gtggctggcgggtGCGGTTGGTTATTGTGCGGCGTTGCCGGGGATTTTGGCAGGGGTGGGTCGGGCGGTGGTTATCGGCCTGGAAGAACGGAATGTCCCTGGTGGCATTCTGCTGACCCTTATCGGTATCTCCATTGTCGGTTTGATCTTCGATCCTAACGTCCATTGCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTa                                                                                                                                                                                                                                             <  2:465979/287‑1 (MQ=255)                                             tGCCCTTTGCGCTGGGTGCTTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTTGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTTACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGa                                                                                                                                                                                                                                        <  1:417814/299‑1 (MQ=255)                                              gCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGAt                                                                                                                                                                                                                                       <  1:336552/299‑1 (MQ=255)                                               cccGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCCTTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTCCTCCGGCGTTTTCGCCATGCCTCCATTTACCGATCAAAACGCACATTCCCTGATTGGCAGCCTGTAC‑T‑TT‑ACTGGGTGGCCTGAATCCTGAAGTCCTCCAAAGCGCCAGTGCCGTGGCGATGCCGGATGTATTTGATg                                                                                                                                                                                                                                      >  2:156351/1‑299 (MQ=255)                                                ccGTGGCGCTGGTTGATTTCGCGATCTTCCCGGTGATTATGTGACTGGTTGGTGTGGCTGTGCTCCTCGGCCTGGAAAAACTGTAGGTCCCTGGTGGCTTTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGc                                                                                                                                                                                                                                     <  2:390253/299‑1 (MQ=255)                                                      cgcgggggtgtttCGCGTCCTTCGCGGTGATTGTGGCACTGTGAGGTCTGGCGTGGATCAACGGCCGTGGAAAACTGAAAGTCCCTGGGTGCATTCTGCTGACCATTATCTGTATCTCCAGTGTCGGTTTGTGCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCgg                                                                                                                                                                                                                               <  2:167432/289‑1 (MQ=255)                                                                       aCCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑Cg                                                                                                                                                                                                             >  1:352313/1‑299 (MQ=255)                                                                              cGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCa                                                                                                                                                                                                      <  1:156351/299‑1 (MQ=255)                                                                                 tGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTCA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGc                                                                                                                                                                                                   >  1:374668/1‑299 (MQ=255)                                                                                           aCTGGTAGGTCTGGCGGTGATCATCGGCCTGGACAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCAGTGAGCGATGAAAACGGCACTTCACTGATTGGCCACCTGGAC‑A‑TTAA‑TGGGCGCGCTGCATCCTGTCGTCCTGCCAAGCGTTCTGGCCCGGGGGATGACGGCAGCACTTGCTGCCCCCCGAACTAACCGTGGCGTGCGCCGGCC‑GGGGGACctgccg                                                                                                                                                                                         >  2:505617/1‑297 (MQ=255)                                                                                             tgttAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTACCTGGTGGCATTCTGCTGACCATGATCGGTATCTCAATTTTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCGGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGa                                                                                                                                                                                       <  1:705957/296‑1 (MQ=255)                                                                                                               tcatcGGCCTGGAAAAACTGAAAGTCCCTGCTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TT‑AGTGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTGCTGGCGCGGTTGATGACGGCGGTATTTGATGCCAACCGACCTATCCGTTCCGT‑CCCCGGCCAGGCCGCCCTGCTGGCGAAAACCGGGCAGatcac                                                                                                                                                                     >  2:149486/1‑298 (MQ=255)                                                                                                                     gCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACgg                                                                                                                                                               <  1:149486/299‑1 (MQ=255)                                                                                                                      ccTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCCGGCGAACCTGATGGATAAAGATGGGCAGATCATCGAAGGt                                                                                                                                                              >  1:244864/1‑299 (MQ=255)                                                                                                                        tGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCt                                                                                                                                                                                                                                                                   <  2:407849/197‑1 (MQ=255)                                                                                                                        tGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCt                                                                                                                                                                                                                                                                   >  1:407849/1‑197 (MQ=255)                                                                                                                                       tCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactg‑accac                                                                                                                                            >  1:389490/1‑299 (MQ=255)                                                                                                                                         ccTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactg‑accactg                                                                                                                                          >  1:559948/1‑299 (MQ=255)                                                                                                                                           tggtggCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGcactg‑accactgac                                                                                                                                        >  1:597662/1‑299 (MQ=255)                                                                                                                                            ggtgtCATTCTGTTGACCATTTTAGGTATCTTATTGGTCGGTGTGATCGTCGTTCCTAGCGTCCATTTCTCTGGCGGTTTGGCCATGCCTTCATTGGTCGATGAGAACGGGAATGCGCTGAGTGTGAGCGTGGACTG‑TA‑A‑TGGGCGGGCTGAAGCCTGTAGTCCTGCCAAGCGTTCTGGCGCGGGTGAGGACGGCAGTATTTGATGCCACCGGAACTATCCGCGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACt                                                                                                                                       <  2:139829/299‑1 (MQ=255)                                                                                                                                             gtggCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTc                                                                                                                                      >  1:602728/1‑299 (MQ=255)                                                                                                                                                            tgtttttcgttttctcTAGTGTGGTTTTGTTCGTGGATTCTTGCGGCCGTTTTTGCGGTGTTTTGTTCAGGCCTCTATTGGGCTATGAAACCGGCAATTCAGTGAGTGGCAGCGTGGGC‑AGTA‑A‑TGGGCGCGCTGAGTACTGTAGTCATGCCAAGCGTGCTGGCGCTGGTGATTACGGCAGTATTTGATGCCAACGGAACTATCCGTGCCGT‑CGCCGGCCCGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGtgtt                                                                                                                       <  2:502436/287‑1 (MQ=255)                                                                                                                                                                   cGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGCAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATGTGATGCCACCGGAACTATCCGTGCCTG‑CGCCGGCCCGGCGAACCTGCCGGATAAAGATGGGCAGATCATCGCAGGTGGGACAGCACGG‑ACCACCAACACAATGACAAGCGTTTTCTCCGGc                                                                                                                >  2:25119/1‑299 (MQ=255)                                                                                                                                                                           cAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTggtggg                                                                                                        >  1:665999/1‑299 (MQ=255)                                                                                                                                                                               tGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCg                                                                                                    <  1:716416/299‑1 (MQ=255)                                                                                                                                                                                gTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGAACATAGACGGGGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGAGTTTTCTCTGGCCTGGTGGGTGCgg                                                                                                   >  1:284250/1‑299 (MQ=255)                                                                                                                                                                                 tCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTCGTCCTGCCCAGCGTTCTGGCCCTGGGGATGACGGCAGTATTGGATGCCCACGGAACTATCCGTCCCGT‑CGCCCGCCAGGCGAACCTGCTGGATAAACATGGGCAGGACATCGACGGTGGGACAGCCCTG‑CCCACACACCCACTGCGCAGCGTTTCTCCTGCCCTggttgggccgca                                                                                                  >  2:698012/1‑293 (MQ=255)                                                                                                                                                                                  cGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCt                                                                                                 <  1:161124/299‑1 (MQ=255)                                                                                                                                                                                              cGATCCAAACGTCCTTTCTTCTGGCGTTTTCCCGATGGCTTGATTGCGCGATGAAACCGTGAATCTACTGGTTGGCAGCCTGGAC‑A‑TTAA‑TGGGGGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGGGAGGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCGAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGta                                                                                     <  2:595863/299‑1 (MQ=255)                                                                                                                                                                                               gatctcagcGTCCATTTCGCCGGCGTGTTCGCCAGTCCTTCATTTGACGATGAAAACGGAAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGtat                                                                                    <  2:352313/291‑1 (MQ=255)                                                                                                                                                                                                    tAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGCGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGTGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGACACTATCACTGCCGT‑CGCCTGCCCGGCGAACCGGCCGGCTAACGATCGGCACATCATCCACGGACGGAAAGCACTGCACATCTG‑CCACATTACCAGCGTTTCTCCTGGCCTGGTGGTCGCCGCCCCGGCacgctgatacagc                                                                               >  2:279983/1‑287 (MQ=255)                                                                                                                                                                                                       cGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAg                                                                            >  1:198700/1‑299 (MQ=255)                                                                                                                                                                                                        gTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAACCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGt                                                                           <  1:289627/299‑1 (MQ=255)                                                                                                                                                                                                         tCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTCTTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACCCCATGAGCAGCGTTTTCTCTGGCCTGGTGGTGGGGGCCCCGCCCGCCGCATACACCCAGTc                                                                          >  2:464901/1‑299 (MQ=255)                                                                                                                                                                                                         tCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGACAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGACCTATCCGTGCCGT‑CGCCGGCCAGGCGCACCTGCTGGATAAAGATGGGCAGATCATCGACGTGGGGAAAGCACTG‑ACCACTGACGCCATGCGCAGCGCTTTCTCTGGCGCGGTGGGTGGGGCCCCGGCACCGGtctcccctgtgta                                                                          >  2:321634/1‑287 (MQ=255)                                                                                                                                                                                                          ccATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGCTGCGACCGTAACTATCCGTGCTGT‑CGCCGGCCAGGCGAACCTGCTGGATCAAGTGGGGCAGATCATTGACGGTGGGCAAGCAATG‑ACCACTGACTCCATGAGAAGCGCTTTCTCTGGCCTGGTGGGTGCGGCTCCGGACGCAGtatcctccgcggtg                                                                         >  2:315191/1‑288 (MQ=255)                                                                                                                                                                                                                ctcCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTgcggcg                                                                   >  1:684702/1‑299 (MQ=255)                                                                                                                                                                                                                ctcCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACACCGAGTCTgcggcg                                                                   >  1:97667/1‑299 (MQ=255)                                                                                                                                                                                                                     gCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAc                                                              <  1:25119/299‑1 (MQ=255)                                                                                                                                                                                                                       gTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAcgg                                                            >  1:109358/1‑299 (MQ=255)                                                                                                                                                                                                                         ttttGCCAGGCCTTCAGTTGGCGATGAAAACGGCAATTCACTGAGTGGCAGGCTGGGC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAcggcg                                                          <  2:109358/299‑1 (MQ=255)                                                                                                                                                                                                                             gCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTAcggcggcgg                                                      >  1:604050/1‑299 (MQ=255)                                                                                                                                                                                                                               cATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCggcgggc                                                    >  1:50865/1‑299 (MQ=255)                                                                                                                                                                                                                                       gatagggCGATGTACGCGGCTAATCACTGATTGGAAGGCTTGAT‑A‑GT‑ACTGGGCGCTCTTAATGCTGTCGTCCTGCCATGGGTTGTGGTGCTGGTGATGGCGGCTGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGAGAAAGAGGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACGCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGGTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAAc                                            <  2:59838/293‑1 (MQ=255)                                                                                                                                                                                                                                                aTGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGAc                                   >  1:374295/1‑299 (MQ=255)                                                                                                                                                                                                                                                          cAATTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCAcc                         >  1:306557/1‑299 (MQ=255)                                                                                                                                                                                                                                                           aaTTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTCAAACCGGTTTGACGGCTATCACCg                        >  1:239997/1‑299 (MQ=255)                                                                                                                                                                                                                                                            aTTCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGt                       >  1:326547/1‑299 (MQ=255)                                                                                                                                                                                                                                                             ttCACTGATTGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAACGCACTG‑ACCACTGACTCCATGAACATCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACACCGAGTCTGCGGCGGGTCCGTCGGCGGGCGTTGCAACCAGTTGTACGGcgttcccgggt                      >  2:686808/1‑289 (MQ=255)                                                                                                                                                                                                                                                                     ttGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCt              <  1:400022/299‑1 (MQ=255)                                                                                                                                                                                                                                                                      tGGCAGCCTGGAC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTg             <  1:505617/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                 aC‑A‑TTAA‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTAGCACCGTTGGCGTGCCGTTCCTCCTGAt  >  1:453628/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                         CGCTAACGGTGTCGGTCTGGTGATTAAAAACCCGCTTGATGGTCTGCCCGTTGCGCTGGGTGATTTCGCGACCTTCCCGGTGATTATGTCACTGGTAGGTCTGGCGGTGATCATCGGCCTGGAAAAACTGAAAGTCCCTGGTGGCATTCTGCTGACCATTATCGGTATCTCAATTGTCGGTTTGATCTTCGATCCTAACGTCCATTTCTCCGGCGTTTTCGCCATGCCTTCATTGAGCGATGAAAACGGCAATTCACTGATTGGCAGCCTGGAC‑A‑TT‑A‑TGGGCGCGCTGAATCCTGTAGTCCTGCCAAGCGTTCTGGCGCTGGTGATGACGGCAGTATTTGATGCCACCGGAACTATCCGTGCCGT‑CGCCGGCCAGGCGAACCTGCTGGATAAAGATGGGCAGATCATCGACGGTGGGAAAGCACTG‑ACCACTGACTCCATGAGCAGCGTTTTCTCTGGCCTGGTGGGTGCGGCTCCGGCAGCGGTATACATCGAGTCTGCGGCGGGTACGGCGGCGGGCGGTAAAACCGGTTTGACGGCTATCACCGTTGGCGTGCTGTTCCTCCTGAT  >  CP006584/4468748‑4469317

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