Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,133,175 G→A T13I (ACC→ATC)  yghR ← putative ATP‑binding protein YghR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,133,1750GA98.8% 290.2 / NA 85T13I (ACC→ATC) yghRputative ATP‑binding protein YghR
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (31/52);  total (33/52)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATTTGGTTTTCATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTC  >  NC_000913/3132928‑3133420
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
agaTTCGGTTTTCCTGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCCGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGCGCAGCCCCAGCTAGAGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTAACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCCGCACCGGAATGTAGTTGGGCTGCGGCGCGATTg                                                                                                                                                                                                                                                      >  2:60154/1‑251 (MQ=255)
agaTGTGGTTTTCATGCTTGTGGTTTTCGGGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCTGAATGTAGTTGGGCTGCGGCGCGATTg                                                                                                                                                                                                                                                      <  2:299723/251‑1 (MQ=255)
       gTTTTCAAGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTAtt                                                                                                                                                                                                                                               <  1:501934/251‑1 (MQ=255)
               gCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagttt                                                                                                                                                                                                                                       <  1:246767/251‑1 (MQ=255)
               gCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtt                                                                                                                                                                                                                                        <  1:255116/250‑1 (MQ=255)
                cTTTGGGTTTTCGCGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtttg                                                                                                                                                                                                                                      >  2:436800/1‑251 (MQ=255)
                      gTTTTCTAGGACTTGGCCGCCAGCCGCCGTTCCAGTCAGACTCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCCCCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACaggttgagtgtg                                                                                                                                                                                                                                <  2:251897/251‑4 (MQ=255)
                      gTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGTCGCCAACCAACGGCAATCGGTTGATTTGGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTGGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCATCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACagtttgagtttg                                                                                                                                                                                                                                <  2:276918/251‑1 (MQ=255)
                        tttCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTTAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGt                                                                                                                                                                                                                               <  1:669970/250‑1 (MQ=255)
                          tCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGCCGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCCCCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGt                                                                                                                                                                                                                            >  2:652215/1‑251 (MQ=255)
                                aTTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATcc                                                                                                                                                                                                                      <  2:384211/251‑1 (MQ=255)
                                aTTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGATCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATcc                                                                                                                                                                                                                      >  1:485646/1‑251 (MQ=255)
                                  ttGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCa                                                                                                                                                                                                                     <  1:721863/250‑1 (MQ=255)
                                    ggCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTc                                                                                                                                                                                                                  <  2:319994/251‑1 (MQ=255)
                                     gccgccTGTCGCCGTGCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCGTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTGGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGGGCCATCACCCCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCg                                                                                                                                                                                                                 <  1:60154/251‑1 (MQ=255)
                                       cgccAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGtt                                                                                                                                                                                                               <  1:593186/251‑1 (MQ=255)
                                       cgacAGTCGCCGTTCGGGCCAGACGGCAACCAACGGCAGTCGTTTGATTTTGCCGCCGTCTTCGCCGGAGAGCAGCCCCAGCTAGCGCCGCTCGGTTTGCCAGTGTGGTTGCAGCGATTTCACCAGGTCGGTGGTCGGGTTTGATGGACCCGGAATATCCACCCTGCCTACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCTAGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGtt                                                                                                                                                                                                               <  1:362156/248‑1 (MQ=255)
                                        gccAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTa                                                                                                                                                                                                              <  1:499897/251‑1 (MQ=255)
                                           aGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAat                                                                                                                                                                                                           >  1:516974/1‑251 (MQ=255)
                                              cGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAatca                                                                                                                                                                                                         <  1:369207/250‑1 (MQ=255)
                                               gCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAatcatc                                                                                                                                                                                                       <  1:547224/251‑1 (MQ=255)
                                                  gTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATcccc                                                                                                                                                                                                    >  2:344293/1‑251 (MQ=255)
                                                     ccagccagACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTAt                                                                                                                                                                                                 >  2:154709/1‑251 (MQ=255)
                                                          cagACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCAtt                                                                                                                                                                                            >  2:8676/1‑251 (MQ=255)
                                                           agACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCAttt                                                                                                                                                                                           >  1:478187/1‑251 (MQ=255)
                                                               gCCAACCAACGGCAAGCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCAATTTGGt                                                                                                                                                                                       <  2:93599/251‑1 (MQ=255)
                                                                 caaccaacGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGCGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTg                                                                                                                                                                                      >  2:304475/1‑250 (MQ=255)
                                                                 caaccaacGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGa                                                                                                                                                                                     <  1:636328/251‑1 (MQ=255)
                                                                   accaacGGCAATCGTTTGATTTTGTCGCCGTCGTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATc                                                                                                                                                                                   <  1:521511/251‑1 (MQ=255)
                                                                      aacGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAgg                                                                                                                                                                                >  1:338300/1‑251 (MQ=255)
                                                                         ggCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGaaa                                                                                                                                                                             <  1:351849/251‑1 (MQ=255)
                                                                                      ttttGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAAt                                                                                                                                                                >  1:170610/1‑251 (MQ=255)
                                                                                          gTCGCCGTCGACGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGGCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAAGCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTGAAATCATCCCCTATCCAGTTTGGTGATCAGGAAAGCGCGGCGTGAATTgcg                                                                                                                                                            <  1:652215/251‑1 (MQ=255)
                                                                                                gTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCACTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGt                                                                                                                                                      >  1:230007/1‑251 (MQ=255)
                                                                                                  cTTCGCCGGAGAGCAGCCCCAGATAGCGCAGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGCATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAGGCGCGGCGTGAATTGCGCGTAGTg                                                                                                                                                     >  2:299147/1‑250 (MQ=255)
                                                                                                        cGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTg                                                                                                                                              <  2:596860/251‑1 (MQ=255)
                                                                                                         ggAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGc                                                                                                                                             <  1:8676/251‑1 (MQ=255)
                                                                                                                        aTAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGCTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGAGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATcggcg                                                                                                                              >  2:726295/1‑251 (MQ=255)
                                                                                                                           gcgcCGCTCGGTTTGCCAGTGTTGTGGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAAGCAGCCCCGGAATGTATTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCCTGCAGAATTTCATcggcgtgg                                                                                                                           <  2:357061/251‑4 (MQ=255)
                                                                                                                               cgcTCGGTTTGGCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTc                                                                                                                       <  2:3702/251‑1 (MQ=255)
                                                                                                                               cgcTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCgggggggc                                                                                                                       >  2:402347/1‑247 (MQ=255)
                                                                                                                               cgcTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTc                                                                                                                       <  1:694514/251‑1 (MQ=255)
                                                                                                                                gcTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGGGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCg                                                                                                                      <  1:66657/251‑1 (MQ=255)
                                                                                                                                     gTTTGCCAGAGCTGTTGCAGCGATTTCACCAGGTCAGTGGTCAGTGTGGCTTTACCGGTGCCATCACACCCGACCACCGCACTCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGCTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGCGTGTCGGCCac                                                                                                                 >  2:304850/1‑250 (MQ=255)
                                                                                                                                          ccAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGccccccccgt                                                                                                            >  2:203002/1‑244 (MQ=255)
                                                                                                                                            aGTGTTGTGGCAGCGAGTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGGGCCACCACACCCGCCCACCGCAAGCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCACCCATTCGTTAAATCATCCCCCATCCATTTTGGTGATCACGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTc                                                                                                         <  2:653383/251‑1 (MQ=255)
                                                                                                                                                     gCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCtt                                                                                                <  1:295517/251‑1 (MQ=255)
                                                                                                                                                         cGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGCCTGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTAGCCCCTcccattcacagatctttttg                                                                                             >  2:680285/1‑233 (MQ=255)
                                                                                                                                                             ttCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGACTGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCATCCCGTTCAAGGCTCttttcggcct                                                                                        >  1:329226/1‑244 (MQ=255)
                                                                                                                                                              tCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCAt                                                                                        <  1:299147/250‑1 (MQ=255)
                                                                                                                                                                  cAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGccc                                                                                                                   <  2:565165/220‑1 (MQ=255)
                                                                                                                                                                  cAGGTCGGTGGTCAGTGCGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGccc                                                                                                                   >  1:565165/1‑220 (MQ=255)
                                                                                                                                                                    ggTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCCAGGATCTTTTTGCCATTGAATGt                                                                                 >  1:529661/1‑251 (MQ=255)
                                                                                                                                                                     gTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTc                                                                                <  1:174947/251‑1 (MQ=255)
                                                                                                                                                                       cGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGAc                                                                                                                                                   <  1:710722/183‑1 (MQ=255)
                                                                                                                                                                       cGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAg                                                                              <  1:178574/251‑1 (MQ=255)
                                                                                                                                                                       cGGTGGTCAGGGTGGATTTACCGGAGCCATCACAACCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGAc                                                                                                                                                   >  2:710722/1‑183 (MQ=255)
                                                                                                                                                                         gtggtCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGtt                                                                            <  1:153477/251‑1 (MQ=255)
                                                                                                                                                                           ggtCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCCTCCAGTTCCAGGATCTTTgtggccttttatgtcagtttg                                                                          >  1:616588/1‑231 (MQ=255)
                                                                                                                                                                            gtCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcg                                                                         <  1:92775/251‑1 (MQ=255)
                                                                                                                                                                             gcAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTgcgg                                                                        <  2:275616/250‑1 (MQ=255)
                                                                                                                                                                              cAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTACAGGTTCAGGAACTTTGCGCCATTGGATGTCAGTTgcggg                                                                       >  1:646583/1‑250 (MQ=255)
                                                                                                                                                                                  gtgGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTTCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAATGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTaa                                                                   <  1:196063/251‑1 (MQ=255)
                                                                                                                                                                                  gtgGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTaa                                                                   <  1:415023/251‑1 (MQ=255)
                                                                                                                                                                                  gtgGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTaa                                                                   <  1:304475/251‑1 (MQ=255)
                                                                                                                                                                                    ggATTTACCGGTGCCATCACACCCGACCCCCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTACTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCTCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAAc                                                                  <  2:690672/250‑1 (MQ=255)
                                                                                                                                                                                       tttACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGaa                                                                                                                                                                              >  1:562135/1‑140 (MQ=255)
                                                                                                                                                                                       tttACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGaa                                                                                                                                                                              <  2:562135/140‑1 (MQ=255)
                                                                                                                                                                                           ccGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTTAAGGATATTTTCGCCCTTGAATGTCAGTTGCGGCGTTACgggggttt                                                          >  1:580862/1‑245 (MQ=255)
                                                                                                                                                                                                  cATCACACCCGACCACCTCAATTAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCACCCCCTATCCAGTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTTGCCCATTCAGTTCAAGGATCTTTTCGCCATTGAATGTCAGTTGTGGCGTAACGGCTATTTTGTCCTg                                                  <  2:364066/251‑1 (MQ=255)
                                                                                                                                                                                                   aTCACACCCGACCACCGCAAGCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTGAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGGCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGt                                                 <  1:696569/251‑1 (MQ=255)
                                                                                                                                                                                                     cacaCCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGGCCCTCCAGGTCAAAGATCTTTTCGCCCTTGAATGTCAGTTGCGGGGTTACGGCGGtttttttcttta                                                >  1:456806/1‑244 (MQ=255)
                                                                                                                                                                                                          ccGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGGTCAAGGATCTTTTCGCCCTTTTATGTCAGGTTCGGCGTAACGGCGATTTTTTTCTGTTACGCtc                                          >  1:177861/1‑250 (MQ=255)
                                                                                                                                                                                                           cGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTCCCGGGACTTTGCGCCATTGAATGTCAGTTGCGGCGGTACGGCGATTTTTTTCTTTAACGcctc                                         >  2:627747/1‑248 (MQ=255)
                                                                                                                                                                                                                  ccGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTtg                                   <  1:680285/251‑1 (MQ=255)
                                                                                                                                                                                                                         cAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCa                            <  2:338300/251‑1 (MQ=255)
                                                                                                                                                                                                                         cAGCCCCGGAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAGGTCCAGGATCTTTTCGCCATTGAATGTTCGTTGCGGCGTAACGGGCGATTTTTCCTGTAAAACTGCCacgtggggggct                            >  2:321042/1‑240 (MQ=255)
                                                                                                                                                                                                                                 gAATGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAAGTCCAGGGTTTTTTCGCCCTTGAATGTGAGTTGCGGCGCAACGGCGGTTTTTTTCTGTTAAGCTGCCCGTtgggggggagggttacg                    >  1:125348/1‑237 (MQ=255)
                                                                                                                                                                                                                                  aaTGTAGTTGGGCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATGCGTTAAAGCATCCCCTAGCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTGCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTAcgc                   <  2:17568/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gCTGCGGCGCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGGTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAGCGGCGATTGTTGCCTGTAACGCTGCCAGTTGGTGGGCAGGTTGACGCGCAAACGCGGt        <  2:485646/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gcggcgCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGTAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGttt      <  1:382486/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gcggcgCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGttt      <  2:616588/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggcgCGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTTCAGCTTAAGGATCTTTTCGCCATTGAATGTCAGTTGCGGCGTGACGGCGGTTTTTTTCTGGAACGGTGCCAGGTGGGGGGCAGGTTGACGAGCACAAGCGGgttgtt   >  1:412326/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ggcgCGATTGTACTATTAACAGTTTGAGTGTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGtt   <  2:580862/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   gcgcGATTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTGTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCGGTAACGCTGCCAGTTGGTGGTCAGGGTTACTCGCAAACGCggggtgttc  <  2:667656/251‑8 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       
AGATTTGGTTTTCATGCTTTGGGTTTTCGAGGATTTGGCCGCCAGTCGCCGTTCCAGCCAGACGCCAACCAACGGCAATCGTTTGATTTTGTCGCCGTCTTCGCCGGAGAGCAGCCCCAGATAGCGCCGCTCGGTTTGCCAGTGTTGTTGCAGCGATTTCACCAGGTCGGTGGTCAGTGTGGATTTACCGGTGCCATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTC  >  NC_000913/3132928‑3133420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: