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

CCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCG  >  NC_000913/4530005‑4530496
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
cgcgcgcCAACGTCAAAGACAGGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGAACTGAGGCGGTCGACCGCCATGACGATCACCAGTCCCATACCGACAAAACCAACCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCAGGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAttt                                                                                                                                                                                                                                                   <  2:547354/248‑1 (MQ=255)
ccTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCAttt                                                                                                                                                                                                                                                   >  1:286715/1‑251 (MQ=255)
          cgacaaaaaaaTGCAGGGCTAGGCCGAACCGCTCAATCACCACTTTGATGGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAACCCAACCCCAACCGTTACCCCTGCCTTAATCGCCGGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTa                                                                                                                                                                                                                                         <  2:350472/241‑1 (MQ=255)
            ccAAAAACATGCAGGGCTAGGCCGAAACGCTAAATCATCCCTTTGATAGGGGGATTGAGGCTGTCGATCGCCATGCCGATCACCAGTCCCATCCCGACAAAACCAACCCCACCCGTTCCCCCGGCCTTAATCGCCTGCAGCCAGGGGATGCGAAAGAAAATGCCGATCACGATCATGATGATAGGAACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAAc                                                                                                                                                                                                                                       <  1:589563/250‑1 (MQ=255)
             cAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCACCCCTTTGATAGGGGGACTGAGGCTGTCGATCGCCAGGACGATCCCCAGTCCCATACCGACAAAACCAACCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACa                                                                                                                                                                                                                                      <  2:460133/251‑1 (MQ=255)
               aaGACATGCAGGGCTAGGCCGAAACGCTCAATCCTCCCTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAAGCCAATCAAGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAACCATAATGACTCATTTTATTGAGTTAAc                                                                                                                                                                                                                                       >  2:53686/1‑248 (MQ=255)
                 aaCATGCCGGGCTGGGCAGAAACGCACAATCATCACTTTGATAGGGGGACTGAGGCTGTAGATCGCCATGACGAGCACCATCCCCATACCGACAAAACCAATCCCACCCGTTACCCCTGCCTTAATCCCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCAAGATGATAGGGACAAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGc                                                                                                                                                                                                                                  <  1:520206/250‑1 (MQ=255)
                    aTGCAAGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCtt                                                                                                                                                                                                                               >  1:195008/1‑251 (MQ=255)
                      gCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTg                                                                                                                                                                                                                              >  2:1844/1‑250 (MQ=255)
                            cTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATttgtt                                                                                                                                                                                                                       >  2:516541/1‑251 (MQ=255)
                              aGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTtgttgt                                                                                                                                                                                                                     >  1:41204/1‑251 (MQ=255)
                                gctgAAACGCTAAAACATCACATTGTGAGTGGGACTGAGTCTGTCGTTCGCCAGGACGACCCCCAGTCCCATACCTACAAAACCAATCCCAACCGTCACCCCGGCCTTAATGTCCTGCAGCCAGGGGATGCGGAGGATAAGGCAAATAACTATAATGAAGAAAGGGACGAAAACGTTCCCTCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTGATTGATTTAACAGACCCTTGAtgttttgttt                                                                                                                                                                                                                   <  2:357802/248‑1 (MQ=255)
                                 ccGAAACGCGCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGACCACCAGTCCCATACCGACAAAACCAATCCCACCCGTTACCCCGGCCTTAATCGCCCGCAGCCAGGGGATGCGGAATATAAGGCCAATCACGATGATGATTATAGGGACGAAAACGGTGCCGCCGATAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGtttt                                                                                                                                                                                                                  <  1:590108/251‑1 (MQ=255)
                                              tcatcaCTTTGATAGGGGGACTGAGGCTGTCGATCGACATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGCTGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGGTTTGTTGTTTTAACGcgtcgtcgt                                                                                                                                                                                                     >  2:501607/1‑251 (MQ=255)
                                               catcaCTTTGATAGGGGGACTGAGGCTGCCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGtt                                                                                                                                                                                                    >  2:489802/1‑251 (MQ=255)
                                                      ttGAAAGGGGGACGGAGGCTGTCGATCGCCAGGACGCCCCCCAGTCCCATACCGACAAAACAAATCCCAACCTTTACCCCTGCCTTAATCGCCGGCAGCCAGGGGAGGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATcc                                                                                                                                                                                               <  2:128315/249‑1 (MQ=255)
                                                       tGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCACAACCGTTACCCCTGCCTTAATCGCCTGCCGCCAGGGGAAGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCtg                                                                                                                                                                                            >  2:19200/1‑251 (MQ=255)
                                                                    gAGGCTGTCTATCCCCATGCAAATCACCGGACCCATACCGACAAACCCAATCCCAGCCGTTCCCCCTGCCTAAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACAATCAGGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGggcggc                                                                                                                                                                               <  1:351765/251‑1 (MQ=255)
                                                                          gTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCAtt                                                                                                                                                                         >  1:145870/1‑251 (MQ=255)
                                                                             gATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTa                                                                                                                                                                        >  1:59380/1‑249 (MQ=255)
                                                                              atcgcccagacgatgactcgTCCTATACCGACAAAACCAATCCCAAACGTTACCCCTGCCTTACTCGCCTGAAGCCAGGGGATGCGGAAGATAAGGCCAATCAAGATAATGCTGATAGGGACGCAAACGGGACCGCCGAGAGACAGGCGGTAATCGAACATAATGCCTAGTTTTATTGAGTTAACAGCGCCTTGATTTGTTGGTTTAACGCGTCGTCGTTGATACCTGTGAGCGGGGCGGCGCCATTAAGc                                                                                                                                                                     >  2:268515/20‑251 (MQ=255)
                                                                                cgccataccGATCACCCGTCCCTTACCGACAAAACAAACCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGAAAAGGCCAATCACGATCAGGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGt                                                                                                                                                                   <  1:501607/243‑1 (MQ=255)
                                                                                cGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGt                                                                                                                                                                   >  1:494971/1‑251 (MQ=255)
                                                                                  ccaggACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAg                                                                                                                                                                 <  1:260134/247‑1 (MQ=255)
                                                                                   caaaaCGAGCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACAACTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGCTCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCGACCATATTGCCTCATTTTATAGAGTTAACAGGGCCTTGATTTTTTGTTTTAACGCGTCGTCGTTGATCCATGTGCGCAGGGCGGCGCCATCAAGCgatgg                                                                                                                                                                >  2:434871/5‑247 (MQ=255)
                                                                                      gacaatccccatgcccAGACCGACAAAACCATTCCCAACCGTCCCCCCCGCCTTAATCGACTGCAGCCAGGGGATGCGGAAGAAAAGGCCAATCACGTTCATGATGAAAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCGTGATTTGTTGTTTTAACGCGTCGTCGTTTATCCCTGTGAGCAGTGCGGCGCCATTAAGCGTAGGGAt                                                                                                                                                             <  2:277153/238‑1 (MQ=255)
                                                                                      gACGATCACCAGTCCCCTACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTACTCGCCTGCCGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGACGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGAAGGGCGGCGCCATTAAGCGTAGGGAt                                                                                                                                                             >  2:238052/1‑251 (MQ=255)
                                                                                      gACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGAt                                                                                                                                                             >  1:603895/1‑251 (MQ=255)
                                                                                          atccccAGTCCCATACCGACAAAACCAACCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCAGGATGATAGGGACGAAACCGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCg                                                                                                                                                         <  1:336661/247‑1 (MQ=255)
                                                                                            cACCAGTCCCAAACCGACAAAACCAATCCCAACCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCATGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTa                                                                                                                                                       <  1:390660/251‑1 (MQ=255)
                                                                                              ccAGCCCCATACCAACAAAACAAACCCCACCCGTTACCCCTGCCGAAACCGCCGGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAAt                                                                                                                                                     <  1:542842/251‑1 (MQ=255)
                                                                                               cAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATc                                                                                                                                                    >  1:618138/1‑251 (MQ=255)
                                                                                                   cccATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGCTGCCGTAATCACTg                                                                                                                                                >  1:26209/1‑251 (MQ=255)
                                                                                                   cccATACCGACAAAACCAATCCCAACCGTTACCACTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTACCGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTg                                                                                                                                                >  1:199082/1‑251 (MQ=255)
                                                                                                    ccATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGt                                                                                                                                               >  1:247994/1‑251 (MQ=255)
                                                                                                         ccGACAAAACCAACCCCACCCGTTCCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGtt                                                                                                                                          <  1:52444/251‑1 (MQ=255)
                                                                                                          cGACAAAACCAATCCCAACCGTTACCCCTGCCTTAAGCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTc                                                                                                                                         <  1:537516/251‑1 (MQ=255)
                                                                                                             cAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGc                                                                                                                                       >  2:197253/1‑250 (MQ=255)
                                                                                                             cAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGc                                                                                                                                       >  2:192251/1‑250 (MQ=255)
                                                                                                              aaaaCCAATCCCAACCTTAGCCCCGGCCTTAATCGCCGGCAGCCAGGGGATGCGAAATATAAGGCCAATCAGGATCAGTAGAATGGGAACGAAAACGTGACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTTAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCAt                                                                                                                                     <  1:268515/251‑1 (MQ=255)
                                                                                                                aaCCAATCCCAACCGTTACCCCTGACTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATgg                                                                                                                                   >  2:132456/1‑251 (MQ=255)
                                                                                                                 accacccccAACCGTCACCCCTGCCATAATCGCGTGGAGCCAGGGGATGCGGAAAATAAGGCCAATCACGATCATGATGAGAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAATAGAGCCTTGAGTTGTTGTTTTAACGCGTCGTCGTTGAGCCCTTTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTTGCATGGa                                                                                                                                  <  2:253840/245‑1 (MQ=255)
                                                                                                                        ccAACCGGTACCCCTGCCTTAATCGCCTGCAGCCAGGGGAAGCGGAAGATAAGACCAACCACGATCATGATTATAGGGAAGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCTTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGAAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACa                                                                                                                           <  1:394156/251‑1 (MQ=255)
                                                                                                                           aCCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAt                                                                                                                         >  2:16326/1‑250 (MQ=255)
                                                                                                                             cgtaaCCCCTGCCTTAATCGCCTGCAGCAAGGGGATGCGGAAGATAAGGCCAAACACGATAATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAg                                                                                                                      <  1:117907/247‑1 (MQ=255)
                                                                                                                                 aCCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTc                                                                                                                   <  1:516541/250‑1 (MQ=255)
                                                                                                                                   ccccGGCCTTAATCGCCTGCAGCCATGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTTCTGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATGTTATTGAGTTAACAGAGCCGTGATTTGTTGTTTGAACGCGTCGTCGTTTATCCCTGTGAGCAGGGCGGCTCCAGTAATCGTAGTGATGCCGTAATCACAGTTGTTTCGCAGGGATGTTGCAATAAGGACCa                                                                                                                 <  2:40766/251‑1 (MQ=255)
                                                                                                                                          cTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGtt                                                                                                          >  2:8424/1‑251 (MQ=255)
                                                                                                                                                      cAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGCGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTCACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAgg                                                                                              >  1:392740/1‑251 (MQ=255)
                                                                                                                                                      cAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCATCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAgg                                                                                              >  2:210782/1‑251 (MQ=255)
                                                                                                                                                         ccAGGGGATGCGGAAGATAAGGCCAATCCCGATAATGAGGAAAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAt                                                                                           <  1:61503/251‑1 (MQ=255)
                                                                                                                                                         ccAGGGGATGCGGAAGATAAGGCAAATCCCGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGGATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAt                                                                                           <  2:247994/251‑1 (MQ=255)
                                                                                                                                                            ggggATGCGGAAGAAAAGACAAATCACGATCATGATGATAGGGCCGAAAACGTTACCGCCGAGATACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGTCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAAGAAGGTCCTTATCGTTACAGTTAATAGGGAtctc                                                                                        <  2:561966/251‑1 (MQ=255)
                                                                                                                                                                  gCGGAAGATAAGCCAAACCACGATAAGGAGGTTAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATTCCTCAGTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACGGTTTGGTCGCATGGAGGTTACAATAAGGTCCATTCCGTTACGGTTAAGAGGGATCTCATTCAg                                                                                  <  2:626367/251‑1 (MQ=255)
                                                                                                                                                                     gAAGATAAGGCCAATCACGACCATGATGATAGGGACGAAAACGTTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACa                                                                               <  1:1336/251‑1 (MQ=255)
                                                                                                                                                                          tAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCAtt                                                                          >  1:283704/1‑251 (MQ=255)
                                                                                                                                                                          tAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCAtt                                                                          >  1:115812/1‑251 (MQ=255)
                                                                                                                                                                           aaGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTg                                                                         >  1:580264/1‑251 (MQ=255)
                                                                                                                                                                            ggaccaatcacaagcagaaggatAGGGAAGAAAAGGTGACCGGCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGGAGGGATGCCGTAATCACTGTTGTTTGGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGGCAGCATTgg                                                                        <  2:207964/231‑1 (MQ=255)
                                                                                                                                                                            atccgcatcccaaTCATGATGATAGGAACGAAATCTGTACCTCCGAGAGACATGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGTGCAGGGCGGCGCCATTAAGCGTAGGTAGTCCGTAATCACTGTTTGTTCGCATGGAGGTTTCAATATGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTgg                                                                        <  1:81890/239‑1 (MQ=255)
                                                                                                                                                                             ggCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTggg                                                                       <  1:129109/251‑1 (MQ=255)
                                                                                                                                                                                  aTCACGATAATGATGATAGGGAAGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTa                                                                  <  1:529575/251‑1 (MQ=255)
                                                                                                                                                                                       gATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGc                                                             >  1:14077/1‑251 (MQ=255)
                                                                                                                                                                                        ataaggatgatAGGGACGAAAAAGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGGTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCt                                                            <  2:614569/246‑1 (MQ=255)
                                                                                                                                                                                           atgatgATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGa                                                          >  1:619348/1‑250 (MQ=255)
                                                                                                                                                                                              atgatAGGGACGCAAACAATACCGCCGAGAGACAGGATGTAATCAAACATAATGACTCCTTTTATTGAGTTAACAGCGCCTTGATTTGTTGTTTTAACGCCTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGAGGCCGTAATCACTGCTGGTTCGCTTGGAGGTTACCATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGACGGTAGTGGCTGAAATg                                                      >  2:296484/1‑251 (MQ=255)
                                                                                                                                                                                                 atAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGccc                                                   <  1:53686/251‑1 (MQ=255)
                                                                                                                                                                                                   aggtaCGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACCTAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTGCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTg                                                 <  1:244233/247‑1 (MQ=255)
                                                                                                                                                                                                     ggACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGtt                                               <  1:8424/251‑1 (MQ=255)
                                                                                                                                                                                                       aataaatcGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtca                                             <  2:4269/244‑1 (MQ=255)
                                                                                                                                                                                                              cGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCATTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAgg                                       >  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)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                    
CCTGCGCCAACGTCAAAGACATGCAGGGCTAGGCCGAAACGCTCAATCATCACTTTGATAGGGGGACTGAGGCTGTCGATCGCCATGACGATCACCAGTCCCATACCGACAAAACCAATCCCAACCGTTACCCCTGCCTTAATCGCCTGCAGCCAGGGGATGCGGAAGATAAGGCCAATCACGATCATGATGATAGGGACGAAAACGGTACCGCCGAGAGACAGGATGTAATCAAACATAATGCCTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCG  >  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: