Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 884,387 A→G 7.8% I239V (ATC→GTC)  mdfA → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913884,3870AG 7.8% 344.4 / 2.2 129I239V (ATC→GTC) mdfAmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (43/70);  new base G (10/2);  total (55/74)
Fisher's exact test for biased strand distribution p-value = 4.03e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 131 total aligned reads displayed.

CGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTG  >  NC_000913/884101‑884674
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                     
                 tGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGCTTa                                                                                                                                                                                                                                                                                    >  2:41718/1‑288 (MQ=255)
                   gcggTGCTTGCGCCGCTACTTGGTCCGCTGGTGGGCGCGGCTGGGATCCATGTGCTGCCCTGGTAGGGGATGTTTGTTTTGGTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTatc                                                                                                                                                                                                                                                                                  <  1:288614/288‑1 (MQ=255)
                   gcgcTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTatc                                                                                                                                                                                                                                                                                  <  2:636606/288‑1 (MQ=255)
                             cTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCg                                                                                                                                                                                                                                                                        <  1:1340301/288‑1 (MQ=255)
                                   tACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGGTTGGTTAGTCTGCCGTTGCTGGCGTGGAGCGCCCCGTCGCCGATTTTCATCCTTACCGGgcagca                                                                                                                                                                                                                                                                 >  1:118171/1‑289 (MQ=255)
                                       tGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGGTGACCTGGGA‑‑GGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           <  2:273339/289‑1 (MQ=255)
                                          tcggcgggtggGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGGGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           <  1:1137804/282‑1 (MQ=255)
                                          tCCGCTGGTGGGGGCGGCGTGGATCCTTGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGACGCATGGGCCGCCATCTCCGTTTTCGGACTGCTTGTCGCCATGCCCATAACACACACGCGTATAGGCGAGGTACTTTCACTTAAATAACTCGGTCCGGACTATAAGCTGATGCTGAAGAAGGGCCGCGTTGTGGGGGGGGTGCTGGCGCTGGGGTTTGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGAATATCATCATTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           <  1:1005332/288‑1 (MQ=255)
                                          tCCGCTGGTGGGCGCGGCTTGGAGCCATGTTCTGCCGTGGGAGGTGATGTTTGTTTTGTTTGACGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           <  2:1352581/288‑1 (MQ=255)
                                          tCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGGCCTGAGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGAAACGAGCCATGCCTGAAACCGCCACACGTAAAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTACCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCAACATTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           <  2:127802/288‑1 (MQ=255)
                                          tCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGTTGGGGTGGGTCGTCCAGTCGCCGATTATCATACTTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           >  2:1183123/1‑288 (MQ=255)
                                          tCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           <  2:81241/288‑1 (MQ=255)
                                          tCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           <  1:347653/288‑1 (MQ=255)
                                          gccGCGGGTGGGCCGGGCGTGGATCCCTGTGCTGCCCTGGGAGGTGAGTTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCATCGGGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGGCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGag                                                                                                                                                                                                                                                           <  1:1104994/287‑1 (MQ=255)
                                           ccGCTGGTGGGGGGGGCGTGGATCCACGTGCTGCCCTGGGATGGGGTGTTTGTTTTGTTTGCCGCAATGGCCGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGagc                                                                                                                                                                                                                                                          <  2:582598/288‑1 (MQ=255)
                                            cGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGCGCCGGGATTCGTTAGTCGGGCGTTGTCGGGGGGGTGCGCCCCGTCGCCGGTTATCATACTTCTCGGCGAGCATTTGagcat                                                                                                                                                                                                                                                        >  1:246460/1‑288 (MQ=255)
                                               tggtggGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGGTTTGTTAGTCTGCCGTTGCTGGCGGGGATCGCCCAGGCGGCGGATATCCTCATTACCGGGGAGGCGTTGAGGAGCTa                                                                                                                                                                                                                                                    >  2:1082154/1‑289 (MQ=255)
                                               tggtggGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCt                                                                                                                                                                                                                                                     <  1:355129/288‑1 (MQ=255)
                                                    ggCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTTGTCTGCCGTTTCTGGCGTGGAGCCCCCAGTCGCCGATTATCCTCCTTACCGGCGAGCAGTTTCGCAGCCTTGAat                                                                                                                                                                                                                                              >  1:119199/1‑289 (MQ=255)
                                                      cgcgGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCt                                                                                                                                                                                                                                                     <  2:1360688/281‑1 (MQ=255)
                                                        cggcgTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATgg                                                                                                                                                                                                                                          >  2:1015779/1‑289 (MQ=255)
                                                                aTCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCGGGGATTCGTTAGTCTGCCGTTGCTTGCGTGGGTCGCCCCGTCGCCGATTTTCATCCTTACCGGCGAGCGGTTTAGCCGCTATTAATATGGCCtgtct                                                                                                                                                                                                                                   >  1:1360688/1‑285 (MQ=255)
                                                                   cATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGGTAGTCTGCCGTTGCCGGCGGGGATCGGCCCAGCGCCGCTTATCATCCTTTCCCGCCGGCCGTTGagacgctccttgtctctcatacaatct                                                                                                                                                                                                                                  >  1:1145040/1‑263 (MQ=255)
                                                                     tgtgCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTTACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGGGGTGTTAGTCTGCAGTTGCTGGCGGGGGTCGCCCAGTCGCCGAGTATCATCATTCACGGCGAGCGGTTGGGTGGCGATGGATAGGGCTtgctggcgggg                                                                                                                                                                                                                             >  1:246266/1‑283 (MQ=255)
                                                                     tgtgCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGGTGGGGTTCGTTAGTCTCCCGGTGCGGGCGTGGGTCGGCCAGTCGCCGAGTATCACCTGTACCGGCGAGCAGTTGAGGCGCCTGTAAACTGGCCTGCTGCAAGTg                                                                                                                                                                                                                             >  1:96435/1‑289 (MQ=255)
                                                                      gtgCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGGGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCGTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGc                                                                                                                                                                                                                            <  1:667560/289‑1 (MQ=255)
                                                                            cccTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGGTTCGTTAGTCTGCCGTTGCTGGCGTGGACTGCCCAGTCGCCGGTTATCATCATTACCGGCGGGCAGTTGAGCTGCTGTGAATATGGCGTGCTGCAAGTGCCttgtt                                                                                                                                                                                                                      >  2:220644/1‑285 (MQ=255)
                                                                                gggAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCgg                                                                                                                                                                                                                  <  2:1176243/289‑1 (MQ=255)
                                                                                       gATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGGTTCGTTAGGCTGCCGTTGCTGGCGTGGGTCGCCCAGGCGCCGATTATTATTCTTACCGGGGGGCAGGTGAGCAGCTCAGACATTGGGCTTCTTCAAGTGCCCATTGTCgggggctta                                                                                                                                                                                                           >  1:551058/1‑284 (MQ=255)
                                                                                            ttgttttgttTGCCGCATTGGCAGCGATCTCCTTTTGCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTg                                                                                                                                                                                                       <  2:119199/288‑1 (MQ=255)
                                                                                              gttttgtttGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCgg                                                                                                                                                                                                    <  1:1015779/289‑1 (MQ=255)
                                                                                              gtttggtttgCCGCATTGGCAGCGATCGCCTTTTTCGGTCTGCAACGGGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCgg                                                                                                                                                                                                    <  1:432806/289‑1 (MQ=255)
                                                                                                tgtgtttgCCGCATTGGCCGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGt                                                                                                                                                                                                  <  1:6149/287‑1 (MQ=255)
                                                                                                 ttgtttgCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGACCGGCCACGCGCCGGTGATCATCATTACCCGCGAGGCGTTGGTGTGCTATGAATATTGCTTGCTGCAAGTGCCTACTTTAGTGGCATTAGTTGGGGGtt                                                                                                                                                                                                 >  1:935425/1‑288 (MQ=255)
                                                                                                              tGGAAGCGATGTCCTTTTTCGGTCGGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAgc                                                                                                                                                                                    <  2:41089/289‑1 (MQ=255)
                                                                                                                      aTCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAACGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGCGCTGGGGTTCGGTAGTTTGCGCTTTGTGGGGGGGAGGGGCCAGTCGGCGATTAGTAACATTTCCGTGGAGACGGTGAGCAGGCGTGGAATCGGGTTTGTTGCAGTGCCTTGGTTGTGGGGGGTAATTTCGGGTAGCTTGCTCTTAGGGCGTCtta                                                                                                                                                                           >  1:1081350/1‑287 (MQ=255)
                                                                                                                             ttttCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTcgcg                                                                                                                                                                    <  2:1161913/289‑1 (MQ=255)
                                                                                                                               ttCGGTCTGCAACGGGCCATGCCGGAAACCGCCACGCGTATAGGCGAGGAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGcgcc                                                                                                                                                                  <  1:1282865/289‑1 (MQ=255)
                                                                                                                                  gacaggCAACGAGCCATGCCTGAAACCGCCACGCGGATGGGCGAGAAACGGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGcgccg                                                                                                                                                                 <  2:752448/282‑1 (MQ=255)
                                                                                                                                    tCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCAc                                                                                                                                                              <  2:433234/288‑1 (MQ=255)
                                                                                                                                     cTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGGGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCCGTTGAGGAGCCATGAATATGGGTTGCTGCAAGTGCCTATTTTTGGGGCGTGTATTGCGGGTAACTTGTTGTTAGCGCGTATGAACTCGcgccgccacg                                                                                                                                                            >  2:109549/1‑285 (MQ=255)
                                                                                                                                      tGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGc                                                                                                                                                                                                                                         <  2:302347/212‑1 (MQ=255)
                                                                                                                                             aGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGa                                                                                                                                                                                                                                                            >  1:583303/1‑188 (MQ=255)
                                                                                                                                             aGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGa                                                                                                                                                                                                                                                            <  2:583303/188‑1 (MQ=255)
                                                                                                                                              gCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGGGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCg                                                                                                                                                    <  1:187453/288‑1 (MQ=255)
                                                                                                                                              gCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCg                                                                                                                                                    <  2:917203/288‑1 (MQ=255)
                                                                                                                                               ccatggcgaaaaCCGCCGGGGGTATAGGCGAGAAACGGCAACTGAAAGAACTCGGTCGTGGCTATAAGCTGGGGCGGAATACCGGCCGCTTGGTGGCGGGGGCGCTGACGCTGGGTTTCGTTAGTCTGCCGTTGCTGGCGTGGATCGGCCAGTCGCCGATTATCATCATGACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGc                                                                                                                                                   <  1:637477/279‑1 (MQ=255)
                                                                                                                                                cATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTg                                                                                                                                                 <  2:827349/289‑1 (MQ=255)
                                                                                                                                                cATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTg                                                                                                                                                 <  2:711312/289‑1 (MQ=255)
                                                                                                                                                        aaaCCGCCACGCGTCTAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGCGCCCTGGCGCGGGGTTTCGTTAGTCTGCAGTTTCTCGCGTGGCTCGCCCCGTCGCCGCTTATCATCACTACCCGCTATCATTTTAGCCGCTGTCAATATGGATTGCTTCAACTTCCCCTTTTTGGTGAGTTTATTTCGGGTAACTTCCTTTACCGTCCTCTTACCTCCCCCCGCCCCGTCCGtttccgtgtttttat                                                                                                                                         >  1:871815/1‑276 (MQ=255)
                                                                                                                                                        aaaCCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGATGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCCCTGGGAATCGGTAGTCTGCCGTTGCTGGCGTGGAACGCCCAGTCGCCGATTACCATCATTACCGGCGAGGAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCTCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATtattat                                                                                                                                         >  2:254725/1‑289 (MQ=255)
                                                                                                                                                                      gcagggcggaAACGGTCGATGCAGGAAATGGGTCGGGATTATAGGCTGGTGCTGAAAGACGGCCGCGGTGTGGCGGGGGTGGTGCCGCTGGCATTCGTTAGTCTGCCGTTGCTGGCATTGGTCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGTAGCTATGCAGATGGCTTGCTGCAACTGCCTATTTGCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGa                                                                                                                          <  1:22167/281‑1 (MQ=255)
                                                                                                                                                                       taGGCGAGCACCTCTCACTGAGAGAGCGCGGTCGTGACTATAATCGGGTGCTGAAGCACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGTGATTCGTTATGCTGCCTGTTCGTGCGTGGTTTGCCCATTCGCCCGTTATCATCCTGAACGGCGAGCAGTTGAACAGCTATGCATATGGCTTGCGGTAAGGGCCTGGTGTCGGGGCCTTCAGTGAGGTTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGGCgct                                                                                                                         <  1:1330859/289‑3 (MQ=255)
                                                                                                                                                                             agaAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGGTTCGTTAGTCTGCCGTTGCTGGCGTGGAGCGCCCAGTCGCCGCTTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGGCTATTTTTGGGGCGTTTATTGCGGGTAACTTGCTTTTTGCGCGTCTTACCTCGCGCCGAACCGTAACTTCCCTGCTTATTATTGGGGGGTGTGCCGTTATGAt                                                                                                                   >  1:907075/1‑289 (MQ=255)
                                                                                                                                                                             agaAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGAGTTTTATGGGCGGCTGGCCGATTATGAt                                                                                                                   >  2:1159009/1‑289 (MQ=255)
                                                                                                                                                                              gaAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGAt                                                                                                                   <  1:415428/288‑1 (MQ=255)
                                                                                                                                                                               aaaCTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCt                                                                                                                                                                                                                                                     >  1:614574/1‑160 (MQ=255)
                                                                                                                                                                               aaaCTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCt                                                                                                                                                                                                                                                     <  2:614574/160‑1 (MQ=255)
                                                                                                                                                                                    gTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTTCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTa                                                                                                            >  2:738577/1‑289 (MQ=255)
                                                                                                                                                                                     tCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGGTTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTTCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCCCGCGCCGCACCGTACGTTCGCTGATTATTCTGGGCGGCTGGCCGCTTATGATTGGTCtgt                                                                                                           >  1:1030052/1‑287 (MQ=255)
                                                                                                                                                                                         tGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGt                                                                                                       <  2:474531/289‑1 (MQ=255)
                                                                                                                                                                                                     gTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGGTCGCCCAGTCGCCGAGTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCCAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGTCCTCGCGCCGCCCCGTACGTTCGCTGATTATTATGGGCGGGTGGCCGATTATGATTGGTTTATTTGTCGCTGCTGCGGc                                                                                          >  1:136769/1‑289 (MQ=255)
                                                                                                                                                                                                     gTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTGAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGTGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCgg                                                                                           <  2:161245/288‑1 (MQ=255)
                                                                                                                                                                                                      tCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCa                                                                                         <  2:246266/289‑1 (MQ=255)
                                                                                                                                                                                                       cGTGACGATAATCGGTGTCTGAAGAACGGCCGCTTGGTGGCGGGGGCGCTGGCGCTGGTATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCa                                                                                         <  1:109549/288‑1 (MQ=255)
                                                                                                                                                                                                              ttaatcgcgTGCTGAATAACGGCCGGAGTTGGTCGGGTGAGCTGGAGCTTGTCTTCGTTGGTCGGGCGTTGCGGGCGGGGTTCGCCCAGTGGCGGATTAGCATCATTCCCGGCGGGCACTTGAGCAGCTACGAATATGGGTTGCTGAAAGTGCCGATTTTCGGTGCGTTGATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTCCGCTGATTATTATGGGCGGCTTTCCGAGTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTAt                                                                                 <  1:1248399/281‑1 (MQ=255)
                                                                                                                                                                                                               taAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCCGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCCTTAATTGCGGGTAACTTGCTGTTTGCGCGTCTTACCTCGCCCCCCACCGGTCGTTCGCTGATTATTTTGGGCGGCTGGCCGGGTTTGATTGGGCTACTGGTCGCTGCTGGGGCAAAGGTTa                                                                                  >  1:1086672/1‑287 (MQ=255)
                                                                                                                                                                                                                aaGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGGTTCGTTAGTCTGCCGTTGCTGGCGTGGGTCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCCGCTATGAATATGGCTTGCTTCAAGGGCCAATTTTCGGGGTGGTAAGTTGGGGTAACTTGACGTTTGCGCGTCTGACATCGCGCCGCCCCGTTCGGTCGCTGATTATTATTGGTGGTTGGCGGTTTATTGTTGGTCTAATGGTCgcgcgttggcgaaggggttac                                                                                 >  1:843861/1‑269 (MQ=255)
                                                                                                                                                                                                                aaGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGGACGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGAACCCGCGCCGCACCGGACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGGTTGGTCTTTTGGGCgctgcttgg                                                                                             >  1:137312/1‑273 (MQ=255)
                                                                                                                                                                                                                     ggTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCgg                                                                                                                                                                                                    <  1:170580/170‑1 (MQ=255)
                                                                                                                                                                                                                     ggTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCgg                                                                                                                                                                                                    >  2:170580/1‑170 (MQ=255)
                                                                                                                                                                                                                       tGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGGAGTTGAGCCGCTATGAATATGGCTGGCTGCAAGTGCATATTTTCGGGGCGGTAAATGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTCCGTTCGCTGAGTATTAAGGGCGGCTGGCCGATTCTGGTTGGTCTTGTCGTCGCCGCTCGCGGAAGGGTTACCTCACCGc                                                                         >  1:818758/1‑288 (MQ=255)
                                                                                                                                                                                                                       tGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTcgcg                                                                                                                                                                    >  2:1337788/1‑199 (MQ=255)
                                                                                                                                                                                                                       tGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTcgcg                                                                                                                                                                    >  2:586578/1‑199 (MQ=255)
                                                                                                                                                                                                                       tGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTcgcg                                                                                                                                                                    <  1:1337788/199‑1 (MQ=255)
                                                                                                                                                                                                                           tGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCAcg                                                                      <  2:843861/288‑1 (MQ=255)
                                                                                                                                                                                                                                 aCGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTAt                                                                 <  2:1086672/287‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTTAGCAGCTATGAATATGGCTTGCTGCAAGGGCCTATTTTCGGGGGGTTAATTGAGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGTATGGTCTATTTGTTGCTGCTGCGGACACGGTTATCTCCTCGCACGCGTATTTa                                                              >  1:168225/1‑289 (MQ=255)
                                                                                                                                                                                                                                  cGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGAGCTCTCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTTCGGCAACGGTTATCTCAGCGCACGGGTAtt                                                                >  2:267408/1‑287 (MQ=255)
                                                                                                                                                                                                                                  cGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGc                                                                                          >  2:650231/1‑261 (MQ=255)
                                                                                                                                                                                                                                  cGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGc                                                                                          <  1:650231/261‑1 (MQ=255)
                                                                                                                                                                                                                                       gCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCTCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGAt                                                         >  2:1142616/1‑289 (MQ=255)
                                                                                                                                                                                                                                       gCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTTCAAGTTCCTATTTTCGGGGCGTTAATTTCGGGTAACTTGCTGTTATCGGTTCTGACCTCGCGCCGCACCGTTCCTTCCCTTATTATTATTGGCGGGCGGCCGGTTATGATTGGTCTATTGGTCGCTGCTGCGGCCAAGGTTAtctctgtccacccggttttttggat                                                         >  1:592205/1‑268 (MQ=255)
                                                                                                                                                                                                                                       gCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGAt                                                         <  2:1202208/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggggCGCTGGCGCTGGGAGTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCCATTTTCGGGGCGTTAATTGCGGGTTACTTTCTGTTAGCGCGTCTGACCTCGGGCCGCACCGTACGTTTGCTGATTATTATGGGGGGCTGGCCGATTTTGGTTGGGCTATTGGTCGCTGCTGCGGCAACGGGTTTGTCATCGCACGGGTATTTATGGAGGAATGTCgggg                                              >  1:483143/1‑288 (MQ=255)
                                                                                                                                                                                                                                                   gggCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCAcgcgg                                                                   >  2:445512/1‑266 (MQ=255)
                                                                                                                                                                                                                                                     gggcgggcgcttgGATTCGTTAGTCTGCCGTTGCTGGCGTAGATGGCCCAGTCGCCGATTATCATCATTACCGGCTAGCAGTTGAGCAGCTATGAATATGGCTTGCTTCAATTGCCTATTTTCGGGGCGTTAATTGCTGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGGACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTa                                            <  1:254725/283‑1 (MQ=255)
                                                                                                                                                                                                                                                     gcgctggcgctgGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCTCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTa                                            <  1:1142616/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gTTAGTCTGCCGGTGCTGGCGTGTATCGCCCAATCGCCAATTAACATCATTACCGGCTAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTTCTTATTTTCTGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCATCTGGCCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCTATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTa                         <  2:1044785/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                        ttAGTTTGCCGTTGTTGGCGTGGATCGCCCGGTGTCCTGTTATCATCATTGCCGGGGAGCCGTTGAGCAGCTATGAATATGGCTTGCTGAAATTGTCTATGTTCGGGGCGTTAATTGCGGGTAACTTGTTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGTTCTATTGGTCACTGCTGCGGCAACGGTTATCTCGTCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTa                         <  1:1183123/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCt                                >  1:127564/1‑280 (MQ=255)
                                                                                                                                                                                                                                                                         tAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCt                                <  2:127564/280‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ccGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTAt                                                                 >  1:1226086/1‑240 (MQ=255)
                                                                                                                                                                                                                                                                                 cGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGc               <  2:136769/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     gCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGTGGCAACGGTTATCTCATTGCACGCGTATTTATGGATGGCTGCCGGGTTAAGTAGTTATGGTTTCGGTATTGGTCTGGTGAAt           >  1:716267/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                       tGGCGTGGATCGCCCAGTCGGCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATg          <  2:871815/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          cGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTTTGGCGAATGCggg      >  1:38279/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                          cGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCggg      >  1:923399/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                            tGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGa     <  2:551058/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               aTCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTg  <  2:907075/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                     
CGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAGTTGAGCAGCTATGAATATGGCTTGCTGCAAGTGCCTATTTTCGGGGCGTTAATTGCGGGTAACTTGCTGTTAGCGCGTCTGACCTCGCGCCGCACCGTACGTTCGCTGATTATTATGGGCGGCTGGCCGATTATGATTGGTCTATTGGTCGCTGCTGCGGCAACGGTTATCTCATCGCACGCGTATTTATGGATGACTGCCGGGTTAAGTATTTATGCTTTCGGTATTGGTCTGGCGAATGCGGGACTG  >  NC_000913/884101‑884674

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