Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,158,814	0	G	A	77.5%	125.1 / 26.2	76	T143T (ACG→ACA)	mdtC	multidrug efflux pump RND permease subunit MdtC
Reads supporting (aligned to +/- strand):  ref base G (5/11);  new base A (29/26);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGCGCGCTGGAAGATGGCACTC  >  NC_000913/2158574‑2159058                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     tGCCACGCCGCTGGAGCGTGCTCTTTGGCGCATGGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATCTTGCAGTTTGATTTTGGCCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGAGGCCCCGCCGCCCGACCT‑ATCGCAAAGCTAGCGCGTCGGAGGCGCCAATTATGATCCTTACGCttccgtcc                                                                                                                                                                                                                                             <  2:586029/250‑8 (MQ=255)  gCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCg                                                                                                                                                                                                                                            <  1:580872/250‑1 (MQ=255)    cACGCCGCTGGAGTGCTCACTTGTGCGCATTGCCGG‑AGGCAGTGACATGACCTCCAGCAGGTCGCTCGGCAGTACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGTGGCGAACCACGCTGCACAAACTTTGCTGCCCAGGGGGATGTCCAGATGCCCGACCT‑ATCGCAAAGCGAAACCGTCGGCTGCGCCAATTATGATCCTCACGCTGACggtcgatt                                                                                                                                                                                                                                         <  2:464765/251‑8 (MQ=255)    cACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCATCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCGCAGTGGGATGCCCAGCCGCCCGACCT‑AGCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACGGTCGTCCGATa                                                                                                                                                                                                                                         <  2:270944/251‑1 (MQ=255)             ggAGCGCTCGCTAGGGCGAATTGCCGGCAG‑CAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGAACGCGTATTATTTTGCAGATTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGATGCACAACGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTAGGATGCTCCCATTATGATCCTCACCCTGCCGTCCGATACTTCTTCGc                                                                                                                                                                                                                                >  2:554222/1‑251 (MQ=255)               aGCGCTCACTTGGGCGCATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGCCCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACAATGACGTCCGATACTTATTCGc                                                                                                                                                                                                                                >  2:273160/1‑249 (MQ=255)                    tCACTTGGGCGCATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTg                                                                                                                                                                                                                          >  1:274607/1‑250 (MQ=255)                       cTTGGGCGCATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAAACCGTCGGATGCGCCAATTATGATCCTCACGCTGCCGTCCGATACTTATTCGCAGGGTGCACt                                                                                                                                                                                                                      >  2:212194/1‑251 (MQ=255)                            gcgcATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACg                                                                                                                                                                                                                 <  1:93930/251‑1 (MQ=255)                              gcATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGAt                                                                                                                                                                                                               <  1:518185/251‑1 (MQ=255)                                  tGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCg                                                                                                                                                                                                           >  1:471884/1‑251 (MQ=255)                                     cGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATATTGCAGTATGATTTTGACCGGGATAGCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑AGCGCAAAGCGAATCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCt                                                                                                                                                                                                        <  1:615460/251‑1 (MQ=255)                                      gg‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTc                                                                                                                                                                                                       >  1:78668/1‑251 (MQ=255)                                            cAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACg                                                                                                                                                                                                   <  1:12346/250‑1 (MQ=255)                                                gAAAGGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGGCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCt                                                                                                                                                                                              <  2:330747/251‑1 (MQ=255)                                                   aTGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGc                                                                                                                                                                                           >  2:339068/1‑251 (MQ=255)                                                     gACCTCCAGCATTTCGCTCGTCAGTACGGGTATTATTTTGCAGTTTTATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTTCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCGCACACTGACGTCCGATACTTATTCGCAGTGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTc                                                                                                                                                                                         <  2:312714/251‑1 (MQ=255)                                                            agcagTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGGCCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGAt                                                                                                                                                                                  <  1:504386/251‑1 (MQ=255)                                                             gcagTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACCCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGAtt                                                                                                                                                                                 >  2:500422/1‑251 (MQ=255)                                                               agTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTc                                                                                                                                                                               >  1:512715/1‑251 (MQ=255)                                                                   cGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCa                                                                                                                                                                            >  1:378169/1‑250 (MQ=255)                                                                    gCTCGGCAGCACGCGTATTATTTTGCAGTTTGATGTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCaaa                                                                                                                                                                          <  2:43522/251‑1 (MQ=255)                                                                        ggCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGCTGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCCGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCg                                                                                                                                                                       >  1:164793/1‑250 (MQ=255)                                                                         gcagcaCGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGAc                                                                                                                                                                     <  1:273160/251‑1 (MQ=255)                                                                          cagcaCGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACg                                                                                                                                                                    <  2:143864/251‑1 (MQ=255)                                                                           agcaCGCGTATTATTTTGCAGTTTGGTTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCGGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACGGTACGATTTCGCCTCGACGCAGCTGGTTCCGACGATTTCACAAAGCGACgg                                                                                                                                                                   <  2:385448/251‑1 (MQ=255)                                                                                       attTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCCAATCGACGGTGTTGGTgatgt                                                                                                                                                       >  1:275155/1‑251 (MQ=255)                                                                                       attTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTgatg                                                                                                                                                        <  2:165495/250‑1 (MQ=255)                                                                                           tGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTgatgtcga                                                                                                                                                    <  1:100953/250‑1 (MQ=255)                                                                                              aGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTgatgtcgatgtc                                                                                                                                                >  1:72693/1‑251 (MQ=255)                                                                                                     ttttGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGc                                                                                                                                          >  2:169359/1‑250 (MQ=255)                                                                                                         gACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTc                                                                                                                                     <  1:174247/251‑1 (MQ=255)                                                                                                          aCCGGGATATCAACGGCGCAGCGCGTGATGTGGAGGCGGCGATCAACGCTGCACAAAGTTTGCTGACCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATGTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAAGCGACGGTGTTGGTGATTTCGATGTCGGAGGCAGCTCa                                                                                                                                    <  2:142864/251‑1 (MQ=255)                                                                                                             gggATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTg                                                                                                                                 >  1:268940/1‑251 (MQ=255)                                                                                                                    cAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCAAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGGGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCccgccg                                                                                                                          >  2:506628/1‑251 (MQ=255)                                                                                                                     aaCGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGt                                                                                                                         >  2:521960/1‑251 (MQ=255)                                                                                                                     aaCGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGt                                                                                                                         >  2:158065/1‑251 (MQ=255)                                                                                                                      aCGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGt                                                                                                                         >  2:488561/1‑250 (MQ=255)                                                                                                                       cGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCACAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTAc                                                                                                                       >  1:406089/1‑251 (MQ=255)                                                                                                                          cgcggcgcGTGATGTTCACGCGGCGATCAACGCTGCACGAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAGGCGAACCCGTCGGATGCGCCAGTTATGATCCTCACACGGACGTCCGAGACTTAGGCGCAGGGTGAGCGGTACGATTTCGCCTGGACGCAGCTGGCTCCGACGATTGCGCAAATCGACGGTGTTGGGGATGTGGATGTCGGAGGCAGCTAACTGCCCGCCGTAcgcg                                                                                                                    <  2:78178/247‑1 (MQ=255)                                                                                                                          cgcAGCGCGTGATGTGCAGGCGGCGATGAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCCGCCGCCCGACCT‑ATCCCAATGCGACCCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGCTACTTATTCGCAGGGTGAACTGTATGATTTCGCCTCGACGCAGCTGTATCCGACGATTTCGCAACTCGACGGTGTTGGTGATGTTGGTGTCGGAGCCAGCGCACTGcacaccggacccc                                                                                                                    >  2:74483/1‑239 (MQ=255)                                                                                                                          cgcAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCCAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTAcgcg                                                                                                                    >  1:283169/1‑251 (MQ=255)                                                                                                                           gcAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGACAGCTCACTGCCCGCCGTACGCGt                                                                                                                   >  1:157402/1‑251 (MQ=255)                                                                                                                              gcgcgagctgtgCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCgg                                                                                                                >  2:622369/9‑251 (MQ=255)                                                                                                                               ggcttGATGTGCGGCCGTCGATCAACGCGGCACAAAGTTTTCTGCCCAGTGGGTTGCCCAGCCGCCCGCCCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACCCTGACGTCCGATACTTATTCGCAGGGTTAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGTGGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCggg                                                                                                               <  2:376355/247‑1 (MQ=255)                                                                                                                                 gttGATGTGCAGGCGGCGAACTGCGATGCAAAAATTTGGCGTCCCAGTGGGGTGCCCAGCTGCCCG‑CCTGACCGCAAAGCGACCCCTTCGGTTGCACTAATTATGACCCTCCCGCGGTCGTCCGATACTTGTTCGCAGGGTGAACTGTACGTGTTCGCCTCGTCGCTGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGTCAGTTCACTTCCTGCCGTACTCGTCGTGCt                                                                                                             <  2:444517/249‑1 (MQ=255)                                                                                                                                   tGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTg                                                                                                            >  2:317791/1‑250 (MQ=255)                                                                                                                                               gaccttCAACGCTGCACAAAGTTTGCTGCCCAGGGGGATGCCCAGCCGCCCGACCT‑ACCGCAAAGCAAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATGGCGCAAATCGACGGTGTTGTTGAGGTCGATGTCGGAGGCAGCCCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGgcgc                                                                                               <  2:382914/246‑1 (MQ=255)                                                                                                                                                gcgATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCt                                                                                              >  1:272411/1‑251 (MQ=255)                                                                                                                                                  gATCAACGCTGCACAAAGTTTGCTGACCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGGCGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGTATCCGCCGGCGCTGt                                                                                            >  2:457824/1‑251 (MQ=255)                                                                                                                                                   aTCAACGCTGCACAAATTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGCCCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTAGGATCCGCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATTTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGtt                                                                                           <  1:1503/251‑1 (MQ=255)                                                                                                                                                   aTCAACGCTGCACAAAGTTTGCTGCCCAGTGGGAGGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGAGCCGCACACTGACGTCCGATTCTTATCCGCAGTGTGAACTGTGCGATTGCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAATTCGACGTTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGtt                                                                                           <  1:554104/251‑1 (MQ=255)                                                                                                                                                    tCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGgcgc                                                                                               >  1:382914/1‑246 (MQ=255)                                                                                                                                                     cAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGACCCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGttt                                                                                          <  1:559522/250‑1 (MQ=255)                                                                                                                                                     cAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGTTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTa                                                                                         <  1:599494/251‑1 (MQ=255)                                                                                                                                                     cAACGCTGCACAAAGTTTGCGGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACGGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTTTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTa                                                                                         <  2:78668/251‑1 (MQ=255)                                                                                                                                                            gCACAAAGTTTGCTGCCCGGTGGGATGCCCAGCCTCCCGCCCT‑ATCGCAAAGCGAACCCGTGGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGGGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAAGCCGCAGGCGCTGTTTAATCAggg                                                                                  <  2:460416/251‑1 (MQ=255)                                                                                                                                                              aaaaaGTTTGCGGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGACCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCg                                                                                <  1:490347/249‑1 (MQ=255)                                                                                                                                                                       gctgcCAAGTGGGAGGCCCAGCCGCCCGACCT‑ACCGCAAAGCGACCCCGCCGGATGCGCCAATTAGGATCCTCACACTGACGGCCGATACTTATTCGCAGGGTGAACTGTCCGATTTCGCCTCGACGCGGCGGGCTCCGACGATGTCGCAAATCGACGGGGTTGGTGATGTGGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTgg                                                                       <  2:231102/251‑1 (MQ=255)                                                                                                                                                                             cAGTGGGATGCCCAGCCGCCCGTCCT‑TACGCAAAGCGAACCCGTCGGATGCCCCAATTAAGACCCCCACGCTGTCGTCCGATTCTTATTCGCAGGGTGACCTGTTCGATTTCGCCTCGACGCAGCTGGCCCCGACGTTTTCGCAAATCGACGGGGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGacgacg                                                                 <  2:273782/251‑1 (MQ=255)                                                                                                                                                                               gTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGacgacga                                                                >  2:115041/1‑249 (MQ=255)                                                                                                                                                                               gTGGGACGCCCAGCCGCCCGACCT‑ATCGCAAAGCGACCCCGCCGGATGCTCCAATTATGGTCCTACCCCTGACGTCCGATCCTTATTCGCAGGGTGAACTGGACGATTTCGCGTCGACGCTGCTGGCGCCGACGATTTAGCAAATCGCCGGTGTGGGGGAGGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTa                                                               <  2:173450/251‑1 (MQ=255)                                                                                                                                                                                   ggttcccAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGAGCCTCACACTGACGTCCTATACTTACTCGCAGGGTTAACTTTGCGATTTCGCCGCGACGCATCTGGCGCCGACTATTGCGCAAATCGACGGTGTTGGTGATGTCGAGGTCGTAGGCAGCTCACGGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCa                                                           <  1:374320/247‑1 (MQ=255)                                                                                                                                                                                            ccgcccgcCCT‑ATCGCAACACGAACCCGGCGGTTGCGCCAATTATTATCCCCACTCTGACGTCCGATACTTATTCCCTGGGTGAACTTGACGATTCCGCCTCGACGCAGCTGGCTCCGAGGATTTCGCATATCGACGGTGTTGTTGATGGCGATGTCGGAGGCAGCTCTCTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCGGGACGACGTACGCACCGCCGTCa                                                  <  2:406089/251‑1 (MQ=255)                                                                                                                                                                                             cgccCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTc                                                   >  1:209499/1‑249 (MQ=255)                                                                                                                                                                                                  gACCT‑ACCGCAAAGCGAACCCGTCGGAGGCGCCAATTATGTCCCTCCCGCTGCCGTCCGATCCTTATTCGCAGGGGGACCTGTACGATTTCGCCCCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATg                                            <  2:195085/251‑1 (MQ=255)                                                                                                                                                                                                          cGCAACGCGACCCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGt                                     >  2:290677/1‑251 (MQ=255)                                                                                                                                                                                                                 cGAACCCGTCGAATGCGCAAATTATGATCCTCACACGGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGTTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGACGCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCATCAGCAATGCCAACGTGCGTaaa                              <  1:21641/251‑1 (MQ=255)                                                                                                                                                                                                                          gggATGCGCCAATTAGGATCCTCACACGGACGTCCGATACTTATTCGCAGGGGGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGgc                     <  1:53310/250‑1 (MQ=255)                                                                                                                                                                                                                                cgcCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGCGCGCTg               >  1:605436/1‑251 (MQ=255)                                                                                                                                                                                                                                 gcCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGCGCGCTgg              >  1:432079/1‑251 (MQ=255)                                                                                                                                                                                                                                 gcCAATTATGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGCGCGCTgg              >  2:84316/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTGATCCTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGCGCGCTGGAAGATgg       >  1:175318/1‑251 (MQ=255)                                                                                                                                                                                                                                          gatccgctcGCTGACGTCCGTTATTTATTCGCAGGGTGAACTGTACTATTTCGCCTCGACGCAGCTGTCTCCGACGATTTCGAAAATCGACGGTGTTGGTGATTTCGATGTCGGAGGCAGCTCACGGCCTGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGCGCGCTGGAAGATGGCa     <  1:221685/243‑1 (MQ=255)                                                                                                                                                                                                                                            cccTCACACTGCCGTCGGGTAGTTATCCGCAGGTTGACCTTTACGATTTCGCGCCGACGCATTTGTCTCCGACGATTTCGCAAATGGACGGTGTTTGGGAGGTCGATGTCGGAGGAAGCTCACTGCCCGCCGATCGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGTTGGACGACGTACGCACCTCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGAGCGCTGGAAGATggttct   <  2:407335/250‑5 (MQ=255)                                                                                                                                                                                                                                              cTCACACTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGCGCGCTGGAAGATGGCACTc  >  1:88578/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     TGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGG‑AGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGGCGGCGATCAACGCTGCACAAAGTTTGCTGCCCAGTGGGATGCCCAGCCGCCCGACCT‑ATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAACCGCAGGGCGCGCTGGAAGATGGCACTC  >  NC_000913/2158574‑2159058

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