Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,193,996	T→G	D70A (GAT→GCT)	thiE ←	thiamine phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,193,996	0	T	G	99.0%	334.8 / NA	99	D70A (GAT→GCT)	thiE	thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (50/48);  total (51/48)
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.87e-01

GCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTC  >  NC_000913/4193751‑4194241                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      gCCGTGCCAGCTGTTCCAGCCCCTGCGTTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGGCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCGGATGCTTGATCGCCAGCCGCCAGTAAGCGTTg                                                                                                                                                                                                                                                  <  1:440021/251‑1 (MQ=255) gCCGTGCCAGCTGTTCCAGCCCCGGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CTTGATGCTTGATCGCCAGCCGCCAGTAAGCGtt                                                                                                                                                                                                                                                   <  1:398816/250‑1 (MQ=255)        cAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAAc                                                                                                                                                                                                                                            >  2:166838/1‑250 (MQ=255)               cccAGCCCCTGCGGGGCAGAAGGCATCTGTGTGGGTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAAgcgcgc                                                                                                                                                                                                                                    <  1:111837/250‑1 (MQ=255)                  aGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATTTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGtt                                                                                                                                                                                                                                 <  1:77602/251‑1 (MQ=255)                   gCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTta                                                                                                                                                                                                                                <  1:402123/251‑1 (MQ=255)                         gCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgc                                                                                                                                                                                                                            <  1:207876/249‑1 (MQ=255)                         gCGGTGCAGAAGGCATCTGTTGGGTTGGCGGCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGGGCTGCCAGCGCGACGTCGAGTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCATGCCTGCCGCGCGGATGGTATTGAGATCGGTGGCTTGCAAAGCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAgcgg                                                                                                                                                                                                                          <  1:229439/251‑1 (MQ=255)                          cGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAAAAATCGCGCGTTATAgcggc                                                                                                                                                                                                                         >  2:171643/1‑251 (MQ=255)                                aGAAGGCATCTGTTTGGTTTGCGTCGGCAACACTTGTCCCAGCGCGATATAAGTGGGGCGGGCTTCCAGCGCGACGCAGATTTCCATA‑CCGTCATGTGTC‑TAAACGCCCAGCCGAAGGCCGGCCGCTCGGATGGCATTGAGATCGGTGGCGTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGAGCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCa                                                                                                                                                                                                                   <  2:521715/251‑1 (MQ=255)                                aGAAGCCATCTGTATCATTTGCGCAGGGACCAAATGTCCTAGCGCGAAACAAGAGTGCCGTGCTGCCCGCGCGACGTAGTTTTCCATA‑TCGACATGCCTA‑GAACCGCCCAGCCGCCCGCCCGCCGAGCGGATGGCACTGAGATCGGTGACTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGCAATAGTTGATAAACAATCGCGCGTTATAGCGGCGGcccc                                                                                                                                                                                                                   >  2:36736/1‑250 (MQ=255)                                       aTCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAt                                                                                                                                                                                                            >  2:95423/1‑251 (MQ=255)                                           gtttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAt                                                                                                                                                                                                            <  1:276497/247‑1 (MQ=255)                                            tttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTTCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGACCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCCATTgccc                                                                                                                                                                                                       >  1:268467/1‑250 (MQ=255)                                            tttggtttgCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGCTTTACATA‑TCGTCATGTGCC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTgccg                                                                                                                                                                                                       >  1:356528/1‑251 (MQ=255)                                                        cGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCCGCGCAATTGCCGCCACGACATCgg                                                                                                                                                                                           >  2:131423/1‑251 (MQ=255)                                                         gggAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGACGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGCCGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCCATTGCCGCCACGACATCGGc                                                                                                                                                                                          >  1:250066/1‑251 (MQ=255)                                                          ggaccacaTGTCCCAGCGCGATATAAGAGGGGCTTGCTGCCAGCTCGACGTCGATTTACATA‑CCGTCATGTGTC‑TAAACGCCTAGCCTCAGGCCTGCCGGGCGGAGGGCTTTGAGATCGGGGGCTTGCAAATCTCCCTGCCCCAAATGGACGCCATCCG‑CTGGCTGCTTGTTCGCCACCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATGGCCGCCACGACATCGGCt                                                                                                                                                                                         <  1:290115/247‑1 (MQ=255)                                                             acacATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTcc                                                                                                                                                                                      >  1:256681/1‑251 (MQ=255)                                                                   gTCCCCGCGCGATATAAGAGGGGCGTGCTGACATCGCGACGTCGCTTTCCATA‑TCGACATGTGTC‑GAAACGACCAGCCGCAAGCCTGCCGCGCGGATGGCATTGCGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAACGCCATTGCCGCCACCACATCGGCTTACACctct                                                                                                                                                                                >  2:530817/1‑251 (MQ=255)                                                                    tCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCtt                                                                                                                                                                               >  1:194181/1‑251 (MQ=255)                                                                     cccAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTc                                                                                                                                                                              <  1:276372/251‑1 (MQ=255)                                                                       cAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCAGGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCAt                                                                                                                                                                            <  1:293433/251‑1 (MQ=255)                                                                        aGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATc                                                                                                                                                                           >  2:101781/1‑251 (MQ=255)                                                                         gcgcTATAAAAGGGGGGCGGGCGGCCAGCGCTCCGTCTATTTCCATA‑TCGTCA‑GGGTCCGAAACGCCCAGCCGCGGGCCTGCCGCGCGGATGGATTTGAGACCGGTGGCTTGCAAATCTTCCTGCCCCAAATGAACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAATGCCGCCACGACATCGGCTTCCTCCTCTTCATcg                                                                                                                                                                          <  2:372245/251‑1 (MQ=255)                                                                            cgATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCgcgc                                                                                                                                                                       >  1:479069/1‑251 (MQ=255)                                                                               tataAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCg                                                                                                                                                                     >  1:484310/1‑250 (MQ=255)                                                                                 aaaGAGGGGCTTCCTGCCAGCGCGCCTGCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGAGGGCATTGAGATCGGTGGCTTGAAAACCTTCCTGCCACAAATGGCCGCCAAACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCTGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCAGCGCGTCGATc                                                                                                                                                                  <  2:523079/250‑1 (MQ=255)                                                                                    gagtggCGTGCTGCCAGCTCGAAGTCTATT‑ACATACGCGTTATGTGTC‑GAAACGCCCAGCCTGAGGCCTGCCGCTGGGATGGTATTGATATCGGGGGCTTGCAAATGTTCCGGCCCCAAATGGAGGCCGTACG‑CCTGATGCTTGATCGCTAGGCTCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCTGCCACGACATCGGCTTCCACCTCGTCAGCGCGCCGATCttt                                                                                                                                                               <  2:526712/247‑1 (MQ=255)                                                                                     ggtggCGTGCTGCCAGCGCTCCGTCGATTTCCATA‑TCGTCATGGGTC‑GAACCGCCCAGCCGCAGGCCTGCCGCGGGGATGGCTTTGAGATGGGTGGCGTGCAAATCTCCCTGCCCCAAATGGCCGCCATACG‑CCTGATGTTTGATCGCCAGCCGCCAGTAAGCTTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCGCCTCTGCATCGCGCCGATCTTTg                                                                                                                                                              <  2:43968/248‑1 (MQ=255)                                                                                       gggCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GCAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATAGGTAGCTTGCAAATCTTCCTGCCCCAACTGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCCGTAAGCGTTGATAAACAATc                                                                                                                                                                                                                                        >  1:512186/1‑175 (MQ=255)                                                                                          cGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCCCGACATCGGCTTCCACCTCTTCCTCGCGCCGATCTTTGCTGCg                                                                                                                                                         >  2:441267/1‑251 (MQ=255)                                                                                            tgctgcCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGCTCTTTGATGCGt                                                                                                                                                        >  1:523536/1‑250 (MQ=255)                                                                                                   aGCGCGACGTCGATTTCCATA‑TAGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTgg                                                                                                                                                  >  1:477910/1‑249 (MQ=255)                                                                                                     cgcgACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGgag                                                                                                                                                >  1:65824/1‑249 (MQ=255)                                                                                                        gACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACCTCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgtac                                                                                                                                           >  1:436811/1‑251 (MQ=255)                                                                                                          ccgtcgattcccagtgcggatggtgtC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATCCG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCGTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgtacg                                                                                                                                          <  2:285760/229‑1 (MQ=255)                                                                                                            tCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAgtacgtac                                                                                                                                       >  2:186458/1‑251 (MQ=255)                                                                                                             cGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGtacgtacg                                                                                                                                      >  1:414760/1‑251 (MQ=255)                                                                                                             cGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGtacgtacg                                                                                                                                      >  1:468393/1‑251 (MQ=255)                                                                                                               aTTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGcc                                                                                                                                    >  1:7376/1‑251 (MQ=255)                                                                                                                 ttCCATA‑TCTTCATGTGTG‑AAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCt                                                                                                                                   <  1:226559/250‑1 (MQ=255)                                                                                                                   cctga‑tCGTCATGTGTG‑GAATTGCCCAGCGGAAGCCCTGCCGCGCGGATGGCATGGACATGGGGGGGTTTACAACCTCCCCGGCCCAAATGGACGCTGTACG‑CCTGATGTTTAAGCGCCAGCCGCCAGTATGCGGTGATAAACAACCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCGCGACATCGGCTTCCTCCTCTTCATCGCGCCGATCTGTGAGGCGTAGTGGGAGAGGACGTACGCCTGCa                                                                                                                                <  2:321137/247‑1 (MQ=255)                                                                                                                           gTCAGGTGGC‑GAAACGCCCAGCCGCAGGCCGGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAACCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGAGCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACa                                                                                                                         <  2:152593/251‑1 (MQ=255)                                                                                                                             cATGTGTG‑GATACGCACAGCAGCAGGCATGCAGCGAGGACGGCATTGAGATCGGTGGCCTGCAAATCTTCCTGCCCCAAATGGCCGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACACTCGCGCGTTATAGCGGCGGCCAAGCACAATTGCCGCCTCGGCATCAGCTTCCACCCCTTCATCGCGCCGCACTTTTGTGCGCAGCTGGAGAGTACCTACGCCTGCCTCCACCAgc                                                                                                                       >  2:368210/1‑250 (MQ=255)                                                                                                                              aTGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGGACGTACGCCTGCATCCAACAGAc                                                                                                                      >  1:480771/1‑251 (MQ=255)                                                                                                                                gtgtC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGt                                                                                                                    >  1:283399/1‑251 (MQ=255)                                                                                                                                      gAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCCTCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCCTCCAACAGACGTTCGAt                                                                                                               >  2:25559/1‑251 (MQ=255)                                                                                                                                        aaCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATAGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTAACTACGACTGCATCCCACAAACGTTCGATcc                                                                                                             >  2:379285/1‑251 (MQ=255)                                                                                                                                        aaCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑ACTGATGCTTGATCGCCAGCCGCCAGTACGCGTTGATAAACAATCCCGCGTTATAGCGGCGGCCCCGCGCAATTGCCGCCACGACATCGGCTTCAACCTCTTCATCGCGCCGATCTTTTATTCCTCGCTGGAGAGTACGTACGCCTGCATCCAACatgagtttgatcc                                                                                                             >  2:283922/1‑239 (MQ=255)                                                                                                                                         aCGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCAc                                                                                                                                                                                                   >  1:240531/1‑164 (MQ=255)                                                                                                                                         aCGCCCAGCCGCAGGCCTGCCGCGCGGAGGGCATGGAGATCGGTGGCTTGCAAATCTTCCGGCCCCAATTGGACGCCATACG‑CCTGATGCTGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCAc                                                                                                                                                                                                   <  2:240531/164‑1 (MQ=255)                                                                                                                                           gCCCAGCCGCAGGCCTGCCGCGCGGATGGCCTTGCGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACt                                                                                                          >  1:440410/1‑251 (MQ=255)                                                                                                                                               aGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACTCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGt                                                                                                        <  1:198296/249‑1 (MQ=255)                                                                                                                                               aGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCCTCCAACAGACGTTCGATCCACTGt                                                                                                        >  2:198296/1‑249 (MQ=255)                                                                                                                                                  cGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCACCAGACGTTCGATCCACTGTACGCt                                                                                                   >  2:232328/1‑251 (MQ=255)                                                                                                                                                      ggCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGt                                                                                                                                                                                                                                  >  1:494110/1‑120 (MQ=255)                                                                                                                                                       gccgcccgtgcggtgggtgtgaacaTCGGGGGCTTGAAAATGTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGTTTTACCGCCGGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTcca                                                                                              <  2:160459/227‑1 (MQ=255)                                                                                                                                                        ccTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTccac                                                                                             <  1:163987/251‑1 (MQ=255)                                                                                                                                                              gcgcgGATGGCATGGAGATCGGGGGCTTGAAAACCTCCCTGCCCCAAATGGACGCCATACG‑CCTGAGGCTTGATGGCCAGCCGCCAGTAAGCGTGGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGATAGTACGGACGCCTGCAGCCAACAGACGTTCGAGCCACTGTACGCTGTCCACCACCgg                                                                                       <  2:99075/251‑1 (MQ=255)                                                                                                                                                                 gggATGGCATTGAGATCGGTGGCTTGAAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGAGCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTa                                                                                    <  1:350038/250‑1 (MQ=255)                                                                                                                                                                 cgGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACt                                                                                                          >  1:188340/1‑229 (MQ=255)                                                                                                                                                                        attggagaTCGGGGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATCCG‑CCTTATGCTTTATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCt                                                                             <  1:248711/247‑1 (MQ=255)                                                                                                                                                                         aTTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTg                                                                            >  1:376790/1‑251 (MQ=255)                                                                                                                                                                            gagaTCGGTGGCTTGCAACTCTTCCCGCCCCAAAGGGGCGCCATACG‑CCGGATGCTGGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGGCCTGAAc                                                                         <  2:233344/251‑1 (MQ=255)                                                                                                                                                                              gaccggtggcgtgaaaacccccccgCCCCAAAGGGACCCCATCCG‑CCGTATGCTGGAACCCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTGATAGCGGCGGCCCAGCGCAATTGCCGCAACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGa                                                                       <  2:13920/227‑1 (MQ=255)                                                                                                                                                                               accgtgggtttGCAAATCTTCCTGCCCCAAATGGACGCCATCCG‑CCTGTTGCTTGATCGCAAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATGGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCGGAACGaa                                                                      <  1:390850/242‑1 (MQ=255)                                                                                                                                                                                 gggTGGCTTCAAAATCTCCCCGCCCAAAAGGTACGCCATACG‑CTTGTTGCTTGATGGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATAGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCAATGTACGCTGTCCACCACCGGGTACAGTCCTGAACGaaaa                                                                    <  2:151389/250‑1 (MQ=255)                                                                                                                                                                                         ggCAAATCGTCCCGCCCCAAATGGACGCCATCCG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACagg                                                            <  2:472872/250‑1 (MQ=255)                                                                                                                                                                                          gCAAATCTTCCGGCCCCAAAGGGACGCCATACG‑CCGGATGCTTAGCCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACagga                                                           <  2:65824/251‑1 (MQ=255)                                                                                                                                                                                           cAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGCGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGAACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGAACagg                                                            >  1:462284/1‑249 (MQ=255)                                                                                                                                                                                           cAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACaggag                                                          <  1:166838/251‑1 (MQ=255)                                                                                                                                                                                           cAAATCTTCCGCCCCCAAATGGACGCCATACG‑CCTGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTGATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTGCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACaggag                                                          <  2:322030/251‑1 (MQ=255)                                                                                                                                                                                                  cccTGCCCCAAATGGACGCCATACG‑CCGGAGGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATc                                                   <  1:441267/250‑1 (MQ=255)                                                                                                                                                                                                  cccTGCCCCAAATGGACACCA‑ACGGCCGGATGGTTGACCGCCGGCCGCCGGTAAGCTTTGATAAACAATCGCTCGTTATAGCGGCGGCCCAGCGCAATTTCCGCCACGACATCGGCTTCCACCTCTTCAGCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACTTACGCCTGCATCCAACAGACGTTCGATCCACTGAACGCTGTCCACCACCGGTTACAGTCCTGAACGAAAAGGTACAGGAGGaaaac                                                    <  2:33187/249‑2 (MQ=255)                                                                                                                                                                                                      gCCCCAAATGGACGCCATACG‑CCTGATGCTTGACCGCCAGCCGCCAGTAAGCGTTGATAACCAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAg                                                 <  1:141512/249‑1 (MQ=255)                                                                                                                                                                                                         ccAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTg                                             >  1:95460/1‑250 (MQ=255)                                                                                                                                                                                                            aaTGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATAc                                         >  1:36656/1‑251 (MQ=255)                                                                                                                                                                                                                gACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAACGGTACAGGCGGAAAATCAGGCTGATACATcc                                     >  1:531955/1‑250 (MQ=255)                                                                                                                                                                                                                 ccGCCAGACG‑CCTGATGCTGGTCCGCCGGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAAGTGCCGCTCCGACATCGGCTTCCCCCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCAc                                    <  2:381359/250‑1 (MQ=255)                                                                                                                                                                                                                    ccATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGCAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCt                                 >  2:250392/1‑251 (MQ=255)                                                                                                                                                                                                                         ccgccggaggctggccggCCAGCCGCCAGTAAGCGTTGAGAAACAAACGCGCGTTATAGCGGCGGCCCGTCGCAATGGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGACGCGTAGCTGGAGAGTACGTACGCCGGCATCCAACATACGTTCGATGCACTGTACGCTGTCCACCGCCGGGTACAGTCCTGAACGAAAAGGTACAGTAGGAAAATCAGGCTGCTACATCACGCTtcc                              <  2:108903/233‑1 (MQ=255)                                                                                                                                                                                                                         cG‑CCTGATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTtcctc                            >  1:132363/1‑251 (MQ=255)                                                                                                                                                                                                                               gattctggATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCCCCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGGATCCAACAGACGTGCGATCCACTGTACGCTGTCCACCTCCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTAc                       <  1:256552/244‑1 (MQ=255)                                                                                                                                                                                                                               gATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTa                        <  1:171643/250‑1 (MQ=255)                                                                                                                                                                                                                               gATGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTAc                       >  1:508523/1‑251 (MQ=255)                                                                                                                                                                                                                               gATGCTTGATCCCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGTTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACGTCACGCTTCCTCCTTAc                       <  2:356528/251‑1 (MQ=255)                                                                                                                                                                                                                                aTGCTTGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACg                      <  1:25559/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggtttgacccccAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGc                     <  2:169238/241‑1 (MQ=255)                                                                                                                                                                                                                                    ttGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAgg                  >  2:307730/1‑251 (MQ=255)                                                                                                                                                                                                                                     tGATCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCaccacc                                                                                         >  1:339119/1‑179 (MQ=255)                                                                                                                                                                                                                                     ggaccGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCGGTCCaccacc                                                                                         <  2:339119/175‑1 (MQ=255)                                                                                                                                                                                                                                        tCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCa                                                  >  1:430174/1‑215 (MQ=255)                                                                                                                                                                                                                                        tCGCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCa                                                  <  2:430174/215‑1 (MQ=255)                                                                                                                                                                                                                                          gCCAGCCGCCAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGAtt            <  1:273266/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ccAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCCATTGCCGCCACGACATCGGCTTCCACATCTTCATCCCGCCGATCTTTGCTGAGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGcc    >  1:280236/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGGACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCct   >  1:267435/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cAGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTAATCGCGCCGATCTTTGCTGCGTAGCTGGAGAGTACGTACGCCTGCATACAACAGCCGTTCGATCCACTGTACGCTGTCCACCACCGGTGACCGTCCTGACCGAAAAGGTACAGGCGGCAAATCAGGCTGAAACCTCACCCTTCCTCCTTACGCCGGTCGATTTCCCCGccc   >  2:180088/1‑250 (MQ=255)                                                                                                                                                                                                                                                    aGTAAGCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCctc  >  2:166178/1‑251 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      GCCGTGCCAGCTGTTCCAGCCCCTGCGGTGCAGAAGGCATCTGTTTGGTTTGCGTCGGGAACACATGTCCCAGCGCGATATAAGAGGGGCGTGCTGCCAGCGCGACGTCGATTTCCATA‑TCGTCATGTGTC‑GAAACGCCCAGCCGCAGGCCTGCCGCGCGGATGGCATTGAGATCGGTGGCTTGCAAATCTTCCTGCCCCAAATGGACGCCATACG‑CCTGATGCTTGATCGCCAGCCGCCAGTAATCGTTGATAAACAATCGCGCGTTATAGCGGCGGCCCAGCGCAATTGCCGCCACGACATCGGCTTCCACCTCTTCATCGCGCCGATCTTTGATGCGTAGCTGGAGAGTACGTACGCCTGCATCCAACAGACGTTCGATCCACTGTACGCTGTCCACCACCGGGTACAGTCCTGAACGAAAAGGTACAGGAGGAAAATCAGGCTGATACATCACGCTTCCTCCTTACGCAGGTAGATTTCTCCGCCTC  >  NC_000913/4193751‑4194241

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