Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006584	4,179,007	C→T	noncoding (1037/1541 nt)	LY180_19945 →	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006584	4,179,007	0	C	T	98.4%	97.3 / NA	63	noncoding (1037/1541 nt)	LY180_19945	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/32);  minor base G (0/1);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTTCAGAGATGAGAAT‑GTGCCTTCGGGAACCGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGACCTCATAAAGTGC  >  CP006584/4178739‑4179262                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                 cGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAg                                                                                                                                                                                                                                    <  1:201991/299‑1 (MQ=9)             aCAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGa                                                                                                                                                                                                                        <  1:575353/299‑1 (MQ=9)                     aGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTgg                                                                                                                                                                                                                <  1:84300/299‑1 (MQ=9)                     aGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTgg                                                                                                                                                                                                                <  1:452334/299‑1 (MQ=9)                           ccTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTa                                                                                                                                                                                                            >  1:563148/1‑297 (MQ=9)                           ccTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTa                                                                                                                                                                                                            <  2:563148/297‑1 (MQ=9)                             ggggtggccAGGCCGTAAACGAGGTCGGATTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGATAAAGCGTTAAGGCGACCGGCCGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAAGTGACGGGTGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGTGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTc                                                                                                                                                                                                        <  2:492773/292‑1 (MQ=11)                                    gcaggacGTAAGTGATGGCGACGGGGAGGGTGTGCCGTTGGGGGGTGGCTTCCGGGGCTAACGCGTTAAGTCGACCGCCGGGGGAGGGCGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTGAATTCGATGCAACGCGAAGAACCTTACCTGGTCTGGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAAc                                                                                                                                                                                                 <  2:613563/293‑1 (MQ=11)                                      aCGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGa                                                                                                                                                                                               >  1:624783/1‑299 (MQ=9)                                               aCGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAAccc                                                                                                                                                                                      <  1:48423/299‑1 (MQ=9)                                                cGAGGTCGACTGGGGGGTTTGGCACTTGAGGCGTTGAGTCCGGGGCTAGCGGGTTAAGTCGACCGCCTGGGGAGTCGGCCCGCAAGGTTAAATCTCAAATGAATTGGCGTGGGCCCTCACAAGCGGTGGAGCATGTGGTTTAATTCGGTGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCt                                                                                                                                                                                     <  2:151505/299‑1 (MQ=11)                                                        aCTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCttt                                                                                                                                                                             <  1:645267/299‑1 (MQ=9)                                                         cTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTttg                                                                                                                                                                            <  1:632912/299‑1 (MQ=9)                                                          ttGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTttgt                                                                                                                                                                           >  1:201387/1‑299 (MQ=9)                                                                  tgtgCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCg                                                                                                                                                                   <  1:112534/299‑1 (MQ=9)                                                                   gtgCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCgg                                                                                                                                                                  <  1:198428/299‑1 (MQ=9)                                                                     gcctggggggCGGGGCTTCCGGAGTTAACGCGTTAAGTCGACCGCCGGGGGAGGACGGCCGCAAGGTTAAAACTAAAATGAATTGACGGTGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTc                                                                                                                                                                <  2:624783/291‑1 (MQ=11)                                                                     gCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTc                                                                                                                                                                >  1:590707/1‑299 (MQ=9)                                                                      ccgttGAGGGGTGGCTTCCGGAGCTATGGCGTTAAGGCAAACGCCTGGGGAGTACGGCCGGAAGGGTAAAACTCAAATTAATTGACGGGTGCCCGCACAAGCGGGGGAGCATGTGGTTTATTTCGATGCAACGCGAAGAACCTTACCTTGTCTTGACATCCAAAGAGATGTCCAGAGATG‑GATTGGTTACTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGGTGGGTTAAGTCCCGGAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTcc                                                                                                                                                               <  2:409571/296‑1 (MQ=11)                                                                         ttGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTccggc                                                                                                                                                            <  1:510523/299‑1 (MQ=9)                                                                          tGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCCTCGGGAACTGTGAGACAGGCGCTGCATGGCTGTCGTCCGCTCGTGTTGTGAAATTTTGGGTTAAGTCCCCACACGAGCGCCACCCCTATCCTTTGTTTCCcgcgtgcccccc                                                                                                                                                           >  2:614862/1‑291 (MQ=11)                                                                                           aGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGAc                                                                                                                                         >  1:250073/1‑299 (MQ=9)                                                                                                      tAAGTCGACCGCCTGGGGAGTACGGCTGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTAGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATa                                                                                                                              <  1:627676/299‑1 (MQ=255)                                                                                                      tAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAAATGACGGGGGCCCGCACAAGCGGGGGAGCATGGGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACGGccgtggata                                                                                                                              >  1:494774/1‑292 (MQ=11)                                                                                                              gggCCGGGGGGGTCCGGCGGTCAGGTTAAACCTAAAAGGATTGGGCGTGGCCCCGCACAAGGGGGGGGGCATGTGGTTTAATTTGGTGCAACGAGAAGAACCTTACCGGGTCTTGAGATCCACAGAACTTTGCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTggag                                                                                                                      <  2:201387/297‑1 (MQ=11)                                                                                                                gCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTggagga                                                                                                                    <  1:91630/299‑1 (MQ=9)                                                                                                                   tGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCg                                                                                                                                                                                                                                <  1:8944/189‑1 (MQ=11)                                                                                                                   tGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCg                                                                                                                                                                                                                                >  2:8944/1‑189 (MQ=11)                                                                                                                     gggAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGGGATAAACTGGAGGAAGGTg                                                                                                               >  1:467374/1‑299 (MQ=9)                                                                                                                        aGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTgggg                                                                                                            >  1:647418/1‑299 (MQ=9)                                                                                                                         gtccGGCAGAAAGGTTAAAGCACAAAGGAATGGGCGGGATCACGCACAGGAGGTGGAGAGTGTGGGTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCGTAGATG‑GATTGGTGCTGTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGa                                                                                                           <  2:293573/296‑1 (MQ=11)                                                                                                                          tACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGAt                                                                                                          <  1:245658/299‑1 (MQ=9)                                                                                                                             gggcggaaGGTTAAACGTCGAAGGAATTGGCGGGGGCCCGCAAGAGCGGGGGGAGATGTGGTTTAATTCGAGGCAACGCGAAAAACCTTACCTGGTCTTGACATCCCCAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACGGTGAGACAGGGGCTGCAGGGCTGTCGTCAGCTCGTGTTGTGAAATGGTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGAc                                                                                                       <  2:611504/293‑1 (MQ=11)                                                                                                                                   aaGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAg                                                                                                 >  1:418847/1‑299 (MQ=9)                                                                                                                                       ttAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcat                                                                                             <  1:16095/299‑1 (MQ=9)                                                                                                                                        tacaaCTCAAAGGAATTGACGGGGGCCCTGCCAAGAGGGTGGGCATGTGGTTTAAGTCGATGCAACCCGAAGAACCTTTCCTGGTGTGGAGATCCAGCGAACTTTCCAGAGATG‑GATTGGTGCCTGCGGGAAATGTGAGACAGGTGCTGCATGGCTTTCGGCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGtcatc                                                                                            <  2:612803/296‑1 (MQ=11)                                                                                                                                          aaaCTCAAATGAATTGACGGGGGCCCGCACAAGCGGTTGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCCGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGAATGGCTTTCGTCACCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGACCGCAACACTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAAGT‑CAAAGGCGTATGCCAGTGAAAAACTGGAGGACGGGGGGGATGAcagcacagatacat                                                                                          >  2:628101/1‑286 (MQ=11)                                                                                                                                           aaCTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCCACGCGAAGCACCTTACCTGGTCTTCACATCCACAGAACTTTCCAGAGATG‑GATTGGTACCTTCGGGAACTGTGAGACAGGTGCTGCATTGCTGTCGTCAGCTCGTGTTGTGAACTGTTGGGTTAAGTCCCGCAACGAACGCCACCCTTACCCTTTTTTGCCCGGCGTCCCGCCGGGAAAC‑TAAAGGGACCTGCCATGGACAAACTGGCGGAAGGTGGGTATACCGCCAAGtcatactc                                                                                         >  2:314948/1‑295 (MQ=11)                                                                                                                                            aCTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTttgttg                                                                                                                                                                         >  1:370599/1‑219 (MQ=255)                                                                                                                                            aCTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTttgttg                                                                                                                                                                         <  2:370599/219‑1 (MQ=255)                                                                                                                                             cTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGc                                                                                       <  1:196358/299‑1 (MQ=9)                                                                                                                                              tCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGgaga                                                                                                                                          <  2:26738/247‑1 (MQ=9)                                                                                                                                              tCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGgaga                                                                                                                                          >  1:26738/1‑247 (MQ=9)                                                                                                                                               cAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCCGCGCGTGTTGTGAAATGTTGGTGTAAGTCCCGCAACGAGCACAACCCTTATCCTTTGTTGCCTCGCGTCCGGCCGGGAACTACACGGCAGAC‑GCCAGTGTTAACCCTGAAGTAGGTGGGACAGCCGTCAAGTCCTCATGGctc                                                                                     >  2:86983/1‑297 (MQ=11)                                                                                                                                                 aaTGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCACAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCtt                                                                                   >  1:287919/1‑299 (MQ=9)                                                                                                                                                            gggggCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGc                                                                        >  1:608094/1‑299 (MQ=9)                                                                                                                                                                      acaAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTg                                                              >  1:382312/1‑299 (MQ=9)                                                                                                                                                                        aaGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCt                                                            >  1:384041/1‑299 (MQ=9)                                                                                                                                                                          ggggTGGGGGATGTGGGTTAAGTCGATGCAACGCGAAGAACCTTACATGGTCTTGACATCCCCAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGGGACAGGTGCTGCATGGCTGTGGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTag                                                          <  2:79546/297‑2 (MQ=11)                                                                                                                                                                            ggTGGAGCATGTGGTTTAATTCGATACAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACaa                                                        >  1:376098/1‑299 (MQ=9)                                                                                                                                                                                     tgtgGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATa                                               >  1:348823/1‑299 (MQ=9)                                                                                                                                                                                             aaTTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAgaga                                       >  1:92120/1‑299 (MQ=9)                                                                                                                                                                                              aTTGGAGGCAACGGGAAGAACCTTACGTGGTCTTGGCATCCACCGAACTTTCCCG‑GAAGGGATTGGTGCCTTCGGGAACTGTGAGACAGTGGCTGCATGGCTGGCGGCAGGTCGTGTTGTGAAATGTGGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGaa                                      <  2:41180/299‑1 (MQ=11)                                                                                                                                                                                                 cGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCCACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACC‑CACAGGAGACTGCCAGTGCTAAACTGGAGGAAGGTGGGTGTGACGTCAAGTCATAATGGCCCTGACCAACAGGGCAACACACCTGCTCCACAGGCACATTCaacgcgacacg                                   >  2:164002/1‑290 (MQ=11)                                                                                                                                                                                                  gATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGa                                  >  1:390620/1‑299 (MQ=9)                                                                                                                                                                                                      cAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTc                              <  1:29595/299‑1 (MQ=9)                                                                                                                                                                                                       aaCGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTCCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGCACT‑CAAAGGAGACTGCCAGTGCTAAACTGGAGGAAGGAGGGGAAGACGTCAAGTCATCATGGCCCTGCAGACCAGGGGCACACACGGGCTACACTGGAGCCAACAAAGACAAACGACCAcg                             >  2:428506/1‑299 (MQ=11)                                                                                                                                                                                                             aagaaCCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCGTGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTgg                                                                                                              <  1:305252/212‑1 (MQ=11)                                                                                                                                                                                                             aagaaCCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCGTGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTgg                                                                                                              >  2:305252/1‑212 (MQ=11)                                                                                                                                                                                                             aagaaCCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGagag                       >  1:437180/1‑299 (MQ=9)                                                                                                                                                                                                                aaCCTTACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCaa                    <  1:222287/299‑1 (MQ=9)                                                                                                                                                                                                                     tACCTGGTCTTGACATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGa               >  1:49090/1‑299 (MQ=9)                                                                                                                                                                                                                                  cATCCACAGAACTTTCCAGAGATG‑GATTGGTGCCTTCGGGAACTGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGACCTCATAAAGTGc  >  1:504330/1‑299 (MQ=9)                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                 CGTGGGGAGCAAACAGGATTAGATACCCTGGTAGTCCACGCCGTAAACGATGTCGACTTGGAGGTTGTGCCCTTGAGGCGTGGCTTCCGGAGCTAACGCGTTAAGTCGACCGCCTGGGGAGTACGGCCGCAAGGTTAAAACTCAAATGAATTGACGGGGGCCCGCACAAGCGGTGGAGCATGTGGTTTAATTCGATGCAACGCGAAGAACCTTACCTGGTCTTGACATCCACAGAACTTTTCAGAGATGAGAAT‑GTGCCTTCGGGAACCGTGAGACAGGTGCTGCATGGCTGTCGTCAGCTCGTGTTGTGAAATGTTGGGTTAAGTCCCGCAACGAGCGCAACCCTTATCCTTTGTTGCCAGCGGTCCGGCCGGGAACT‑CAAAGGAGACTGCCAGTGATAAACTGGAGGAAGGTGGGGATGACGTCAAGTCATCATGGCCCTTACGACCAGGGCTACACACGTGCTACAATGGCGCATACAAAGAGAAGCGACCTCGCGAGAGCAAGCGGACCTCATAAAGTGC  >  CP006584/4178739‑4179262

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: ‑