Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,440,768 G→A L150L (CTG→TTG)  paeA ← putative cadaverine/putrescine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,440,7680GA98.8% 260.2 / ‑4.6 85L150L (CTG→TTG) paeAputative cadaverine/putrescine exporter
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (41/42);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.59e-01

CTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCC  >  NC_000913/4440528‑4441012
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
cTTCGTGGAGCTCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTAAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCa                                                                                                                                                                                                                                             >  2:245757/1‑250 (MQ=255)
cTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCa                                                                                                                                                                                                                                             >  1:516746/1‑250 (MQ=255)
   cGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAAc                                                                                                                                                                                                                                          <  1:287288/250‑1 (MQ=255)
      ggAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTAc                                                                                                                                                                                                                                      <  1:549979/251‑1 (MQ=255)
            cAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGc                                                                                                                                                                                                                                  >  1:372148/1‑249 (MQ=255)
             aaaCCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGc                                                                                                                                                                                                                                  >  1:220425/1‑248 (MQ=255)
              aaCCAAATCACGTTGTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCGCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGaa                                                                                                                                                                                                                              <  1:294765/251‑1 (MQ=255)
                 cAAATCACGTTTTAACGCGGTGTCAGTGAAGACGGAACGGTACGGGATTCAAGCTCAAAGGCGTTTTCAATCAGCTCGTGTTACTGTTTAAGTCACACGCCCGCAAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTATACGTACCATTGGCAGTTTGGAACTACGGAAGGTTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAg                                                                                                                                                                                                                             >  2:601653/1‑249 (MQ=255)
                     ccACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCCCAAAGACGCTTTCAACCAGCTCTTGTTCCTGGGTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATGGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCTTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATc                                                                                                                                                                                                                       <  1:610474/250‑1 (MQ=255)
                      cACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCg                                                                                                                                                                                                                      <  1:360527/251‑1 (MQ=255)
                                   ggTGTCATTGAAGACGGAACGGTACGGGATTCCAGCCCAAAGACTTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACg                                                                                                                                                                                                         <  2:24014/251‑1 (MQ=255)
                                     tgtCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGc                                                                                                                                                                                                        >  1:357619/1‑250 (MQ=255)
                                         aTTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATAAGCTCGTGTTCCTGTTTAAGTAACACGCCAGCCAGCGCACCGGCTTCCACTACAGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCAGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACACAAAGc                                                                                                                                                                                                   >  1:556422/1‑251 (MQ=255)
                                           tgatctgggcaCGGCACGGGAATCCAGCTTGAAGACGTTTTCAATCGGCTCGTGCTCTTGTATACGTAACACGCCGGCCAGAGGACCAGCTTCACCTACGGCGTAGATGTCATCAGCAGCGATGTCATCTTTACGTACTATTAGCAGTTTGAAAATACGAAAGATTATGTTCGCCAGGCCGTTGAAGAATACCACCACCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGAGGAAACGCAAAGcc                                                                                                                                                                                                  >  2:101182/11‑250 (MQ=255)
                                            gAAGACGGAACGGTAAGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGATCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATAGGGTTGATGATACGCAAAGCcac                                                                                                                                                                                                >  2:310689/1‑251 (MQ=255)
                                                                  ccAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCg                                                                                                                                                                          >  1:285904/1‑251 (MQ=255)
                                                                     gctcccagccGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCCGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACCGGCAGACGCGCATCGGGTTGAAGATACGCAAAGCCACAGCTTCTGGCGCAATCATCCCgttg                                                                                                                                                                       >  2:54866/10‑248 (MQ=255)
                                                                      cTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGc                                                                                                                                                                      >  2:617774/1‑251 (MQ=255)
                                                                         aaaGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACCCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGCAGCGCATCGGGTTGATGATACGCAACGCAACAGCTTCTGGCGCAATCATACCGATGCGtt                                                                                                                                                                   >  2:296104/1‑251 (MQ=255)
                                                                           aGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGCTGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAACCGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTc                                                                                                                                                                 >  1:188953/1‑251 (MQ=255)
                                                                                 tttCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGt                                                                                                                                                            <  2:197483/250‑1 (MQ=255)
                                                                                 tttCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACAATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGtt                                                                                                                                                           >  1:384825/1‑251 (MQ=255)
                                                                                     aaTCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGAGACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAAt                                                                                                                                                       >  2:447166/1‑251 (MQ=255)
                                                                                      accagtgcgtctccttGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAAAAGTGATGTCATCTTTACGTACCATTGGCAGGGTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTTCAAACGGACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACTGCTTCTGGCTCAATCATCCCGATGCGGTTCGGGTTGAAATc                                                                                                                                                      <  2:556422/236‑1 (MQ=255)
                                                                                             cGTGTTCCTTTTTACGTAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAAc                                                                                                                                               <  2:188891/251‑1 (MQ=255)
                                                                                             cGTGTTCCTGTTTACGTAACCCGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAAc                                                                                                                                               <  1:501178/251‑1 (MQ=255)
                                                                                               tgtTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAAc                                                                                                                                               <  2:2898/249‑1 (MQ=255)
                                                                                                  tCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGAt                                                                                                                                          >  1:484288/1‑251 (MQ=255)
                                                                                                      gTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAAc                                                                                                                                      >  1:103000/1‑251 (MQ=255)
                                                                                                          aCGTAACCCCCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGAGAAACATg                                                                                                                                   <  2:368587/250‑1 (MQ=255)
                                                                                                             aaacccgcccgccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTTGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCAGCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGcc                                                                                                                                 <  2:47046/244‑1 (MQ=255)
                                                                                                                cacGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGTTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCa                                                                                                                            <  2:613324/251‑1 (MQ=255)
                                                                                                                      cgccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACg                                                                                                                      <  1:220235/251‑1 (MQ=255)
                                                                                                                      cgccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCGTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTGCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAAct                                                                                                                      <  2:481021/251‑2 (MQ=255)
                                                                                                                       gccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACga                                                                                                                     <  2:145161/251‑1 (MQ=255)
                                                                                                                        ccAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACg                                                                                                                      >  2:42494/1‑249 (MQ=255)
                                                                                                                          aGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAGGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGaa                                                                                                                  <  1:473573/251‑1 (MQ=255)
                                                                                                                             gcACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGaa                                                                                                                  >  1:462849/1‑248 (MQ=255)
                                                                                                                                ccGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAagagag                                                                                                            <  1:223142/251‑1 (MQ=255)
                                                                                                                                  agCTTCAACTACCGAGTAGATGTCATCAGAAGTGATGTCATCCTTACGTACCATAGGCAGTTTAAAAGTACGGAAGTTTCTGTTCGCCAGGCCGATGACGATCCACACCCACGGGGTGCAAACGTACAGGCAGGAGCGCATCGGGCTGATGATCCGCCAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCGGCCAACAGGATAAACATGCCCGGTACTAACGGAAAAgagaga                                                                                                           >  2:99837/2‑250 (MQ=255)
                                                                                                                                     ttCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGCGAGAATa                                                                                                        >  2:74215/1‑250 (MQ=255)
                                                                                                                                       ccACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATaaa                                                                                                      <  2:479708/250‑1 (MQ=255)
                                                                                                                                        cACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCGCAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGc                                                                                                    >  1:557800/1‑251 (MQ=255)
                                                                                                                                               gcgTAGATGTCATCAGAAGGGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGACGCGCATCGGGTTGATGATACGCACAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTAACTACCGAGAAAGAGAGAATAAAGCTCAGTTg                                                                                             >  2:497250/1‑251 (MQ=255)
                                                                                                                                                 gTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGcc                                                                                           >  2:461416/1‑250 (MQ=255)
                                                                                                                                                 gTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCt                                                                                           >  1:346980/1‑251 (MQ=255)
                                                                                                                                                   aGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGACGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATAAGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAATAGAGAATAAAGCTCAGTTGCTc                                                                                          >  1:598240/1‑250 (MQ=255)
                                                                                                                                                   aGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCa                                                                                         <  2:174440/251‑1 (MQ=255)
                                                                                                                                                    gATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCag                                                                                        <  1:310689/251‑1 (MQ=255)
                                                                                                                                                       gTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCagaga                                                                                     <  2:516746/251‑1 (MQ=255)
                                                                                                                                                         catcaGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCagagag                                                                                    <  1:74215/250‑1 (MQ=255)
                                                                                                                                                           tcaGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGctc                                                                                 >  2:218626/1‑251 (MQ=255)
                                                                                                                                                                   gATGTCATCGTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACa                                                                         <  2:357619/251‑1 (MQ=255)
                                                                                                                                                                        cATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACatat                                                                      >  2:13990/1‑249 (MQ=255)
                                                                                                                                                                        cATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAg                                                                    <  2:7410/251‑1 (MQ=255)
                                                                                                                                                                            tttACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCggg                                                                >  1:9303/1‑251 (MQ=255)
                                                                                                                                                                                cGTACCATTGGCAGCTTGCAACTACGGAACATTCTGTTCGCCAGGCCGTTGAAGAACCACACCAAAGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGCTGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCCGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGgaga                                                             >  2:294259/1‑250 (MQ=255)
                                                                                                                                                                                   aCCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACa                                                          <  2:188953/250‑1 (MQ=255)
                                                                                                                                                                                      aTTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCt                                                      >  1:449327/1‑251 (MQ=255)
                                                                                                                                                                                         ggCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCtgtg                                                   <  1:245757/251‑1 (MQ=255)
                                                                                                                                                                                             gTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGaaaa                                               <  1:536054/251‑1 (MQ=255)
                                                                                                                                                                                                 gcaaaTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAACCATGCCAGTAACTAACGAGAAAGAGAGAATATAGCTCAGTTGCTCAGAGACCTCAGCCGCCATAGAGCGGGAGAAAAGGCTGTGAAAAGCTg                                           >  2:501493/3‑251 (MQ=255)
                                                                                                                                                                                                  aaaaTACGGTAGGTTATGTTCGCAAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACTTCCAGGCCGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGCATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGTAGAACAGGCTGTGAAAAGCTgg                                          >  2:620895/1‑251 (MQ=255)
                                                                                                                                                                                                  aaaaTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTgg                                          <  1:8114/251‑1 (MQ=255)
                                                                                                                                                                                                  aaaaTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGAGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCCAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGCAAAGCTgg                                          >  2:21791/1‑251 (MQ=255)
                                                                                                                                                                                                     aTACGGAAGATTATGTTCGCCAGGCCGTTGCAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATCAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTgg                                          >  2:580397/1‑248 (MQ=255)
                                                                                                                                                                                                     aTACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGgag                                        <  1:580397/250‑1 (MQ=255)
                                                                                                                                                                                                       aCGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGaaa                                     <  1:218626/251‑1 (MQ=255)
                                                                                                                                                                                                            aGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGcc                                <  1:498168/251‑1 (MQ=255)
                                                                                                                                                                                                              aTTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCa                                                                                                                                >  1:543469/1‑153 (MQ=255)
                                                                                                                                                                                                              aTTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCa                                                                                                                                <  2:543469/153‑1 (MQ=255)
                                                                                                                                                                                                              aTTATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGc                              >  1:100820/1‑251 (MQ=255)
                                                                                                                                                                                                               ttATGTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGcc                             >  2:24426/1‑250 (MQ=255)
                                                                                                                                                                                                                   gtccGCCAGGACGTCGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCTGCGAGTTGATGATACGCTAAACAACAGCTTCTGCCGCAATCTTACCGATGCGTTTCGCGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTTACGAGAAATAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCGGGAGAAAATACCGCATCa                          <  2:129469/247‑1 (MQ=255)
                                                                                                                                                                                                                   gTTCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCAc                         <  1:21791/251‑1 (MQ=255)
                                                                                                                                                                                                                   gTTCGCCAGGCCGTTGAAGAACCACAACAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAACGCTGGCGAAAATGCCGCATCAc                         >  2:577524/1‑251 (MQ=255)
                                                                                                                                                                                                                     tCGCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCAcc                        <  1:347452/250‑1 (MQ=255)
                                                                                                                                                                                                                       gCCAGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGGTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACcgac                     <  2:372148/251‑1 (MQ=255)
                                                                                                                                                                                                                       cccAGGCCGTTGAAGAACCGCACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGGTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATAGAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACcgac                     <  2:128083/250‑1 (MQ=255)
                                                                                                                                                                                                                          aGGCCGTTGAAGAACCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAAGGCCGCATCACCGACGAt                  <  2:220425/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gaagaaCCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATAccgccg           <  2:462849/250‑1 (MQ=255)
                                                                                                                                                                                                                                   aagtaCCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCAGCGGGTTGATGATCCGCAAAGCCAAAGCCTCTGGCGAAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCAAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTGGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCGGTGAAAAGCTGGAGAAAATGCCGCATCACCGCCGATACCGCCgag         <  2:24341/247‑1 (MQ=255)
                                                                                                                                                                                                                                   aagaaCCACACCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGAAATCATCCCGATGCGCTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACCTATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATAccgccg           >  1:58893/1‑249 (MQ=255)
                                                                                                                                                                                                                                     gaaCCACACCAACGGGGTGCAAACGTCCAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGaa       <  1:209946/251‑1 (MQ=255)
                                                                                                                                                                                                                                          acacCAACGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGACAATGCCGCATCACCGACGATACCGCCGAGAATCGcc  >  2:489017/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    
CTTCGTGGAGATCAAACCAAATCACGTTTTCACGCGGTGTCATTGAAGACGGAACGGTACGGGATTCCAGCTCAAAGACGTTTTCAATCAGCTCGTGTTCCTGTTTACGTAACACGCCCGCCAGCGCACCGGCTTCCACTACCGCGTAGATGTCATCAGAAGTGATGTCATCTTTACGTACCATTGGCAGTTTGAAAATACGGAAGATTATGTTCGCCAGGCCGTTGAAGAACCACACCAGCGGGGTGCAAACGTACAGGCAGAAGCGCATCGGGTTGATGATACGCAAAGCCACAGCTTCTGGCGCAATCATACCGATGCGTTTCGGGGTTAAATCCGCAAACAGGATAAACATGCCAGTCACTAACGAGAAAGAGAGAATAAAGCTCAGTTGCTCAGAGAGCTCTGCCGACATATAGCGGGAGAACAGGCTGTGAAAAGCTGGAGAAAATGCCGCATCACCGACGATACCGCCGAGAATCGCC  >  NC_000913/4440528‑4441012

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: