Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,530,252 C→A intergenic (‑9/+3) sgcC ← / ← sgcB putative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,2520CA97.6% 252.6 / NA 84intergenic (‑9/+3)sgcC/sgcBputative PTS enzyme IIC component SgcC/putative PTS enzyme IIB component SgcB
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (36/45);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

CCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCG  >  NC_000913/4530005‑4530496
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
cgcgcgcCAACGTCAAAGACAGGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGAACTGAGGCGGTCGACCGCCATGACGATCACCAGTCCCATACCGACAAAACCAACCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCAGGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAttt                                                                                                                                                                                                                                                   <  2:547354/248‑1 (MQ=255)
ccTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAttt                                                                                                                                                                                                                                                   >  1:286715/1‑251 (MQ=255)
          cgacaaaaaaaTGCAGGGCTAGGCCGAACCGCTCAATCACCACTTTGATGGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAACCCAACCCCAACCGTTA‑CCCCTGCCTTAATCGCCGGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTa                                                                                                                                                                                                                                         <  2:350472/241‑1 (MQ=255)
            ccAAAAACATGCAGGGCTAGGCCGAAACGCTAAATCATCCCTTTGATAGGGGGATTGAGGCTGTCGATCGCCATGCCGATCACCAGTCCCATCCCGACAAAACCAACCCCACCCGTTC‑CCCCGGCCTTAATCGCCTGCAGCCAGGGGATGCGAAAGAAAATGCCGATCACGATCATGATGATAGGAACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAAc                                                                                                                                                                                                                                       <  1:589563/250‑1 (MQ=255)
             cAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCACCCCTTTGATAGGGGGACTGAGGCTGTCGATCGCCAGGACGATCCCCAGTCCCATACCGACAAAACCAACCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACa                                                                                                                                                                                                                                      <  2:460133/251‑1 (MQ=255)
               aaGACATGCAGGGCTAGGCCGAAACGCTCAATCCTCCCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAAGCCAATCAAGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAACCATAATGACTCATTTTATTGAGTTAAc                                                                                                                                                                                                                                       >  2:53686/1‑248 (MQ=255)
                 aaCATGCCGGGCTGGGCAGAAACGCACAATCATCACTTTGATAGGGGGACTGAGGCTGTAGATCGCCATGACGAGCACCATCCCCATACCGACAAAACCAATCCCACCCGTTA‑CCCCTGCCTTAATCCCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCAAGATGATAGGGACAAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGc                                                                                                                                                                                                                                  <  1:520206/250‑1 (MQ=255)
                    aTGCAAGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCtt                                                                                                                                                                                                                               >  1:195008/1‑251 (MQ=255)
                      gCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTg                                                                                                                                                                                                                              >  2:1844/1‑250 (MQ=255)
                            cTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATttgtt                                                                                                                                                                                                                       >  2:516541/1‑251 (MQ=255)
                              aGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTtgttgt                                                                                                                                                                                                                     >  1:41204/1‑251 (MQ=255)
                                gctgAAACGCTAAAACATCACATTGTGAGTGGGACTGAGTCTGTCGTTCGCCAGGACGACCCCCAGTCCCATACCTACAAAACCAATCCCAACCGTCA‑CCCCGGCCTTAATGTCCTGCAGCCAGGGGATGCGGAGGATAAGGCAAATAACTATAATGAAGAAAGGGACGAAAACGTT‑CCCTCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTGATTGATTTAACAGACCCTTGAtgttttgttt                                                                                                                                                                                                                   <  2:357802/248‑1 (MQ=255)
                                 ccGAAACGCGCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGACCACCAGTCCCATACCGACAAAACCAATCCCACCCGTTA‑CCCCGGCCTTAATCGCCCGCAGCCAGGGGATGCGGAATATAAGGCCAATCACGATGATGATTATAGGGACGAAAACGGT‑GCCGCCGATAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGtttt                                                                                                                                                                                                                  <  1:590108/251‑1 (MQ=255)
                                              tcatcaCTTTGATAGGGGGACTGAGGCTGTCGATCGACATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGCTGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGGTTTGTTGTTTTAACGcgtcgtcgt                                                                                                                                                                                                     >  2:501607/1‑251 (MQ=255)
                                               catcaCTTTGATAGGGGGACTGAGGCTGCCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGtt                                                                                                                                                                                                    >  2:489802/1‑251 (MQ=255)
                                                      ttGAAAGGGGGACGGAGGCTGTCGATCGCCAGGACGCCCCCCAGTCCCATACCGACAAAACAAATCCCAACCTTTA‑CCCCTGCCTTAATCGCCGGCAGCCAGGGGAGGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATcc                                                                                                                                                                                               <  2:128315/249‑1 (MQ=255)
                                                       tGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCACAACCGTTA‑CCCCTGCCTTAATCGCCTGCCGCCAGGGGAAGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCtg                                                                                                                                                                                            >  2:19200/1‑251 (MQ=255)
                                                                    gAGGCTGTCTATCCCCATGCAAATCACCGGACCCATACCGACAAACCCAATCCCAGCCGTTC‑CCCCTGCCTAAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACAATCAGGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGggcggc                                                                                                                                                                               <  1:351765/251‑1 (MQ=255)
                                                                          gTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCAtt                                                                                                                                                                         >  1:145870/1‑251 (MQ=255)
                                                                             gATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTa                                                                                                                                                                        >  1:59380/1‑249 (MQ=255)
                                                                              atcgcccagacgatgactcgTCCTATACCGACAAAACCAATCCCAAACGTTA‑CCCCTGCCTTACTCGCCTGAAGCCAGGGGATGCGGAAGATAAGGCCAATCAAGATAATGCTGATAGGGACGCAAACGGG‑ACCGCCGAGAGACAGGCGGTAATCGAACATAATGCCTAGTTTTATTGAGTTAACAGCGCCTTGATTTGTTGGTTTAACGCGTCGTCGTTGATACCTGTGAGCGGGGCGGCGCCATTAAGc                                                                                                                                                                     >  2:268515/20‑251 (MQ=255)
                                                                                cgccataccGATCACCCGTCCCTTACCGACAAAACAAACCCCAACCGTTC‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAAAAGGCCAATCACGATCAGGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGt                                                                                                                                                                   <  1:501607/243‑1 (MQ=255)
                                                                                cGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGt                                                                                                                                                                   >  1:494971/1‑251 (MQ=255)
                                                                                  ccaggACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAg                                                                                                                                                                 <  1:260134/247‑1 (MQ=255)
                                                                                   caaaaCGAGCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CAACTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGCTCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCGACCATATTGCCTCATTTTATAGAGTTAACAGGGCCTTGATTTTTTGTTTTAACGCGTCGTCGTTGATCCATGTGCGCAGGGCGGCGCCATCAAGCgatgg                                                                                                                                                                >  2:434871/5‑247 (MQ=255)
                                                                                      gacaatccccatgcccAGACCGACAAAACCATTCCCAACCGTCC‑CCCCCGCCTTAATCGACTGCAGCCAGGGGATGCGGAAGAAAAGGCCAATCACGTTCATGATGAAAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCGTGATTTGTTGTTTTAACGCGTCGTCGTTTATCCCTGTGAGCAGTGCGGCGCCATTAAGCGTAGGGAt                                                                                                                                                             <  2:277153/238‑1 (MQ=255)
                                                                                      gACGATCACCAGTCCCCTACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTACTCGCCTGCCGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGACGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGAAGGGCGGCGCCATTAAGCGTAGGGAt                                                                                                                                                             >  2:238052/1‑251 (MQ=255)
                                                                                      gACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGAt                                                                                                                                                             >  1:603895/1‑251 (MQ=255)
                                                                                          atccccAGTCCCATACCGACAAAACCAACCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCAGGATGATAGGGACGAAACCGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCg                                                                                                                                                         <  1:336661/247‑1 (MQ=255)
                                                                                            cACCAGTCCCAAACCGACAAAACCAATCCCAACCGTTC‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCATGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTa                                                                                                                                                       <  1:390660/251‑1 (MQ=255)
                                                                                              ccAGCCCCATACCAACAAAACAAACCCCACCCGTTA‑CCCCTGCCGAAACCGCCGGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAAt                                                                                                                                                     <  1:542842/251‑1 (MQ=255)
                                                                                               cAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATc                                                                                                                                                    >  1:618138/1‑251 (MQ=255)
                                                                                                   cccATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGCTGCCGTAATCACTg                                                                                                                                                >  1:26209/1‑251 (MQ=255)
                                                                                                   cccATACCGACAAAACCAATCCCAACCGTTA‑CCACTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTACCGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTg                                                                                                                                                >  1:199082/1‑251 (MQ=255)
                                                                                                    ccATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGt                                                                                                                                               >  1:247994/1‑251 (MQ=255)
                                                                                                         ccGACAAAACCAACCCCACCCGTTC‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGtt                                                                                                                                          <  1:52444/251‑1 (MQ=255)
                                                                                                          cGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAAGCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTc                                                                                                                                         <  1:537516/251‑1 (MQ=255)
                                                                                                             cAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGc                                                                                                                                       >  2:197253/1‑250 (MQ=255)
                                                                                                             cAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGc                                                                                                                                       >  2:192251/1‑250 (MQ=255)
                                                                                                              aaaaCCAATCCCAACC‑TTAGCCCCGGCCTTAATCGCCGGCAGCCAGGGGATGCGAAATATAAGGCCAATCAGGATCAGTAGAATGGGAACGAAAACGTG‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTTAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCAt                                                                                                                                     <  1:268515/251‑1 (MQ=255)
                                                                                                                aaCCAATCCCAACCGTTA‑CCCCTGACTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATgg                                                                                                                                   >  2:132456/1‑251 (MQ=255)
                                                                                                                 accacccccAACCGTCA‑CCCCTGCCATAATCGCGTGGAGCCAGGGGATGCGGAAAATAAGGCCAATCACGATCATGATGAGAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAATAGAGCCTTGAGTTGTTGTTTTAACGCGTCGTCGTTGAGCCCTTTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTTGCATGGa                                                                                                                                  <  2:253840/245‑1 (MQ=255)
                                                                                                                        ccAACCGGTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGAAGCGGAAGATAAGACCAACCACGATCATGATTATAGGGAAGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCTTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGAAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACa                                                                                                                           <  1:394156/251‑1 (MQ=255)
                                                                                                                           aCCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAt                                                                                                                         >  2:16326/1‑250 (MQ=255)
                                                                                                                             cgtaa‑CCCCTGCCTTAATCGCCTGCAGCAAGGGGATGCGGAAGATAAGGCCAAACACGATAATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAg                                                                                                                      <  1:117907/247‑1 (MQ=255)
                                                                                                                                 a‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTc                                                                                                                   <  1:516541/250‑1 (MQ=255)
                                                                                                                                   ccccGGCCTTAATCGCCTGCAGCCATGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑TCTGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATGTTATTGAGTTAACAGAGCCGTGATTTGTTGTTTGAACGCGTCGTCGTTTATCCCTGTGAGCAGGGCGGCTCCAGTAATCGTAGTGATGCCGTAATCACAGTTGTTTCGCAGGGATGTTGCAATAAGGACCa                                                                                                                 <  2:40766/251‑1 (MQ=255)
                                                                                                                                          cTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGtt                                                                                                          >  2:8424/1‑251 (MQ=255)
                                                                                                                                                      cAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGCGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTCACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAgg                                                                                              >  1:392740/1‑251 (MQ=255)
                                                                                                                                                      cAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCATCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAgg                                                                                              >  2:210782/1‑251 (MQ=255)
                                                                                                                                                         ccAGGGGATGCGGAAGATAAGGCCAATCCCGATAATGAGGAAAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAt                                                                                           <  1:61503/251‑1 (MQ=255)
                                                                                                                                                         ccAGGGGATGCGGAAGATAAGGCAAATCCCGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGGATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAt                                                                                           <  2:247994/251‑1 (MQ=255)
                                                                                                                                                            ggggATGCGGAAGAAAAGACAAATCACGATCATGATGATAGGGCCGAAAACGTT‑ACCGCCGAGATACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGTCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAGAAGGTCCTTATCGTTACAGTTAATAGGGAtctc                                                                                        <  2:561966/251‑1 (MQ=255)
                                                                                                                                                                  gCGGAAGATAAGCCAAACCACGATAAGGAGGTTAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATTCCTCAGTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACGGTTTGGTCGCATGGAGGTTACAATAAGGTCCATTCCGTTACGGTTAAGAGGGATCTCATTCAg                                                                                  <  2:626367/251‑1 (MQ=255)
                                                                                                                                                                     gAAGATAAGGCCAATCACGACCATGATGATAGGGACGAAAACGTT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACa                                                                               <  1:1336/251‑1 (MQ=255)
                                                                                                                                                                          tAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCAtt                                                                          >  1:283704/1‑251 (MQ=255)
                                                                                                                                                                          tAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCAtt                                                                          >  1:115812/1‑251 (MQ=255)
                                                                                                                                                                           aaGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTg                                                                         >  1:580264/1‑251 (MQ=255)
                                                                                                                                                                            ggaccaatcacaagcagaaggatAGGGAAGAAAA‑GGTGACCGGCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGGAGGGATGCCGTAATCACTGTTGTTTGGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGGCAGCATTgg                                                                        <  2:207964/231‑1 (MQ=255)
                                                                                                                                                                            atccgcatcccaaTCATGATGATAGGAACGAAATCTGT‑ACCTCCGAGAGACATGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGTGCAGGGCGGCGCCATTAAGCGTAGGTAGTCCGTAATCACTGTTTGTTCGCATGGAGGTTTCAATATGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTgg                                                                        <  1:81890/239‑1 (MQ=255)
                                                                                                                                                                             ggCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTggg                                                                       <  1:129109/251‑1 (MQ=255)
                                                                                                                                                                                  aTCACGATAATGATGATAGGGAAGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTa                                                                  <  1:529575/251‑1 (MQ=255)
                                                                                                                                                                                       gATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGc                                                             >  1:14077/1‑251 (MQ=255)
                                                                                                                                                                                        ataaggatgatAGGGACGAAAAAGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGGTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCt                                                            <  2:614569/246‑1 (MQ=255)
                                                                                                                                                                                           atgatgATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGa                                                          >  1:619348/1‑250 (MQ=255)
                                                                                                                                                                                              atgatAGGGACGCAAACAAT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGACTCCTTTTATTGAGTTAACAGCGCCTTGATTTGTTGTTTTAACGCCTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGAGGCCGTAATCACTGCTGGTTCGCTTGGAGGTTACCATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGACGGTAGTGGCTGAAATg                                                      >  2:296484/1‑251 (MQ=255)
                                                                                                                                                                                                 atAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGccc                                                   <  1:53686/251‑1 (MQ=255)
                                                                                                                                                                                                   aggtaCGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACCTAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTGCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTg                                                 <  1:244233/247‑1 (MQ=255)
                                                                                                                                                                                                     ggACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGtt                                               <  1:8424/251‑1 (MQ=255)
                                                                                                                                                                                                       aataaatcGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtca                                             <  2:4269/244‑1 (MQ=255)
                                                                                                                                                                                                              cGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAgg                                       >  1:143705/1‑250 (MQ=255)
                                                                                                                                                                                                                        cgAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCa                             <  1:98061/251‑1 (MQ=255)
                                                                                                                                                                                                                         gagagaCAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGc                              >  1:324878/1‑249 (MQ=255)
                                                                                                                                                                                                                              aCAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgt                       <  1:252017/251‑1 (MQ=255)
                                                                                                                                                                                                                               cAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgtg                      >  1:29351/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGGCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAGAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTTCAGCTTGTGTGCAAt                <  2:41204/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaTCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCGTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGGTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATc               <  2:195008/250‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaCATAATGCATCATTTTATTGGGTTAACAGAGCCATGCACTGATGCTCTAACGCGTCGTCGCTGATCCCTGTGAGCAGGCCGGCGACCTTAAGCGTACGGATGCCGTAATCAATGTTGGTTCGCATGGAGGTTACAATACGATCCATACCGTTAAAGATAAGAGGGATCTCATACAGACAGCATTGGGCGGTAGTGGCTGAAATGCACTGTTCAGTTAGGAACTCTTGCCGCTAGTGTGCAATCATGGt          >  2:364112/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaCATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGtt         >  2:458692/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aCATAATGCCTCATTTTATTGAGTTAACAGGGCCGTGATTTGTTGTTTTAACGCGTCGTCGTTGATACCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTg        <  2:119070/251‑1 (MQ=255)
                                                                                                                                                                                                                                             aCATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCGCTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGATGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATTCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTg        <  1:419048/251‑1 (MQ=255)
                                                                                                                                                                                                                                              cATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTg        >  1:168605/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  aTGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTc   <  1:56037/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    gCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGAGCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCg  <  2:145870/250‑1 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
CCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTA‑CCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGT‑ACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCG  >  NC_000913/4530005‑4530496

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