Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,104,554 T→C Y141C (TAC→TGC)  cpxA ← sensor histidine kinase CpxA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,104,5540TC98.8% 293.4 / NA 80Y141C (TAC→TGC) cpxAsensor histidine kinase CpxA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (46/33);  minor base A (1/0);  total (47/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCGGCAAGGAATTCCTGTGGCCCCGCTTCCAGTTCCGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATC  >  NC_000913/4104314‑4104798
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
gCGGCAAGGAATTCCTGTGGCCCCGCTTCCAGTTCCGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGTCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGGAATTTAATGAAATCGGATGGAGTACTTCTGGCCGGACGAAGCAGATAAAGTGGGCAATTATCTTc                                                                                                                                                                                                                                            <  2:227001/251‑1 (MQ=255)
  ggCAAGGAATTCCTGTGGCCCCGCTTCCAGTTCCGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGc                                                                                                                                                                                                                                          >  1:149028/1‑251 (MQ=255)
   gaaaGGACTTCCTGTGGCCCCGCTTCCAGTTCCGGGTCCTGGCGTGCGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCAGGTTTTGCCAGACGCCAGGCAAACCACAACAGCAGCGGCGTACTGACCAAAATGGGGACAATCAGCAGTAATAGCGGGCGGGCAAACCGTAAGGTAATGAAATCGGTTTGAGAACTGCTGGCCGGACGAATCAGATAAAGATGGCAATTTTCTTCGcc                                                                                                                                                                                                                                         >  2:213021/3‑251 (MQ=255)
   gCAAGGAATTCCTGTGGCCCCGCTTCCAGTTCCGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGcc                                                                                                                                                                                                                                         >  1:91651/1‑251 (MQ=255)
   gCAAGGAATTCCTGTGGCCCCGCTTCCAGTTCCGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGcc                                                                                                                                                                                                                                         >  2:50828/1‑251 (MQ=255)
              cTGTGGCCCCGCTTCCAGTTCCGGGTGCTGGCGTACGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCCGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAACTCGGATTGAGAACTGCTGGCCGGCCGAATCAGATAAAGTTGGCAATAATCTTCGCCATCACGCACgg                                                                                                                                                                                                                              >  2:70543/1‑251 (MQ=255)
                                ttCCGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGtt                                                                                                                                                                                                            >  2:127799/1‑251 (MQ=255)
                                  ccGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCCGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTcc                                                                                                                                                                                                          >  1:31727/1‑251 (MQ=255)
                                  ccGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTcc                                                                                                                                                                                                          <  1:375944/251‑1 (MQ=255)
                                    gggTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTAGCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCAc                                                                                                                                                                                                        <  1:50828/251‑1 (MQ=255)
                                     ggTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCAcg                                                                                                                                                                                                       <  1:25770/251‑1 (MQ=255)
                                               tAAGTTTCCCGGGGCAACTTCATCGGCAGCGTTGTTCAGCTTACGCGCCGGTGTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGGTAACGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGAGAAAGTTGGCAATTATCTTCGCCAGCACGCACGGAGAACGTACCGACCAGTTCCACGCTGCCATAc                                                                                                                                                                                              <  2:80943/250‑1 (MQ=255)
                                                  gTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACtttt                                                                                                                                                                                          <  1:364687/251‑1 (MQ=255)
                                                    ttCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGCACAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATCCTTTttc                                                                                                                                                                                        >  2:340010/1‑251 (MQ=255)
                                                     tCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGCAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTttct                                                                                                                                                                                       >  2:132821/1‑251 (MQ=255)
                                                          gggCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGTAGGACGGACCGACCAATTCCACGCGGTCATACTTTTTCTTCTg                                                                                                                                                                                   <  1:413029/250‑1 (MQ=255)
                                                          gggCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTg                                                                                                                                                                                   <  1:264097/250‑1 (MQ=255)
                                                          gggCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGc                                                                                                                                                                                  >  1:394384/1‑251 (MQ=255)
                                                            gCAACTTCATCGGCAGCGTTTTTCAGCGTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCgg                                                                                                                                                                                <  2:261943/251‑1 (MQ=255)
                                                               aCTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGACAGACTCCAGGCCAACCACAACAGCAGCGGCGAACTAACCAACATGGCGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCGGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACAGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGgatg                                                                                                                                                                             >  1:462063/1‑251 (MQ=255)
                                                                   cATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATc                                                                                                                                                                          <  1:132821/250‑1 (MQ=255)
                                                                         cAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGtt                                                                                                                                                                    >  1:228334/1‑250 (MQ=255)
                                                                           gCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATc                                                                                                                                                                 <  1:255847/251‑1 (MQ=255)
                                                                            cGTTTTTCAGCTTACGCGCCGGTGTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATCGTGACAATCAGCAGTAAGAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTTGCCGGACGAATCAGATAAAGTTGGCAATTAGCTTCGCCATCACGCACGGAGAACGGACCTACCAGTTCCACGCGTCCATACTGGTTCCTCTGCGTATGATCGGCGTTATCg                                                                                                                                                                <  2:462063/251‑1 (MQ=255)
                                                                               ttttCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTTAGAACTGCTGGCCGGACGAAGCAGATAAAGTTTGCAATTATCGTCGCCATCACGCTCGTAGAACGGACCGACCAGTTCCTCGCGGCCATACTTTTTCTGCTGCGGATGATCGGCGTTATCGGc                                                                                                                                                              <  1:174843/250‑1 (MQ=255)
                                                                               ttttCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGc                                                                                                                                                              <  2:256945/250‑1 (MQ=255)
                                                                                 ttCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTg                                                                                                                                                           >  1:201207/1‑251 (MQ=255)
                                                                                   cAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGAc                                                                                                                                                         >  1:25939/1‑251 (MQ=255)
                                                                                         aCGCGCCGGTTTTGCCAGACTCCATGCCAACCACAACAGCAGCGGCGTACTGACCAACAGGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAAGCAGATAAATGTGGCAATTATCTTCGCCATCACGCACGGAGAACGTACCGACCAGTTCCACGCGGCCATACTTTTTCGTCTGCGGATTATCGGCGTTATCGGCCTGACCAATa                                                                                                                                                    <  2:232363/250‑1 (MQ=255)
                                                                                          cgcgCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTAGCGGCCTGACCAATaaa                                                                                                                                                  <  1:341090/251‑1 (MQ=255)
                                                                                              ccGGTTTTGCCAGACGCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAAGGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTAGCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTAGCGGCCTGACCAATAAAGTTa                                                                                                                                              <  2:407517/251‑1 (MQ=255)
                                                                                                 gTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCACACAGTAAGTTACTGACATCGGATTGAGAACTGCTGGCCGGACGAATCAGATCAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCTGCCATACTTTTTCTTCTGCGGATACTCGGCGTTATCCGCCTTCCCAATAACGTTACgc                                                                                                                                           >  2:132531/1‑250 (MQ=255)
                                                                                                 gTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGAGTTATCGGCCTGACCAATAAAGTTACg                                                                                                                                            >  1:457826/1‑250 (MQ=255)
                                                                                                       ccAGACTCCAGGCCCCCCAAAACAGCAGCGGCGTACTGCCCAACATGGCGACAATCCGCAGTAATAGCGGGCGGTCCAACAGCAAGATAATGAAATCGGACTGAGAACCGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGCCCCTACTTTTTCTTCTGCGGACGATCGGCGTTATCGGCTAGACCAATACAGTTACGAATGAtt                                                                                                                                     >  2:52267/1‑250 (MQ=255)
                                                                                                        cAGACTCCAGGCCAAACACAACAGCAGAGGCGTACTGACCACCATGGTGCCAATCAGCAGTAATAGCGGGCGGTCAACCAGTAAGTTAATGAACTCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATCATCGGCGTTATCGGCCTGACCAATAAAGTTACGacggctct                                                                                                                                    >  2:474860/1‑244 (MQ=255)
                                                                                                           aCTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGc                                                                                                                                  >  2:110926/1‑250 (MQ=255)
                                                                                                              ccAGGCCAACCACAACAGAAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGGTACTGAAATCGGATTGAGAACTGCTGGACGGACGAATCAGATAAAGTTGGCCATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGAGGCCATACTTTTTCTTCTGCGGATGCTCGGCGGTAACGGCCTGACCAATACAGTTaggactggacagcattt                                                                                                                              >  2:136751/1‑235 (MQ=255)
                                                                                                                  gCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCt                                                                                                                          <  1:482501/251‑1 (MQ=255)
                                                                                                                        cacaACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTc                                                                                                                    <  2:394384/251‑1 (MQ=255)
                                                                                                                            acaGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGGTACGAATGATCTGCATTTCGCTGCGTTCAgcg                                                                                                                <  2:228334/251‑1 (MQ=255)
                                                                                                                                    ggcgTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACCGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTa                                                                                                                                                                   >  2:326845/1‑192 (MQ=255)
                                                                                                                                    ggcgTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTc                                                                                                                    >  1:209689/1‑239 (MQ=255)
                                                                                                                                    ggcgTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTAGCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTc                                                                                                                    <  2:209689/239‑1 (MQ=255)
                                                                                                                                    ggcgTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGGTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTa                                                                                                                                                                   <  1:326845/192‑1 (MQ=255)
                                                                                                                                      ggTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCAcgc                                                                                                      <  2:31727/250‑1 (MQ=255)
                                                                                                                                      cgTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTGCAGCGCCGATCAcgc                                                                                                      <  1:70543/251‑1 (MQ=255)
                                                                                                                                         aCTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGc                                                                                                   >  1:94893/1‑251 (MQ=255)
                                                                                                                                         aCTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGc                                                                                                   >  1:82270/1‑251 (MQ=255)
                                                                                                                                            gACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAAGCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCtt                                                                                                <  2:149028/251‑1 (MQ=255)
                                                                                                                                                aaCATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATc                                                                                                          >  1:348579/1‑237 (MQ=255)
                                                                                                                                                aaCATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATc                                                                                                          <  2:348579/237‑1 (MQ=255)
                                                                                                                                                      gTGACAATCAGCAGTCATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGACATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGGTCTGAGCCGATCACGCGGCCTTCggtggttac                                                                                      >  2:381273/1‑248 (MQ=255)
                                                                                                                                                               agcagTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACa                                                                              >  1:377663/1‑250 (MQ=255)
                                                                                                                                                                gcagTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACaaa                                                                            <  2:82270/251‑1 (MQ=255)
                                                                                                                                                                   gTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGc                                                                         >  1:303208/1‑251 (MQ=255)
                                                                                                                                                                      aTAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGACCACTGCTGGCCAGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGACTGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGATGt                                                                      >  1:199128/1‑251 (MQ=255)
                                                                                                                                                                        aGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCAACACGCACGGAGAACGGACCGACCAGTCCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCGGACCAATAAAGTTACGAATGATCTGCATTTCGCGGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTc                                                                     <  2:331231/250‑1 (MQ=255)
                                                                                                                                                                           gggcggTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTgg                                                                 >  1:379658/1‑251 (MQ=255)
                                                                                                                                                                                  cAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATCACGTTCCGAATGATCTGCATTTCGCTGCGTTCAGCGCCGCTCACGCGGCCTTCGGTGGTCACCAATCACAAACGCTGTCCTGGCGGTGcc                                                          >  2:103529/1‑251 (MQ=255)
                                                                                                                                                                                         tAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTAt                                                   >  1:272093/1‑251 (MQ=255)
                                                                                                                                                                                                   aaaTCGGATTGAGCACTGCTGGCCGGACGACTCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTCCGAATGATCTGCCTTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACCAACGCTGTCCTGGCGGTGCCCACTTTTCCATCGcccc                                         >  2:458399/1‑250 (MQ=255)
                                                                                                                                                                                                        ggATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGCTCCACGCGGCCCTACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATCACAAACGCTGTCCCGGCGGTGCCCACTTATCAATCCCCCGGAAc                                     >  1:22420/1‑250 (MQ=255)
                                                                                                                                                                                                         gATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAg                                   <  2:379658/251‑1 (MQ=255)
                                                                                                                                                                                                            tGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGAc                                 <  2:91651/250‑1 (MQ=255)
                                                                                                                                                                                                             gagaACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCAAGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGGGCCCACTCATCAATCGCCCGGAACAGAcgc                               >  2:252564/1‑251 (MQ=255)
                                                                                                                                                                                                              agaACTGCTGGCCGGACGAATCAGATAAAGTTGTCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCCCGCGGCCTTCGGGGGTCGCCAATAACAAACGCTGTCCGGGCGGTGCCCACTTATCAATCGCCCGGAACAGAcgcc                              <  2:134550/251‑1 (MQ=255)
                                                                                                                                                                                                               gaACTGCTGGCCGGACGAATCAGATAAATTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAAGAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGc                                                                                                   <  2:299633/181‑1 (MQ=255)
                                                                                                                                                                                                               gaACTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCCTCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATCAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGc                                                                                                   >  1:299633/1‑181 (MQ=255)
                                                                                                                                                                                                                aaCTGCTGGCCGGAAGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTCTCAATCGCCCGGAACAGAcgccgc                            >  1:213292/1‑251 (MQ=255)
                                                                                                                                                                                                                 aCTGCTGGCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACgccgcc                           >  1:6871/1‑251 (MQ=255)
                                                                                                                                                                                                                      tgGCCGGACGCATCAGATAAAGTTGGACATTATCTTCGCCATCACGCACGGAGAACGGACCGAACAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATACAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGGGGTCACCAATAACAAACGCTGTCCTGGCGGAGCCCACTTATCCATCGCCCGGAACAGACGCCGccacccc                      >  2:192312/1‑249 (MQ=255)
                                                                                                                                                                                                                       ggCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCcaccacc                     >  1:172864/1‑250 (MQ=255)
                                                                                                                                                                                                                       ggCCGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACAGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCGGACCAATACAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGGGCCCCCTTATCAATACCCCGGAACAGACGCCGCccccacc                     >  2:40388/1‑250 (MQ=255)
                                                                                                                                                                                                                         ccGGACGAATCCGATAACGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCATTAAAGTTCCGAATGATCTGAATTTAGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGCCCTGGAGGTGACCACCTATCTATCACCCCGAACCGCCGCCGCCACGACAtc                   >  2:294360/1‑250 (MQ=255)
                                                                                                                                                                                                                         ccGGACGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCCTACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGCCGCCGCcaccacctt                   >  2:181110/1‑248 (MQ=255)
                                                                                                                                                                                                                             aCGAATCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAAt               >  1:186127/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTCAGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTCTCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGtt           >  2:3029/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aGATAAAGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTggg        >  1:278716/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaaGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCgg     <  2:26056/250‑1 (MQ=255)
                                                                                                                                                                                                                                          aGTTGGCAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATAGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATc  >  1:344998/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
GCGGCAAGGAATTCCTGTGGCCCCGCTTCCAGTTCCGGGTGCTGGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATC  >  NC_000913/4104314‑4104798

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: