Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,854,297 G→A E126K (GAA→AAA)  emrD → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,854,2970GA99.4% 552.1 / NA 163E126K (GAA→AAA) emrDmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (79/82);  total (79/84)
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.98e-01

Only 100 of 167 total aligned reads displayed.

GGGCGGTGCAGAGCGTAATGGGCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATGAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGCTTTTCGGTAACAGCGG  >  NC_000913/3854049‑3854544
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
ggtcGGTGCAGAGCGTAATGGGCGCTTTTCTGCTTACTTACGGTGTCTCACAGCTGTTTTCTGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGGCGGAATGTCCATTTGTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATaaa                                                                                                                                                                                                                                                        <  2:435343/248‑1 (MQ=255)
  gcgtgtcaggtcctagtGGGCACTTATCTGCTGCCTTACGGTGTCCCACAGCTGTTTTATGGCCCGTTTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATTTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACg                                                                                                                                                                                                                                                      <  2:823185/235‑1 (MQ=255)
  gCGGTGCGGGGCTTAATGGCTCCTTATCTGCTGCCTTACGGTTTCTCACAGCTGTTTCATGGCCCGATGTCCGCCCGCGTGGGCCGCCGACCGGTTATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGCCCTCCAGGTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACg                                                                                                                                                                                                                                                      <  2:1446172/251‑1 (MQ=255)
    gttGCAGAGCGTAATGGCCTCATATCTGCTGTCTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGa                                                                                                                                                                                                                                                    <  1:986518/249‑1 (MQ=255)
      tGCAGAGCGTAATGGGCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGCCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGCAATGGACGGGACTGTACCGTGGGGTTTATACAAACGGa‑cc                                                                                                                                                                                                                                                 >  2:353339/1‑250 (MQ=255)
        cAGAGCGTAATGGGCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGGGATTTATATAAACGGA‑caca                                                                                                                                                                                                                                               >  2:668077/1‑251 (MQ=255)
        cAGAGCGTAATGGGCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGCTTTATATAAACGGA‑caca                                                                                                                                                                                                                                               >  1:1159366/1‑251 (MQ=255)
             cGTAATGGGCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTgc                                                                                                                                                                                                                                          >  1:625773/1‑251 (MQ=255)
              gtatttggtccTTACCTGCTTACTTACGGTTTCTCTCCGCTGTTTTTTGTCCCTATTTCCGACCTCGTGGGCCGCCGTCCGGTGATCCTCCTCGGGATGTCCATTTTTATTCTGGCAACGCTGGTCGCGCTCACGGCCTCCAGTTTGACGTTGTTTATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTgcg                                                                                                                                                                                                                                         <  2:1536638/241‑1 (MQ=255)
              gTAATGGGCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGAACTTCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGAGGGGGACCGGCGGTGGGGGCGTAATGGCGCGGTGTTTACCGGGAGATTTAATTAAACgggacaggttgcg                                                                                                                                                                                                                                          >  2:1582266/1‑240 (MQ=255)
                  tGGGCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGCTTTACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCCACGGCGTTGGCGGCGCAATGGCTCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGCTGCGCCAt                                                                                                                                                                                                                                     >  2:1436306/1‑251 (MQ=255)
                   gggCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGCATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGAGTAATGGCGCGTACTTTACCGCGTGAGTTATATAATCGGA‑CACAGTTGCGCCATg                                                                                                                                                                                                                                    >  1:284226/1‑251 (MQ=255)
                    tgtataagggaccggtcTTACGTTGTCTCACAGCTGTTTTATGCCCCGATTTCCGACCGCTTGGGCCGCCGACCGGTGATCCTCGTCGGATTGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGc                                                                                                                                                                                                                                   <  2:1164457/235‑1 (MQ=255)
                             ctgctgACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGGCCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAAAAGc                                                                                                                                                                                                                            >  2:1141590/1‑249 (MQ=255)
                             ctgctgACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCCCCGACCGGTGATCCTCGTCGGCATGTCCATTTTTATGCTGGCAACGCCGCGCGCGGTTACGACCTCAACTTCGACGGTGTTGATTCCCGCCAGAGCGGCGCAGGGGACGGGCACTGGCCTTGGCGGCGTATTGGCGCGTACTTTACCGCGAGCATTATAGAAACGGACCACA‑TTGCGCCCTGCCCACACCCt                                                                                                                                                                                                                          >  2:1078744/1‑251 (MQ=255)
                                ctgACTTACGGTGTCTCACATCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGtt                                                                                                                                                                                                                       <  1:323430/251‑1 (MQ=255)
                                             tctcACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATTGGGCtt                                                                                                                                                                                                          >  1:1339305/1‑251 (MQ=255)
                                                cacaGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATtctc                                                                                                                                                                                                       >  1:417766/1‑251 (MQ=255)
                                                      gtttttATGGCCCGGTTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGt                                                                                                                                                                                                 <  2:79669/249‑1 (MQ=255)
                                                         tttGTGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATTTCCATTTTTCTGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACAt                                                                                                                                                                                                                 <  1:154418/232‑1 (MQ=255)
                                                          ttttGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGt                                                                                                                                                                                             <  2:494223/248‑1 (MQ=255)
                                                          ttATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCATGTTAAAAATGGGGAt                                                                                                                                                                                                           >  1:1164223/1‑237 (MQ=255)
                                                                cccGTTTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCg                                                                                                                                                                                       <  2:829534/251‑1 (MQ=255)
                                                                  cGATTTCCGCCCGCGTGCGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGCCCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCa                                                                                                                                                                                     <  2:1380107/251‑1 (MQ=255)
                                                                  cGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAAAGGA‑CACAGTTGGGCCATGCCAACAGCCTGTTTAAAATGGGGATTCTCTGCAGTCCGTTGCTCGCa                                                                                                                                                                                     >  1:1002423/1‑251 (MQ=255)
                                                                    aTTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGc                                                                                                                                                                                                                            <  1:1712102/210‑1 (MQ=255)
                                                                     tttCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACCGt                                                                                                                                                                                                                                             >  2:1727956/1‑192 (MQ=255)
                                                                      ttCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGTCGGCGTCACGGGGCGTACTTTACCGCGACATTGATATAAACGGG‑CACAGTGGCGCCCTGCTAACAGACTGTTAAACATGGGGATTCGAGTCAGTCCGTTGTTCGCCCCGc                                                                                                                                                                                 >  2:152800/1‑251 (MQ=255)
                                                                       tccatccacgTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCTGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCt                                                                                                                                                                                <  1:883475/243‑1 (MQ=255)
                                                                          gACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGCTGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGCTTCCCGTCAGTCCGTTGCTCGCACCGCTAAt                                                                                                                                                                             >  2:1702618/1‑251 (MQ=255)
                                                                              gcgTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTTGTCGCGGTCACGACCTCCCGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATcggc                                                                                                                                                                         <  1:1431759/251‑1 (MQ=255)
                                                                              gcgTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATcggc                                                                                                                                                                         <  1:667931/251‑1 (MQ=255)
                                                                                 tGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGt                                                                                                                                                                      >  1:1003083/1‑251 (MQ=255)
                                                                                           aCCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGCCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATa                                                                                                                                                            >  1:536545/1‑251 (MQ=255)
                                                                                           aCCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGATGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGAGTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTTGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATa                                                                                                                                                            >  1:985245/1‑251 (MQ=255)
                                                                                            ccGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGCTGGCGGCGTAATGGCGCGCACTTTACCGCGAGATTAATATACACGGA‑CACAGTTGCGCCATGCTACCAGCCTGTTAAACATGGGGATTCCCGTc                                                                                                                                                                                                    >  2:1712445/1‑210 (MQ=255)
                                                                                              ggTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACAt                                                                                                                                                                                                                 >  1:827637/1‑195 (MQ=255)
                                                                                               gTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTTCTCGCAcc                                                                                                                                                                                   <  1:644868/224‑1 (MQ=255)
                                                                                                      tcgtctGAATGTCCATTTTTATGCTAGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGACGTAATGGCGCGTACTTTACCGCGAGTTTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCTGCGGTCTGCTGGATACGATGTGGAAc                                                                                                                                                 <  1:1148082/251‑1 (MQ=255)
                                                                                                       cgtcgGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGGGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACt                                                                                                                                                <  2:147599/251‑1 (MQ=255)
                                                                                                             aaTGTCCGTTTGTGTGCCGGCACCTCTGGTCGCGGTCCCGCCATCCAGTTGGCCGGTGTTGAGTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGgcgcg                                                                                                                                          <  2:983145/251‑1 (MQ=255)
                                                                                                               tGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTGGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCTTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTTGCACCGCTAATCGGCGGCCTGCTGGATACGATGTGGAACGGgcg                                                                                                                                            >  2:416565/1‑247 (MQ=255)
                                                                                                                 tCCATTTTTTTGCTGGAAACGCGTGTCGCGGTCACGACCTCCAGTTTTCCGGTGTTGTTTGCCGCCAGCGCAGTGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTg                                                                                                                                      <  2:1003994/251‑1 (MQ=255)
                                                                                                                  ccATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCAcc                                                                                                                                                                                   <  1:1721014/205‑1 (MQ=255)
                                                                                                                    aTTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGCTa                                                                                                                                   >  1:1508760/1‑251 (MQ=255)
                                                                                                                     tttttATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCg                                                                                                                                                                                                      <  1:1237951/183‑1 (MQ=255)
                                                                                                                     tttttATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCa                                                                                                                                                                                     <  1:188942/200‑1 (MQ=255)
                                                                                                                       tttATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCGGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATAc                                                                                                                                                           >  2:645002/1‑224 (MQ=255)
                                                                                                                              tGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGGGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTttg                                                                                                                         <  2:1509058/251‑1 (MQ=255)
                                                                                                                              tGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGAGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTTGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCGGTTTACTCTTTttg                                                                                                                         >  2:639922/1‑251 (MQ=255)
                                                                                                                               gacgacgagggTCGCGGCCACGACCCCCAAGGTGACGGTATTGATTGCCGCCAGCGCGATGCAAGGGGTGGGCACCGGCGTTGGGGGCATAATGGCGCTTACTTTACCGCGAGAATTATATAAACGGA‑CACAGTTGGGCCATGCTAACAGCCTGTTAAACATGGAGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGTCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTttgt                                                                                                                        <  1:931516/242‑1 (MQ=255)
                                                                                                                                      cTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCt                                                                                                                 <  1:438631/251‑1 (MQ=255)
                                                                                                                                      cTGGTCGCGGTCACGACCTCCAGCTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGCAATGGCGCGTACTTTACCGCGAGATTTATATATACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTCAAACACGGGGATTCTCGTCCGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGGGGAACTGGCGCGCCTGTCATCCCTTTTTGTGGGtgtt                                                                                                                 >  2:931707/1‑248 (MQ=255)
                                                                                                                                        ggTCGCGTTCACGACCTCCATTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCttt                                                                                                               <  1:1581950/251‑1 (MQ=255)
                                                                                                                                          tCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATAc                                                                                                                                                           >  2:1721362/1‑205 (MQ=255)
                                                                                                                                             cgGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGAt                                                                                                                                                             >  2:189004/1‑200 (MQ=255)
                                                                                                                                             cgGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCt                                                                                                          >  1:1238557/1‑251 (MQ=255)
                                                                                                                                                      ccTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGt                                                                                                                                                                                                                                             >  2:883478/1‑111 (MQ=255)
                                                                                                                                                      ccTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGAcc                                                                                                 >  1:1032061/1‑251 (MQ=255)
                                                                                                                                                       cTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATctct                                                                                                                              >  1:1708215/1‑221 (MQ=255)
                                                                                                                                                                    gtgtTGATTGCCGCCAGCGCGATGCAGGGGATGGGCCCCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATtct                                                                                                                                                                                                        <  1:686568/134‑1 (MQ=255)
                                                                                                                                                                       ttGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAAAAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGACGGTACTCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGTCCTTCAGTCTTGCCCGcgt                                                                                >  2:1234157/1‑249 (MQ=255)
                                                                                                                                                                        tGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTtgtg                                                                                                            <  2:605929/222‑1 (MQ=255)
                                                                                                                                                                          aTTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGAt                                                                                                                                                                                                           <  1:415491/125‑1 (MQ=255)
                                                                                                                                                                          aTTCCCGCCAGCGAGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCTTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACTATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGAt                                                                             <  2:1592881/251‑1 (MQ=255)
                                                                                                                                                                            ggccaccAGCGCGATGCAGGGGCTGGGGCCCGGCGTTGGCAGCGTAATGGCGCGTACTTTCCCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTTGTGTGACCTTCAGTATGGCCCGCTGGATGc                                                                           <  1:836404/250‑1 (MQ=255)
                                                                                                                                                                             gccgccAGCGCGCCGCCGGGGCTGGGCACCCGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGCTTTATATAACCGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGCATCCCGGTCCGTCCGTCGCCCGCACCGCTACTCGCCGCTCCGCTGCACACGATTTGGCACCGCCCCGCCAGTCATCCCCTTCTCCTGGTTCTTTGTCCCGGGTTGCCCTCCTATATGGCCCCCTGGCTGcg                                                                          >  2:1189114/1‑250 (MQ=255)
                                                                                                                                                                             gccgccAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTTGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGcc                                                                          >  2:1338050/1‑251 (MQ=255)
                                                                                                                                                                             gccgccAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGACGGATGcc                                                                          >  1:1094648/1‑251 (MQ=255)
                                                                                                                                                                             gccgccAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTCGTGTGACCTTCAGTATGGCCCGCTGGagtcc                                                                          >  2:836590/1‑247 (MQ=255)
                                                                                                                                                                                  cAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCCCCGCTAATCGGCGGTCTGCTGGATACGACGTGGAACTGGCGCGCCTGTTTTCTCTTTTTGTTGGTTCCTTGTGCTGGTTGGACCTTCAGGATGGCCCGCTGGATGCGGGgaa                                                                     >  2:1023745/1‑251 (MQ=255)
                                                                                                                                                                                         aTGCCGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGc                                                                           <  2:973188/238‑1 (MQ=255)
                                                                                                                                                                                            cAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGGTGGATGCCGGAAACGCGTCCGGt                                                           >  1:1592550/1‑251 (MQ=255)
                                                                                                                                                                                            cAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGGTGGCGGAAACGCGTTCggg                                                           >  1:1445001/1‑250 (MQ=255)
                                                                                                                                                                                             aGGGGATGGGCACCGGCGTTGGCGGCGGCATGGCGCGTACTTTAGCGCGAGATTTATATAAACGGA‑CACAGTGGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGTCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTc                                                          <  2:1238811/251‑1 (MQ=255)
                                                                                                                                                                                               gggATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGa                                                        <  1:1023535/251‑1 (MQ=255)
                                                                                                                                                                                                ggATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGccc                                                                                    >  1:605794/1‑222 (MQ=255)
                                                                                                                                                                                                ggATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGAt                                                       <  1:958179/251‑1 (MQ=255)
                                                                                                                                                                                                  aTGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCCCGTCAGTCCGTTGCTCGCAcc                                                                                                                                                                                   >  2:415587/1‑125 (MQ=255)
                                                                                                                                                                                                       cACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCACGCTGTATGCCGGAAACGCGTTCGGTCGATTCAAcgc                                                >  1:1676606/1‑251 (MQ=255)
                                                                                                                                                                                                          cggatttGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTAAATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTGAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCa                                             <  2:536667/246‑1 (MQ=255)
                                                                                                                                                                                                          cGGCGTTGTCGGCGTCATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTTGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCa                                             <  1:353261/251‑1 (MQ=255)
                                                                                                                                                                                                          cGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTTCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCa                                             <  2:1347864/251‑1 (MQ=255)
                                                                                                                                                                                                          cGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCa                                             <  1:831025/251‑1 (MQ=255)
                                                                                                                                                                                                            gCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCAcg                                           >  1:1277826/1‑251 (MQ=255)
                                                                                                                                                                                                             cGTTGGCGGCGTAATGGCCCGTTCTTTACCGCGAGTTTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCAcgc                                          <  2:1083825/251‑1 (MQ=255)
                                                                                                                                                                                                               ttGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGGAACGCGTCCGGTCGGTGCACCGCGAACgcgc                                        >  2:384201/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggcggcGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTTTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACCCCCCt                                      >  2:958377/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggcggcGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCAc                                                   >  1:972990/1‑238 (MQ=255)
                                                                                                                                                                                                                 ggcggcGTAATGGCGCGTACTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCt                                      >  1:1347584/1‑251 (MQ=255)
                                                                                                                                                                                                                      cgTAATGGCGCGTACTTTACTGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGCACTGGCGCGTCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTa                                 >  1:148404/1‑251 (MQ=255)
                                                                                                                                                                                                                      cgTAATGGCGCGTACCTTACCGCGCGAGTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGCCGCGCCTGTTATCTCTTTTTGTTGGCTCTTTGTGCTGGTGTGACCTTCAGCATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGACGCACAGCGCACGCACCTGCTTa                                 >  1:1388854/1‑251 (MQ=255)
                                                                                                                                                                                                                        tattggggcttccTTGACCGTGAGAGTCATATGTACCGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGCGTTCTCGTCAGTCCGTTGCTCGCACCCCTACTCGGCGGTCTGCTTGATACGATGTGGAACTGGCTCGCCTGTTTTCTCTTTTTGTTTGTTCTGTGTGCTGGTTTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGATCGATGCACCGCGCACGCGCCTGCTTAcc                               <  2:1098432/239‑1 (MQ=255)
                                                                                                                                                                                                                                    cTTTACCGCGAGATTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACg                   >  2:1269623/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tttACCGCGAGATGTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGGCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGc                  <  2:998382/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tttACCAGGAGATTTATATAATCGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTGTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGc                  <  1:1078516/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gaTTTATATAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGGTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGCTGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGCTTTTCGGTaa        >  1:1401715/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   tatataAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGCTGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCCCCGCGCACGCGCCTGCTTACCAGTTATAAAACGCGTTTCGGTAACAGc    >  1:99722/1‑251 (MQ=255)
                                                                                                                                                                                                                                                    atatAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAAc                                                                                                                                                 <  2:330538/109‑1 (MQ=255)
                                                                                                                                                                                                                                                    atatAAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGCTTTTCGGTAACAGCg   <  1:176505/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        
GGGCGGTGCAGAGCGTAATGGGCGCTTATCTGCTGACTTACGGTGTCTCACAGCTGTTTTATGGCCCGATTTCCGACCGCGTGGGCCGCCGACCGGTGATCCTCGTCGGAATGTCCATTTTTATGCTGGCAACGCTGGTCGCGGTCACGACCTCCAGTTTGACGGTGTTGATTGCCGCCAGCGCGATGCAGGGGATGGGCACCGGCGTTGGCGGCGTAATGGCGCGTACTTTACCGCGAGATTTATATGAACGGA‑CACAGTTGCGCCATGCTAACAGCCTGTTAAACATGGGGATTCTCGTCAGTCCGTTGCTCGCACCGCTAATCGGCGGTCTGCTGGATACGATGTGGAACTGGCGCGCCTGTTATCTCTTTTTGTTGGTTCTTTGTGCTGGTGTGACCTTCAGTATGGCCCGCTGGATGCCGGAAACGCGTCCGGTCGATGCACCGCGCACGCGCCTGCTTACCAGTTATAAAACGCTTTTCGGTAACAGCGG  >  NC_000913/3854049‑3854544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑