Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,332,791	0	A	C	27.3%	163.1 / 15.7	77	V161G (GTG→GGG)	basS	sensor histidine kinase BasS
Reads supporting (aligned to +/- strand):  ref base A (54/2);  new base C (2/19);  total (56/21)
Fisher's exact test for biased strand distribution p-value = 8.12e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCTGGCAATAGCAGGGTTTGCTGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTGATAGCAGATAAAT  >  NC_000913/4332553‑4333030                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                ctctGGCAATAGCAGGGTTTGCTGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTg                                                                                                                                                                                                                                     >  1:358452/1‑251 (MQ=255)   ctGGCAATAGCAGGGTTTGCTGTCGCTGGTCGAGCATGGAACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGc                                                                                                                                                                                                                                   >  1:467337/1‑251 (MQ=255)       aaaaagcggggTTTGCGGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCAAACGAGGGGAAAATAACATCTTCCAGCAGTTTTCCATGTTGATAATACCCGA‑AAAAAAATGACTGCCCGG‑CACGCGCCATTTGCAGCAGCGGGGAGACGCTCTCCATCATCGGACCAAGCCGTGCCACTACCGGTGTTAAATAAAT‑GTG‑A‑GGGTTTTCCCCCAGCAGTCCCAAATGAAACC‑GCCCCCCCGCCAGTGGCGTTc                                                                                                                                                                                                                               <  1:502190/243‑1 (MQ=255)        aaTAGCAGGGTTTGCTGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCg                                                                                                                                                                                                                              >  1:614546/1‑251 (MQ=255)         aTAGCAGGGTTTGCTGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGc                                                                                                                                                                                                                             >  1:476548/1‑251 (MQ=255)              gggttttGCTGTCGCGGGCCGAGCGGGGTACTG‑GTTCCGTCATAGGGGGGGAGAATGACATCTCCCAGCATTTTTACATGTTAATAATTACGGG‑AAAAAAATGACTGCCGGG‑CCCGCGCCAGTTGCAGCAGCGGGGAGCCGCTCTCCATCATCTGACCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGAAAAC‑CCACCCCCGCCAGTGGCGTTCGCAGTTc                                                                                                                                                                                                                        <  1:77410/247‑1 (MQ=255)              aGGGTTTGCTGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTc                                                                                                                                                                                                                        >  2:40026/1‑251 (MQ=255)               gggTTTGCTGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                                                                                                                       >  1:88555/1‑251 (MQ=255)                 gTTTGCTGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GGG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCCCCCCCGCCAGTGGCGTTCGAAGTTCgcg                                                                                                                                                                                                                     >  2:411793/1‑249 (MQ=255)                       tGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGAc                                                                                                                                                                                                               >  1:24080/1‑251 (MQ=255)                          cGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTc                                                                                                                                                                                                            >  1:593397/1‑251 (MQ=255)                              ggTCGAACATGGTACTGCGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGATTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCAGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGGTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCg                                                                                                                                                                                                        >  2:590047/1‑251 (MQ=255)                                    gCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaac                                                                                                                                                                                                  >  1:618112/1‑251 (MQ=255)                                      aTGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                                                                                                                                 >  1:458128/1‑250 (MQ=255)                                      aTGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                                                                                                                                 >  2:320605/1‑250 (MQ=255)                                         gTACTGAGTTC‑GTAATACGGGGGGAAATTGACACCTCCCAGAATTTTTACTTGTTAATAATTACGGA‑AAGAAAATGCCTGTCCGG‑CACGCGCCAGTTGCAGCGGCGGGGAGCCCCTCCCCTTCCTCTGATAAAGCCGTGCCACAACCGGTGCTACATCAAT‑GTG‑ATGCTTTTCCGCC‑AGCAGTCCCAGATGCAACC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCt                                                                                                                                                                                             <  1:546071/251‑1 (MQ=255)                                           aCTGAGTTC‑TTCATACGGGGGGAGAATGCCATCTTCAAGAAGTTTTCCATGTGAATAATCACCGA‑AAAAAAATGCCTGCCCGG‑CCCGCGCAATTGGCAGCAGCTGGGAGCCGCTCTCCATCATCTGACCAAGCCGTGCCCATACGGGTGCTACATAAAT‑GTG‑AGGCTTTTTCGCC‑AGCAGTCCCAGAGGAAAAC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                                                                                                                           <  1:490147/251‑1 (MQ=255)                                                  tC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATc                                                                                                                                                                                    >  1:579512/1‑251 (MQ=255)                                                       aaaagggggggAGAATGACATTTTCCAGCATTTTTACTTTTTGAAATTTACGGG‑AAGAAAATGACTGTCGGG‑CACGCGCAAGTTGCAGCGGGGGGGAGCCTCCCTCCTTCATCTGATCAAGCCGTCCCCCTAAGGGTGCTCCACAAAT‑GTG‑AGGCTTTTCCCCC‑AGGAGTCCCAGATGCAAAC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                                                                                                                <  1:533190/244‑1 (MQ=255)                                                        aTACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGTGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                                                                                                                               >  2:216417/1‑251 (MQ=255)                                                        aTACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                                                                                                                               >  2:507196/1‑251 (MQ=255)                                                            gAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCt                                                                                                                                                                           >  1:258632/1‑251 (MQ=255)                                                                ggAGAATGACATCTTCCAGCAGTTTTACATGATGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                                                                                                                       >  1:637103/1‑251 (MQ=255)                                                                    aaggccctctcccAGCATTTTTACTTTTTGATAATCACGGG‑AAAAAAATGACTGTCCGG‑CCCGCGCCAGTTGCAGCAGCTGGGAGACCCTCCCCATCTTCTGATCAAGCCGGGCCACTACCGGTGCTACA‑CAATTGTG‑ATGCTTTTTCGCC‑AGCAGTTCCAGATGCAACC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                                                                                                                   <  1:371693/240‑1 (MQ=255)                                                                      tGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCTTGCTGGTCAGGCg                                                                                                                                                                  >  2:61008/1‑250 (MQ=255)                                                                       gACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGAc                                                                                                                                                                >  2:147039/1‑251 (MQ=255)                                                                        aCATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGACGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACATTTCGTTATCCAGCGTGCTGGTCAGGCGACt                                                                                                                                                               >  2:84713/1‑251 (MQ=255)                                                                         cATCTCCCAGCATTTTTACATGTTGA‑AATTTCCGGGAAAAAAATGACTGCCCGG‑CCCGCGCAATGTGCACCGGCTGGGAGACCCTCTCCATCTTCTGATCAGGCCGTGCCATTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTCCCAGATGAAAAC‑CCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTg                                                                                                                                                              <  1:16499/251‑1 (MQ=255)                                                                                   cagTTTTACATGTTGCTAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGt                                                                                                                                                    >  2:553891/1‑251 (MQ=255)                                                                                   cagTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGt                                                                                                                                                    >  1:289681/1‑251 (MQ=255)                                                                                      ttttACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAA‑GCACCCCCGCCAGTGGCGtt                                                                                                                                                                                                                                >  2:412692/1‑172 (MQ=255)                                                                                          aCATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGc                                                                                                                                             >  1:397190/1‑251 (MQ=255)                                                                                             tGTTAATAATTCCCGG‑AAGAAAAGGACTGCCCGG‑CCCGCGCCAGTTGCAGCAGCGGGGAGACGCCCCCCATCTTCTGCTCAACCCGTCCCACTATCGGTGCCACAAAAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCGAGAGCCAAAC‑CCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCg                                                                                                                                           <  1:29883/250‑1 (MQ=255)                                                                                              tttGACAATTACGGG‑AAAAAAATGATTGTCCGG‑CCCGCGCCAGTGGCAGCGGCGGGGAGACGCTCCCCATCATCTGACCAAGCCGTGCCACTAGGGGTGCAAAACTAAG‑GTG‑AGGCTTTTCCGCC‑AGCAGTTCCAGATAAAACC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGa                                                                                                                                          <  1:514929/249‑1 (MQ=255)                                                                                                    aaTTACCGG‑AAGAAAAGGACTGCCCGG‑CCCGCGCCAGTGGCAGCAGCGGGGAGCCCCCCTCCATCATCTAATCAAGCCGTGCCACTAACGGTGCTACATAAAT‑GTG‑ATGCGTTTCCCCC‑AGCAGTCCCAGATGCAAAC‑CCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAaccacc                                                                                                                                   <  1:229828/251‑1 (MQ=255)                                                                                                       tACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCt                                                                                                                                >  2:56363/1‑251 (MQ=255)                                                                                                         ccGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCCACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGCTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGACAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTCACCAGATCGTTTAACCTCGAAACCACCGCTTc                                                                                                                              >  1:568474/1‑251 (MQ=255)                                                                                                           gg‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCgt                                                                                                                            >  1:337659/1‑250 (MQ=255)                                                                                                              aagaaAATGACTGTCCGGGCA‑GCGCCAGTTGCAGCCGATGGGAGACGGACTCCCTCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTC‑ATGCGTTTTCGCC‑AGAAGTTCCAGATGCAAAC‑GCACCCTCGCCCGTGGCGGTCGCAATTCGTGCGCGACGTCAGCGGTCAACAGCCATTCGTCATCCGGCGTGCTGGTCAGGCGACTGACCAGTCCGGTTGACGCCGAAACCCCCGCTTCCAtt                                                                                                                          >  2:261872/1‑251 (MQ=255)                                                                                                              aagaaAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAtt                                                                                                                          >  1:481668/1‑251 (MQ=255)                                                                                                                 aaaattgccggcccgccccgcgCCAGTTGCGGCGGCGGGGAGCCCCCCCCCATCATCGGACCAAGCCGTGCCACTACGGGTGCTACACCAAT‑GGG‑TTGCGTTTCCGCC‑AGCATTCCCAGATGCAACC‑CCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc                                                                                                                        <  1:507196/233‑1 (MQ=255)                                                                                                                 aaaaTGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg                                                                                                                       >  2:409346/1‑251 (MQ=255)                                                                                                                    aTGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCgtgg                                                                                                                    >  2:258148/1‑248 (MQ=255)                                                                                                                       aCTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTg                                                                                                                 >  1:61494/1‑251 (MQ=255)                                                                                                                       aCTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTg                                                                                                                 >  1:592038/1‑251 (MQ=255)                                                                                                                              ctgaaCGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATACGTGTTCGCC‑AGCAGTTACAGATGCCAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCAACGTCAGCGGTAAACAACATTTCGTTCTCACGCGTGCTGCTACGGCGACTGAACAGATCGTTTCACGACGACACCACCGCTTCGATTTCGAGGGAGGCGCtg                                                                                                           >  2:244008/5‑251 (MQ=255)                                                                                                                               g‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                                                                                                         >  1:33833/1‑251 (MQ=255)                                                                                                                               g‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCtgtg                                                                                                         >  1:236577/1‑251 (MQ=255)                                                                                                                               g‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GAG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGCGCTGGTCAGGCGACTGACCAGACCGCTTAACGCCGAAACCACCGCTTCGATTTCGAGGGGGGCGCtgtg                                                                                                         >  1:311180/1‑251 (MQ=255)                                                                                                                                     cgcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCCAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGAAGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGCCTGACCAGATCGTTTAACGCCGAAAACACCGCTTCGATTTCGAGGGTGGCGCTGTGAATgg                                                                                                    >  2:145693/1‑251 (MQ=255)                                                                                                                                       ccAGTTGCAGCAGCTGGGAGGCGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCAGGTCAGGCGCCTGACCAGATCGTTTAACGCCGAACCCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                                                                                   >  2:196128/1‑250 (MQ=255)                                                                                                                                         aGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAAt                                                                                                >  1:541853/1‑251 (MQ=255)                                                                                                                                                     ggggggccactccccatcttcgtaccAAGCGGTGCCACTAACGGTGCAACATAAAG‑GTT‑ATGTGTTTCCGCC‑AGCAGTCCCAGATGCAACC‑CCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGtt                                                                                    <  2:458128/226‑1 (MQ=255)                                                                                                                                                      gggAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTg                                                                                   >  2:367989/1‑251 (MQ=255)                                                                                                                                                        gagaCGCTCCCCATCTTCTAACCAGGCCGTGCAACTAACGGTGCTACATCAAT‑GTG‑GTGCTTTTCCGCC‑AGCATTTCCAGATGAAAAC‑CCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTc                                                                                 <  1:97374/251‑1 (MQ=255)                                                                                                                                                               tctcCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTg                                                                          >  1:508416/1‑251 (MQ=255)                                                                                                                                                                 cccctcgtcTGAACAAGCCGTCCCACTAGCGGTGCTACATCAAT‑GGG‑GTGCTTTTCCGCC‑AGCAGTCCCAGATGAAAAC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcg                                                                        <  1:367989/244‑1 (MQ=255)                                                                                                                                                                  ccATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgc                                                                       >  1:305552/1‑251 (MQ=255)                                                                                                                                                                    atcatcTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGCAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTgcgcgc                                                                     >  2:13961/1‑251 (MQ=255)                                                                                                                                                                         cTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTcc                                                                 >  1:5602/1‑250 (MQ=255)                                                                                                                                                                               aaGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCtt                                                          >  1:69453/1‑251 (MQ=255)                                                                                                                                                                               aaGCCGGGCCACTAGGGGTGTTACATCAA‑‑GTGGATGCTTTTTCGCC‑AGCGGTTCTAGAGGCAAAC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCtt                                                          <  1:411793/251‑1 (MQ=255)                                                                                                                                                                                           tagcGGTGCTACATCAAT‑GGG‑ATGCGTTTTCCCC‑AGCAGTTCGAGA‑GCAAACCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTcc                                               <  1:553561/247‑1 (MQ=255)                                                                                                                                                                                           tAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCg                                              >  1:159728/1‑251 (MQ=255)                                                                                                                                                                                               ggggtgttactcaagtgtgttggtttttCCCC‑AGCATTTCCAGAGGAAACC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGcc                                            <  1:503848/227‑1 (MQ=255)                                                                                                                                                                                               ggTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGcc                                            >  1:203466/1‑249 (MQ=255)                                                                                                                                                                                               ggTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACAGCTTCGATTTCGCGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTGTTGCAGCCCCGcc                                            >  2:430090/1‑249 (MQ=255)                                                                                                                                                                                                  gcaacaccAAT‑GTG‑ATGTGTTTCCCCC‑AGCAGTCCCAGAGGAAAAC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAgcgg                                       <  1:155491/244‑1 (MQ=255)                                                                                                                                                                                                          aaT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCC‑GCCAGTGGCGTTCGCCGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAACCAGCGTGCTGGTCAGGCGACTGGACAGGTCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGc                                                                                                   >  2:557775/1‑182 (MQ=255)                                                                                                                                                                                                            t‑gtg‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATg                             >  1:84253/1‑251 (MQ=255)                                                                                                                                                                                                                    gCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGt                       >  1:314971/1‑251 (MQ=255)                                                                                                                                                                                                                        ttgcGCC‑AGCAGTCCCAGATGCACAC‑CCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTAcc                    <  2:113766/247‑1 (MQ=255)                                                                                                                                                                                                                           cGCC‑AGCAGTCCCAGACGCAACC‑GCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGTGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCt                <  1:216417/251‑1 (MQ=255)                                                                                                                                                                                                                             cc‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTgc              >  2:37221/1‑250 (MQ=255)                                                                                                                                                                                                                                        aaGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTGATAGCAGATaa    <  1:452492/250‑1 (MQ=255)                                                                                                                                                                                                                                           aTGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGCGGGCGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCCGCTCCGCCCGCGGGCGGGTGGTGCGGCGTACCGCCTggttgacggtaaat  >  2:29303/1‑237 (MQ=255)                                                                                                                                                                                                                                           aTGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTGATAGCCGATAAAt  >  2:258171/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                CTCTGGCAATAGCAGGGTTTGCTGTCGCTGGTCGAGCATGGTACTGAGTTC‑GTCATACGAGGGGAGAATGACATCTTCCAGCAGTTTTACATGTTGATAATTACCGG‑AAGAAAATGACTGTCCGG‑CACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAAT‑GTG‑ATGCGTTTTCGCC‑AGCAGTTCCAGATGCAAAC‑GCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCGTTAAGTTGTCGGCGGTGCGCGCTTCCAGCTCTTTTTGCAGCTCCGCCAGCGGGCGGGTGATGCGGCGTACCGCCTGATAGCAGATAAAT  >  NC_000913/4332553‑4333030

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