Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	4,305,151	G→C	P51P (CCG→CCC)	LY180_20645 →	rhsA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	4,305,151	0	G	C	100.0%	132.8 / NA	45	P51P (CCG→CCC)	LY180_20645	rhsA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (23/22);  total (23/22)

GCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAG  >  CP006584/4304859‑4305313                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                    gcgcCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCcctgct                                                                                                                                                                <  1:346663/298‑1 (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:659144/1‑298 (MQ=255)                                                                  gtTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTG‑CCGTTCATCctct                                                                                              <  1:355264/298‑1 (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)                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                    GCGCCGCACGGCGCTTTTTTTGTGCCTGTAATTTGCCTTTTATAAATTCTGAAAAATTCCGGGGTGTTAACTCATTATCATTTCCTGTTG‑TTGACTGGCGCTAGGCTTCTGGTTCTTATTTATTTCAAAAGGATTTTACGATGAGCGGAAAACCGGCGGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTCCAGGGTTCAGCCGGGGTACGTATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGCCCCGCTG‑CCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAG  >  CP006584/4304859‑4305313

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