Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 3,944,733 T→C I75T (ATT→ACT)  mdtM → multidrug resistance protein MdtM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,944,7330TC98.3% 1942.3 / 8.8 525I75T (ATT→ACT) mdtMmultidrug resistance protein MdtM
Reads supporting (aligned to +/- strand):  ref base T (2/7);  new base C (255/257);  total (258/267)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 526 total aligned reads displayed.

CGTATCTGTCGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCGTTGGTTATGT  >  DH1/3944583‑3944880
                                                                                                                                                      |                                                                                                                                                   
  tATCTGTCGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTg                                                                                                                                                   >  4:1666032/1‑151 (MQ=255)
   aTCTGTCGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTgg                                                                                                                                                  >  2:755970/1‑151 (MQ=255)
   aTCTGTCGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTg                                                                                                                                                   <  7:1330593/150‑1 (MQ=255)
        tCGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGc                                                                                                                                              >  5:253773/1‑150 (MQ=255)
         cGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGc                                                                                                                                              >  2:499586/1‑149 (MQ=255)
         cGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGc                                                                                                                                              >  2:1456482/1‑149 (MQ=255)
         cGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAg                                                                                                                                            <  4:1310600/151‑1 (MQ=255)
         cGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAg                                                                                                                                            <  3:1833916/151‑1 (MQ=255)
          gACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAgg                                                                                                                                           <  7:1070038/151‑1 (MQ=255)
          gACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAgg                                                                                                                                           >  7:255980/1‑151 (MQ=255)
           aCGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGc                                                                                                                                          >  6:1633423/1‑151 (MQ=255)
           aCGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGc                                                                                                                                          >  7:667386/1‑151 (MQ=255)
           aCGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGc                                                                                                                                          >  2:1022648/1‑151 (MQ=255)
             ggATCTGATGCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCg                                                                                                                                        <  8:215470/151‑1 (MQ=255)
             ggATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCg                                                                                                                                        >  1:1933044/1‑151 (MQ=255)
             ggATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTCGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGACCTCGCCGCAGGCCg                                                                                                                                        >  4:1203123/1‑151 (MQ=255)
              gATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCgg                                                                                                                                       >  8:546403/1‑151 (MQ=255)
                   gATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGt                                                                                                                                      <  7:910273/147‑1 (MQ=255)
                       cAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGGCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTa                                                                                                                              <  1:1777034/151‑1 (MQ=255)
                             gggATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCgggg                                                                                                                        <  6:1191359/151‑1 (MQ=255)
                              ggATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGgc                                                                                                                       <  4:692013/151‑1 (MQ=255)
                               gATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGgcg                                                                                                                      >  2:378358/1‑151 (MQ=255)
                                aTCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCATTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCCGGgcgc                                                                                                                     >  1:306807/1‑151 (MQ=255)
                                 tCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCt                                                                                                                    >  6:404557/1‑151 (MQ=255)
                                      aaTGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAAttt                                                                                                               <  6:1804711/151‑1 (MQ=255)
                                      aaTGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAAtt                                                                                                                <  1:1021078/150‑1 (MQ=255)
                                         gtgGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTa                                                                                                            <  1:1585759/151‑1 (MQ=255)
                                          tgGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTAc                                                                                                           >  3:1505926/1‑151 (MQ=255)
                                                gTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGc                                                                                                     <  1:1257549/151‑1 (MQ=255)
                                                 tGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGcc                                                                                                    >  8:406860/1‑151 (MQ=255)
                                                  gATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCt                                                                                                   <  3:1774760/151‑1 (MQ=255)
                                                   aTTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTg                                                                                                  >  1:120788/1‑151 (MQ=255)
                                                      ttAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCTCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTgcg                                                                                                <  7:203923/150‑1 (MQ=255)
                                                      ttAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTgcg                                                                                                <  6:738656/150‑1 (MQ=255)
                                                       tAATGCCGACGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGcgcc                                                                                              >  8:1049163/1‑151 (MQ=255)
                                                        aaTGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGcgccg                                                                                             >  1:39583/1‑151 (MQ=255)
                                                        aaTGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGcgccg                                                                                             <  2:1738237/151‑1 (MQ=255)
                                                        aaTGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGcgcc                                                                                              <  4:366865/150‑1 (MQ=255)
                                                        aaTGCCGATGTCAGGCTGGCCCCTGCTGCCGTCAGTCTCTCTCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGcgcc                                                                                              <  8:754898/150‑1 (MQ=255)
                                                           gCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTATACCCTTGCCTGCGCcgcg                                                                                           <  6:602990/150‑1 (MQ=255)
                                                              gATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACa                                                                                        >  6:224180/1‑150 (MQ=255)
                                                               aTGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAAt                                                                                      <  7:1631768/151‑1 (MQ=255)
                                                               aTGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAAt                                                                                      >  6:989952/1‑151 (MQ=255)
                                                               aTGTCAGGCTGGCCGCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCTCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAAt                                                                                      <  7:1811790/151‑1 (MQ=255)
                                                               aTGTCAGGCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAAt                                                                                      <  5:370588/151‑1 (MQ=255)
                                                                ggTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATg                                                                                     <  6:1713799/150‑1 (MQ=255)
                                                                 gTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGt                                                                                    >  8:1103443/1‑151 (MQ=255)
                                                                   cAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGtt                                                                                   <  2:1211966/150‑1 (MQ=255)
                                                                   cAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTc                                                                                  >  8:1193240/1‑151 (MQ=255)
                                                                   cAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCGGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGACCTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTc                                                                                  <  3:1426439/151‑1 (MQ=255)
                                                                   cAGGCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTc                                                                                  <  5:623326/151‑1 (MQ=255)
                                                                    aGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTca                                                                                 <  7:1718775/151‑1 (MQ=255)
                                                                      tCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTcac                                                                                >  2:1806954/1‑150 (MQ=255)
                                                                         ggCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACaa                                                                              <  5:1145317/149‑1 (MQ=255)
                                                                           ccccTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTc                                                                          <  8:1815468/151‑1 (MQ=255)
                                                                            cccTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCt                                                                         <  4:777995/151‑1 (MQ=255)
                                                                              ctgctgCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCCCACCGTCTAt                                                                       >  2:1932358/1‑151 (MQ=255)
                                                                              ctgctgCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTAt                                                                       >  5:1314664/1‑151 (MQ=255)
                                                                               tgctgcCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATg                                                                      >  8:1338461/1‑151 (MQ=255)
                                                                               tgctgcCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCAGCGACAATGTTCACAACGTCTATg                                                                      >  3:1171824/1‑151 (MQ=255)
                                                                                gctgcCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATg                                                                      >  3:601010/1‑150 (MQ=255)
                                                                                  tgcCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGac                                                                    >  5:1393676/1‑150 (MQ=255)
                                                                                    ccttCAGTCTCTTTCTTGCTCGCTGTATGGCGGTACAGTGGCTTCTGGGGCCGCTTTCCGACATAACTGGCCGCAGGCCGGTGCTTATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGAcaca                                                                 <  2:889341/148‑1 (MQ=255)
                                                                                     cGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGacac                                                                  >  2:991508/1‑149 (MQ=255)
                                                                                        cAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGttt                                                             >  8:807768/1‑151 (MQ=255)
                                                                                            ctctATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCtt                                                          >  4:300242/1‑150 (MQ=255)
                                                                                             tctATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTAt                                                        >  7:435305/1‑151 (MQ=255)
                                                                                                aTCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATcgc                                                     >  3:40521/1‑151 (MQ=255)
                                                                                                aTCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATcg                                                      <  3:1242643/150‑1 (MQ=255)
                                                                                                  cTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATcgcgc                                                   >  8:1416070/1‑151 (MQ=255)
                                                                                                     gCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGt                                                 <  8:667919/150‑1 (MQ=255)
                                                                                                        ggCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGc                                               <  4:1547735/149‑1 (MQ=255)
                                                                                                          gggTTTGGCGTTACAGTGGCTGCTGTTGCCGCGTTCGGACAGAACTGGCCTCAGGCCGGTGCTGATTACAGTGGCGCTAATTTTTACCCTTGCCTGCTCCGCGACAATGTTCACAACGTCTATGACAGAGTTTCTTATCGCGCGTGCAAtt                                           <  6:1409978/150‑1 (MQ=255)
                                                                                                           ggTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTc                                          <  5:988453/151‑1 (MQ=255)
                                                                                                           ggTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTc                                          >  5:1738963/1‑151 (MQ=255)
                                                                                                            gTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCa                                         <  4:978776/151‑1 (MQ=255)
                                                                                                              aTGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAg                                        <  2:1525269/150‑1 (MQ=255)
                                                                                                                ggCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGc                                     <  4:1448462/151‑1 (MQ=255)
                                                                                                                   gTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCAcc                                  >  2:1621668/1‑151 (MQ=255)
                                                                                                                     tACAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCAcc                                  >  4:559690/1‑149 (MQ=255)
                                                                                                                       cAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTa                              <  5:1554746/151‑1 (MQ=255)
                                                                                                                       cAGTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGt                               >  2:328240/1‑150 (MQ=255)
                                                                                                                         gTGGCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATc                            <  1:1470071/151‑1 (MQ=255)
                                                                                                                           ggCTGCTGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCt                           <  2:601388/150‑1 (MQ=255)
                                                                                                                            gctgctGGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTAAAGGGGAGCTAATTTTTAACCTTGCCTGCGCAGCCCCAATGTTAACAAAGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGt                         >  1:696603/1‑151 (MQ=255)
                                                                                                                             ctgctgGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGtt                        >  3:1647305/1‑151 (MQ=255)
                                                                                                                             ctgctgGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGtt                        >  1:969421/1‑151 (MQ=255)
                                                                                                                              tgctgGGGCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGTGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGttt                       >  2:269405/1‑151 (MQ=255)
                                                                                                                                   gggCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTg                  >  7:1481756/1‑151 (MQ=255)
                                                                                                                                   gggCCGCTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTg                  >  1:1662923/1‑151 (MQ=255)
                                                                                                                                      ccGCTGTCCGACAGAACGGGCCGCAGGGCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCa               <  8:1643549/151‑1 (MQ=255)
                                                                                                                                         cTTTCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCg            >  5:1891116/1‑151 (MQ=255)
                                                                                                                                           ttCCGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCGtt          >  8:1881499/1‑151 (MQ=255)
                                                                                                                                              cGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACAGAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCGTTGGt       <  5:78879/151‑1 (MQ=255)
                                                                                                                                              cGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCGTTGGt       <  4:1925080/151‑1 (MQ=255)
                                                                                                                                              cGACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAAGTCAGGGCACCAGTATCTGTTTTATTGCCACCGTTGGt       >  4:107336/1‑151 (MQ=255)
                                                                                                                                               gACAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCGTTGGtt      >  3:832860/1‑151 (MQ=255)
                                                                                                                                                 cAGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCGTTGGTTAt    <  1:496970/151‑1 (MQ=255)
                                                                                                                                                  aGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCGTTGGTTATg   >  1:842799/1‑151 (MQ=255)
                                                                                                                                                  aGAACTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGGATCTGTTTTATTGCCACCGTTGGTTATg   >  1:1390822/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CGTATCTGTCGACGGATCTGATCCAGCCTGGGATCATTAATGTGGTACGTGATTTTAATGCCGATGTCAGTCTGGCCCCTGCTGCCGTCAGTCTCTATCTTGCTGGCGGTATGGCGTTACAGTGGCTGCTGGGGCCGCTTTCCGACAGAATTGGCCGCAGGCCGGTGCTGATTACCGGGGCGCTAATTTTTACCCTTGCCTGCGCCGCGACAATGTTCACAACGTCTATGACACAGTTTCTTATCGCGCGTGCAATTCAGGGCACCAGTATCTGTTTTATTGCCACCGTTGGTTATGT  >  DH1/3944583‑3944880

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: