Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	4,305,124	G→A	G42G (GGG→GGA)	LY180_20645 →	rhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	4,305,124	0	G	A	96.2%	137.4 / NA	54	G42G (GGG→GGA)	LY180_20645	rhsA
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (23/27);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 2.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAG  >  CP006584/4304831‑4305313                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                               agcaAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGa                                                                                                                                                                                            >  1:454934/1‑298 (MQ=255)   caAATAAGACCTGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACg                                                                                                                                                                                          <  1:68060/298‑1 (MQ=255)     aaTAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTc                                                                                                                                                                                        <  1:290267/298‑1 (MQ=255)         aGACCGGATTTCCGAAAACAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGGGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGACAACCGGCGTCGCGGACGGGCGACATGACGAAGTATGGCGGTAGCATTACCCAGGGTTCACCCGGGGTAGCTATTGGGGCCCCCACCGGCGTGGCCTGGTCGGTGTGCACCGGCGGGGTGGCGTcgtgc                                                                                                                                                                                    >  2:605092/1‑294 (MQ=255)                    gcGAAAAGAGCGCCGAGTGGGGCTTTTTTTGTCCCGGAAATTGGCCTGTTTGAATTTTTGTAAAGTTCCTGGGTGTTATTCCATTTTCATTTCCTGTTG‑TTTACTGGCGCTAGGTTTTTGGTTTTTATTTATTTCAGAAGGATTTTACGAGGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCaa                                                                                                                                                                         <  2:116314/297‑1 (MQ=18)                        cgatctcgccgcACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATccc                                                                                                                                                                     <  1:202597/292‑2 (MQ=255)                          aaaGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAACAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTATGGTTCTTATTTATTTCTAAAGGATTTTACGATGCGCGGACAACCGGCGGCGCGTCAGGGCGACCTGCCGCAGTCTGCCGGGACCATTGTCCAGGGTTCAGCGCGGTTACGTATTGGTGCCCCAACCGCCGTGGCCTGTCTGGGGTGCCCCGGCGGAGTGCCGTCCGGCACTCAGGCCAATcaccg                                                                                                                                                                   >  2:294512/1‑294 (MQ=255)                             gcgcCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCcctgct                                                                                                                                                                <  1:346663/298‑1 (MQ=255)                               gccgcACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGGTCGGTGTGCCCCGGCGGGGTGGCCTCCGGCCCTCCCGTCAAcccccccaccg                                                                                                                                                              >  2:202597/1‑290 (MQ=255)                                    aCGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCTTTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCa                                                                                                                                                         >  1:581277/1‑298 (MQ=255)                                               tttGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTcc                                                                                                                                              <  1:318490/298‑1 (MQ=255)                                                  gtgCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCgg                                                                                                                                           >  1:259385/1‑298 (MQ=255)                                                     catttAATTTTTCTTTTTTTAATTCTGAAAAATTCCGGGGTGTTAACTCAGTATCATTTCCGGTTG‑TTGACTGGCGCTCGGCTTCTGGTTCTTATTTATGTCAAAAGGATTTTACGATGAGCGGAAGAGCGGGGGCGCGTCCGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGa                                                                                                                                        <  2:157785/294‑1 (MQ=18)                                                        gTACTTTGCTTTTTAGAAATTGTGAAAATGTCCGGGGTGTTAATTCTTTTTCGTTTCC‑GTTGTTTGAGTGGCGTTAGGCTTCTGGTTCTGACTTATTTCAAAAGGATTTTACGATGGGCGGAAAACCGGCGGCGCGTAAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAc                                                                                                                                     <  2:122188/298‑1 (MQ=17)                                                         tAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAcc                                                                                                                                    <  1:434329/298‑1 (MQ=255)                                                          aaTTTTGCTTCTATCAATTGTGAAAAATGCCTGGGTGTTAACTCATTCTCATTTCTTGTGT‑TTGAATGGGGCTCGGTGTCTGGTTCTTATGTAATGCAAAAGGATGTGACGATTGGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTAGGGCTGTAGCATTGTCCATGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGGGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCg                                                                                                                                   <  2:359350/298‑1 (MQ=17)                                                               gatttttATAAATCTTAAAAAATTCCGGGGGGTTAACTCTTTCTCTTTTCCGGTTG‑TTGACTGGGGCTAGGCTTCTGGTTCTTATTTATTTAAAAAGGATTGTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATc                                                                                                                              <  2:293888/295‑1 (MQ=18)                                                                 cTTTTGTAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAg                                                                                                                                                      <  2:541478/272‑1 (MQ=255)                                                                 cTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTT‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGc                                                                                                                            <  1:412976/298‑1 (MQ=255)                                                                 cTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTCTTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGAATTGTCCAGGTTTCAGCCGGGGTCCGCACCGGTGCCCCCACCGGCCTGGCCTGTTCGGTGTGCCCCCGCGGCGTGACGCCCGGCCATCCGGCCACTCCCCTGCTCGCTGCACACGTCCTTCCCTGTGAAACCGaacccac                                                                                                                            >  2:611510/1‑292 (MQ=21)                                                                 cTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAg                                                                                                                                                      >  1:541478/1‑272 (MQ=255)                                                                   tttATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGccc                                                                                                                          <  1:454634/298‑1 (MQ=255)                                                                     tataAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTg                                                                                                                        >  1:223099/1‑298 (MQ=255)                                                                              tGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGCTTTTACGATGAGCGGAAACCCGGCGGCGCGTCAGGCCGACATGACGCAGTATGGCGGTAGCATTGTCCCGGGTTCAGCCGGTGACCGTATTCGTGCCCCCACCGGCCTGGCCCGTCCGCTGTGCCCCCCCGGAGTGCAGCCCCGCCATCCGGTCCACCCCCTGCCCTGCACCACGCCCCTCCCCGGCGAAACCAACACCGCCCTgcccggcccg                                                                                                               >  2:5636/1‑298 (MQ=18)                                                                                  aaaTTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑c                                                                                                          <  1:123479/298‑1 (MQ=255)                                                                                        cGGGTTGTTAACTCATTTTCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTTGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAACCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCa                                                                                                    <  1:605092/298‑1 (MQ=255)                                                                                              gtTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCctct                                                                                              <  1:355264/298‑1 (MQ=255)                                                                                              gtTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCctct                                                                                              >  1:659144/1‑298 (MQ=255)                                                                                              gtTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCctct                                                                                              >  1:350306/1‑298 (MQ=255)                                                                                                tacccctttgtCATTTTCTGTTG‑TTGACGGGCGCTTGGCTTGTGTTTCTTATTTAGTTCAAAAGGATTTTCCGGTGGGCGGAAACCCGGGGGGGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCCTTGTCCAGGGTGCAGCCGGGGTACGTATGGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTcc                                                                                            <  2:350306/288‑1 (MQ=14)                                                                                                  aCTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCg                                                                                          >  1:300984/1‑298 (MQ=255)                                                                                                     cATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGCACGTATTGGTGCCCCCACCGGCGTGGCCTGCTCGGTGTGCCCCGGCGGAGTGCCGTACGGCCATCCGGTCCACCCCCTGCTCGGTGCAAAGGTCCTTCCCGGCGCACCCCACCTCGACCTCCCCGGCCCCCAG‑CCCTTCATCTTCCACCGccc                                                                                       >  2:221432/1‑296 (MQ=18)                                                                                                       ttATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGATATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCt                                                                                     >  1:436896/1‑298 (MQ=255)                                                                                                            aTTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGATGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTAcagc                                                                                >  1:324904/1‑298 (MQ=255)                                                                                                                         tgACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGCACGTATTGGTGCCACCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGCGTGACGTCCGGCCATCCGGTCAATCCCCTGCCCGGTGCAAAGGGCCTTCCCGTGGCACACGACATAGCCCCACCCGGCCCGCTGCCCGTTC‑TCCCCCCCCGCCCCCACCCCACTTAcccagct                                                                    >  2:269231/1‑293 (MQ=17)                                                                                                                             tGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACaaaaa                                                                >  1:315570/1‑298 (MQ=255)                                                                                                                              ggCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAAc                                                               <  1:468600/298‑1 (MQ=255)                                                                                                                                     ggCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcgcg                                                        >  1:279958/1‑298 (MQ=255)                                                                                                                                                cTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAg                                             <  1:658613/298‑1 (MQ=255)                                                                                                                                                  tatttattTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGcc                                           >  1:226549/1‑298 (MQ=255)                                                                                                                                                   atttatttCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACCCAAACGCCCGCGCCGGTGGGGAGCCt                                          >  1:428183/1‑298 (MQ=255)                                                                                                                                                       atttCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGGGCAAAGGTCCTTCCCGGGGCAACCGACATCGCCCTGCCCGGCCCGCTG‑CCCTTCATCCTCTCCCGCCCCTACAGCAGCTACCGGCCAAAACCGCCCGCCCCGGTGGGGGCCCCCGGc                                      >  2:612333/1‑298 (MQ=11)                                                                                                                                                         ttCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGGACGTATTGGTGCCCCCACCGGCCTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCCGTCAATCCCCTGCCCGGTGCAAAGCTCCTTCCCGGTGACACCGCCATCGCCCTCCCCGGCCCCCCG‑CCGCCCCCCCCCCCCCGCACCTCCACCGCCTCGCGGCCCAAACCGCCCGCCCCGGCGGCGGGCCCCGGccc                                    >  2:191400/1‑298 (MQ=11)                                                                                                                                                           cAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCg                                  <  1:368217/298‑1 (MQ=25)                                                                                                                                                            aaaaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCgg                                 >  1:483309/1‑298 (MQ=25)                                                                                                                                                             aaaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGTCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGc                                <  1:394347/298‑1 (MQ=25)                                                                                                                                                             aaaGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCCGTTACCGGACAAAAACGCCCGCCCCCGCGGGGAGCCTCGGGCCCGGc                                >  1:35059/1‑298 (MQ=255)                                                                                                                                                                  aTTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCCGCCCGCCG‑CCGTTCATCCTCCCCCGCACCCACAGCAGTCACCGGCCAAAACCGCCCGCCCACGCGGGGCGCCCCGCCCCCgcgccgca                           >  2:269975/1‑291 (MQ=11)                                                                                                                                                                       aCGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGc                      <  2:279958/298‑1 (MQ=14)                                                                                                                                                                         gATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCt                    <  1:425964/298‑1 (MQ=14)                                                                                                                                                                                  aaagCCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATc           <  1:588322/298‑1 (MQ=14)                                                                                                                                                                                         gcggcgCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgc                                                                                                                       <  1:648200/184‑1 (MQ=18)                                                                                                                                                                                         gcggcgCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTgc                                                                                                                       >  2:648200/1‑184 (MQ=18)                                                                                                                                                                                           ggcgCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAg  <  1:612333/298‑1 (MQ=11)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                               AGCAAATAAGACCGGATTTCCGAAAAAAGCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAG  >  CP006584/4304831‑4305313

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