Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 773,145 T→C 5.4% L35P (CTC→CCC)  cydB → cytochrome d terminal oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913773,1450TC 5.4% 592.0 / 13.3 177L35P (CTC→CCC) cydBcytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base T (84/83);  new base C (5/5);  total (89/88)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

Only 100 of 179 total aligned reads displayed.

ATCTTCTCAATGGTGCTGATTTGCGGCCTGTATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATT  >  NC_000913/772871‑773416
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                               
gccttctCACTGGTGCTTGTTTGGTGCCTGTATACCCTGTTCCTGTTGGCAGAAAGGTTCTTAATTGTCGAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTACCACGACTACTCAGCCGGCACGCTAAGACAGTGGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCt                                                                                                                                                                                                                                                                                 <  4:368412/273‑1 (MQ=255)
  cttctCAATGGTGCTTGTTTGCGGCCTGGATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCa                                                                                                                                                                                                                                                                               <  1:124920/275‑1 (MQ=255)
        aaTGGTGCTGATTTGCGGCCTGTATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATTATCGATTATGAAGTATTCCGTTTTATCTGGTTGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTTACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGtt                                                                                                                                                                                                                                                                         >  2:264131/1‑275 (MQ=255)
            gTGCTGATTTGCGGCCTGTATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCt                                                                                                                                                                                                                                                                     <  2:79026/275‑1 (MQ=255)
               cTGATTTGCGGCCTGTATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATAGGGGTGGGCATGACTACCCGtttacccg                                                                                                                                                                                                                                                                   >  1:2855/1‑269 (MQ=255)
                tGATTTGCGGCCTGTATACCCTGTTCCTGGTGGCATAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGt                                                                                                                                                                                                                                                                 <  2:286956/275‑1 (MQ=255)
                       cGGCCTGTATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTGTTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCt                                                                                                                                                                                                                                                                     >  2:220207/1‑264 (MQ=255)
                       cGGCCTGTATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACg                                                                                                                                                                                                                                                          <  2:256730/275‑1 (MQ=255)
                           cTGTATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTGTTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGt                                                                                                                                                                                                                                                                 <  4:220206/264‑1 (MQ=255)
                                      gTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATNGATTATGAAGTATTGCGTTTTATCTGGTTGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGAGTGNTTCGACATTGGGGTGGGGCTGCTCACCAGTNNCCTCGGGCGTCACGCCAACggagctag                                                                                                                                                                                                                                             >  3:65244/1‑266 (MQ=255)
                                        tCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                                                                                                                                                         >  2:359959/1‑275 (MQ=255)
                                                gCAGAATTGTTCTTAATCGTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTaa                                                                                                                                                                                                                                 <  2:477291/275‑1 (MQ=255)
                                                           ctataTGTTCAAGGTGGAACGCCCCGGCCCAAGCAGCCGGAAAACCGGTCGCTATAACTTTTAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTGAAATGATCGATTATGAAGTATTGCGTTTTATCGGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGca                                                                                                                                                                                                                       <  4:375601/270‑1 (MQ=255)
                                                                tGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTCTTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTAGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTCACGACACCGAGAGTCGAATTATAATTAACTCCAGTgaccaacca                                                                                                                                                                                                                 >  1:160310/1‑267 (MQ=255)
                                                                   tCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCGCCCGTTTCCTCGGTCGCAACGCCACCGAGCGTCGAATTATGATTAACCCCATTGCACCCCACTgg                                                                                                                                                                                                              >  4:78735/1‑275 (MQ=255)
                                                                                                                     tttggtcAGTCTTCCACGACTACTCAGCCGGCAAGGTAAGACAGGAGTCGTCCAATGATCGTTTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtgctgc                                                                                                                                                            <  4:506817/269‑1 (MQ=255)
                                                                                                                         aGCAGTCTTCCACGTCTACTCAGCCGGCACGCTAAGACAGGACTCGTCAAATGATCGATTATGAAGTATTGCGTTTTCACGGGTGGCTGCTGGTAGGCCTTCTGCTGATTGGTTATGCAGTCACTGACGGTTTCGACATGGGGGTTGGCATGCTCACCCGTTTCCTAGGTCGTACCGATACCAAGCGTCGAATTATGATTAACTCCATTGCACCAACATGGGACGGTAACCAGGTTTGGCTGCACACCGCGGGCGGCCGACTCTTTGCTTCCTg                                                                                                                                                         >  4:442799/1‑274 (MQ=255)
                                                                                                                         aGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATCGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGCACGGTAACCAGGTTTGGCTGACCCCCGCGGCCGGCCCCCTCTTTGCTGCCTg                                                                                                                                                         >  4:281880/1‑274 (MQ=255)
                                                                                                                              tCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTtg                                                                                                                                                                >  2:506818/1‑262 (MQ=255)
                                                                                                                              tCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCGTTGCTGCCTGGCCGa                                                                                                                                                    >  1:100623/1‑274 (MQ=255)
                                                                                                                               cTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCg                                                                                                                                                                                                                                             >  2:430647/1‑184 (MQ=255)
                                                                                                                                    aCGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAccc                                                                                                                                                                                                                                                                            <  7:547340/148‑1 (MQ=255)
                                                                                                                                        ctactCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGttt                                                                                                                                                                                                                                                                        <  6:475335/148‑1 (MQ=255)
                                                                                                                                              aGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCgg                                                                                                                                                                                                                                                                  >  8:394028/1‑148 (MQ=255)
                                                                                                                                                  ggCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTtct                                                                                                                               <  1:65244/275‑1 (MQ=255)
                                                                                                                                                    cACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGAGCATGCCCAccccca                                                                                                                                                                                                                                                                         >  5:950285/1‑132 (MQ=255)
                                                                                                                                                      cGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACg                                                                                                                                                                                                                                                          <  8:793178/148‑1 (MQ=255)
                                                                                                                                                        cTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGac                                                                                                                                                                                                                                                        >  7:756569/1‑148 (MQ=255)
                                                                                                                                                             aCAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGa                                                                                                                                                                                                                                                   >  7:31716/1‑148 (MQ=255)
                                                                                                                                                                 gAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGCGTTGGCATGTTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGTTCACCGAGGTCTGCGCACTCTTTGCTGCCTTGCCGATGGTCTCTGCCGCTGCGTTTTCCTTCTTCTAtgtg                                                                                                                 >  2:321246/1‑274 (MQ=255)
                                                                                                                                                                  aGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCgg                                                                                                                                                                                                                                                                  >  8:1073711/1‑128 (MQ=255)
                                                                                                                                                                  aGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCgg                                                                                                                                                                                                                                                                  <  7:1073711/128‑1 (MQ=255)
                                                                                                                                                                   gtcgtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGCCACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGTGACGGTAACCAGGTTTGGCTGATCACCTCGGGCGGCGCACTCTTTGCTGCCTGGCCGAAGGTCTATGCCCCTGCGTTCTCCGGCCTCTAttgt                                                                                                                 >  4:504272/1‑269 (MQ=255)
                                                                                                                                                                   gtcgtcAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCga                                                                                                             <  2:136564/275‑1 (MQ=255)
                                                                                                                                                                        cAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTa                                                                                                                                                                                                                                        <  8:329022/148‑1 (MQ=255)
                                                                                                                                                                          aaTGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATg                                                                                                                                                                                                                                      <  5:288938/148‑1 (MQ=255)
                                                                                                                                                                               tCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGt                                                                                                                                                                                                                                                                 <  5:688800/116‑1 (MQ=255)
                                                                                                                                                                               tCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGa                                                                                                                                                                                                                                     <  8:471464/144‑1 (MQ=255)
                                                                                                                                                                                cGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTa                                                                                                                                                                                                                                                             <  8:926698/119‑1 (MQ=255)
                                                                                                                                                                                cGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTa                                                                                                                                                                                                                                                             >  7:926698/1‑119 (MQ=255)
                                                                                                                                                                                cGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGa                                                                                                                                                                                                                                     <  6:384988/143‑1 (MQ=255)
                                                                                                                                                                                  aTTATGAAGTATTGCGTTTTATCTGGTTGCTGCTGGTTGGCGTTCTGCTGATTGGTGTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGaa                                                                                                                                                                                                                                           <  6:597898/135‑1 (MQ=255)
                                                                                                                                                                                          gTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGcac                                                                                                                                                                                                                      >  5:565085/1‑148 (MQ=255)
                                                                                                                                                                                              tGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGt                                                                                                                                                                                                                                                                 >  7:803076/1‑101 (MQ=255)
                                                                                                                                                                                              tGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATg                                                                                                                                                                                                                                      >  7:264242/1‑128 (MQ=255)
                                                                                                                                                                                              tGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATg                                                                                                                                                                                                                                      <  8:264242/128‑1 (MQ=255)
                                                                                                                                                                                              tGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCa                                                                                                                                                                                                                            <  6:61098/138‑1 (MQ=255)
                                                                                                                                                                                                 gTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTg                                                                                                                                                                                                               >  7:284363/1‑148 (MQ=255)
                                                                                                                                                                                                   tttATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCATCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTggg                                                                                                                                                                                                             <  8:789356/148‑1 (MQ=255)
                                                                                                                                                                                                   tttATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTAt                                                                                                                                                                                                                                       <  6:852059/122‑1 (MQ=255)
                                                                                                                                                                                                            tggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAg                                                                                                                                                                                                                                                  <  7:647480/102‑1 (MQ=255)
                                                                                                                                                                                                            tggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAg                                                                                                                                                                                                                                                  >  8:647480/1‑102 (MQ=255)
                                                                                                                                                                                                            tggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAAcc                                                                                                                                                                                                    <  8:1041099/148‑1 (MQ=255)
                                                                                                                                                                                                            tggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAAcc                                                                                                                                                                                                    <  5:841955/148‑1 (MQ=255)
                                                                                                                                                                                                             ggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGt                                                                                                                                                                                                                                                                 <  6:789516/86‑1 (MQ=255)
                                                                                                                                                                                                             ggCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGt                                                                                                                                                                                                                                                                 >  5:789516/1‑86 (MQ=255)
                                                                                                                                                                                                               ctgctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACCCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGCTGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGtt                                                                   >  2:374797/1‑274 (MQ=255)
                                                                                                                                                                                                                 gctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCg                                                                                                                                                                                                                                                >  5:640105/1‑99 (MQ=255)
                                                                                                                                                                                                                 gctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGtt                                                                                                                                                                                               <  7:394028/148‑1 (MQ=255)
                                                                                                                                                                                                                  ctgGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGttt                                                                                                                                                                                              <  8:304006/148‑1 (MQ=255)
                                                                                                                                                                                                                    ggTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTgg                                                                                                                                                                                            <  6:802960/148‑1 (MQ=255)
                                                                                                                                                                                                                    ggTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTgg                                                                                                                                                                                            <  8:284363/148‑1 (MQ=255)
                                                                                                                                                                                                                    ggTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTgg                                                                                                                                                                                            >  7:678962/1‑148 (MQ=255)
                                                                                                                                                                                                                      ttGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCa                                                                                                                                                                                                   <  8:603305/139‑1 (MQ=255)
                                                                                                                                                                                                                         gCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTggg                                                                                                                                                                                                             >  8:347057/1‑126 (MQ=255)
                                                                                                                                                                                                                         gCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTggg                                                                                                                                                                                                             <  7:347057/126‑1 (MQ=255)
                                                                                                                                                                                                                             tCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTaa                                                                                                                                                                                                                                 <  6:680379/102‑1 (MQ=255)
                                                                                                                                                                                                                               tgctgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGt                                                                                                                                                                                                                                               >  4:109077/1‑86 (MQ=255)
                                                                                                                                                                                                                                gctgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGa                                                                                                                                                                                                                                     <  5:601811/95‑1 (MQ=255)
                                                                                                                                                                                                                                gctgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGa                                                                                                                                                                                        <  8:131473/140‑1 (MQ=255)
                                                                                                                                                                                                                                  tgATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACAcc                                                                                                                                                                                                                                                     >  5:453357/1‑77 (MQ=255)
                                                                                                                                                                                                                                     ttGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCgcgggcgg                                                                                                                                                                           <  8:678962/148‑1 (MQ=255)
                                                                                                                                                                                                                                      tGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGt                                                                                                                                                                                                                                               <  7:916621/79‑1 (MQ=255)
                                                                                                                                                                                                                                         ttttGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAAtt                                                                                                                                                                                                                                         >  8:8079/1‑82 (MQ=255)
                                                                                                                                                                                                                                            ttcctgcagtGACGGTTTCGGAAGGGGGGGGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTGGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAAccc                                     <  4:359958/266‑1 (MQ=255)
                                                                                                                                                                                                                                               aGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCtt                                                                                                                                                                  >  5:927084/1‑147 (MQ=255)
                                                                                                                                                                                                                                                 tCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGt                                                                                                                                                                                                                                                                          >  7:444883/1‑41 (MQ=255)
                                                                                                                                                                                                                                                     tGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGcc                                                                                                                                                           >  7:856840/1‑148 (MQ=255)
                                                                                                                                                                                                                                                      gACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTa                                                                                                                                                                                                       <  6:949782/103‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCGACATGGGGGTGGGCATGCTCCCCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACcgc                                                                                                                                                                                 <  8:700109/118‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACcgc                                                                                                                                                                                 >  7:700109/1‑118 (MQ=255)
                                                                                                                                                                                                                                                             tCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGAt                                                                                                                                                   <  6:454727/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACcg                                                                                                                                                                                  <  8:124599/114‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACcg                                                                                                                                                                                  >  7:124599/1‑114 (MQ=255)
                                                                                                                                                                                                                                                                   tGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACCCTGGGCCGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCCTTGCAGCCTGGCTGATGGTCTACGCCCCTGCGTTCTCCGGCTTCTATGTGGCGATGACCATCGTGCTCGCGACTTTGTTCCTCCTTCAGGGTGCTTTTTACTACCGCTCCCTTTTTGAAGTCCCCCCATGGCGTAAACTGTggtgttgg              >  4:408891/1‑269 (MQ=255)
                                                                                                                                                                                                                                                                    gggggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGAt                                                                                                                                                                                                                                    >  5:156742/1‑60 (MQ=255)
                                                                                                                                                                                                                                                                    gggggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGAt                                                                                                                                                                                                                                    >  7:105924/1‑60 (MQ=255)
                                                                                                                                                                                                                                                                    gggggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGAt                                                                                                                                                                                                                                    <  6:156742/60‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ggggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCg                                                                                                                                                     <  6:41151/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ggggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCg                                                                                                                                                     >  5:41151/1‑138 (MQ=255)
                                                                                                                                                                                                                                                                     cgggTGAGCATGCCCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTg                                                                                                                                                                                         >  7:737057/2‑102 (MQ=255)
                                                                                                                                                                                                                                                                      gggTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTcc                                                                                                                                                                                                                             <  7:285267/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGtt                                                                                                                                                                                               >  5:958323/1‑93 (MQ=255)
                                                                                                                                                                                                                                                                        gTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGtt                                                                                                                                                                                               <  6:958323/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                          gggCATGCTCACCCGTTTCCTCGGTCGTAACGAAACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCt       <  1:5247/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ggCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCCAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCtt      >  2:414057/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                              aTGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTggg                                                                                                                                                                                                             <  8:42839/73‑1 (MQ=255)
                                                                                                                                                                                                                                                                              aTGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTggg                                                                                                                                                                                                             >  7:42839/1‑73 (MQ=255)
                                                                                                                                                                                                                                                                               tGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGtt                                                                                                                                                                                               >  8:871538/1‑86 (MQ=255)
                                                                                                                                                                                                                                                                                  tCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTtctc                                                                                                                              <  8:739680/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                               
ATCTTCTCAATGGTGCTGATTTGCGGCCTGTATACCCTGTTCCTGGTGGCAGAATTGTTCTTAATGTTCAAGTTTGCACGCCTCGGCCCAAGCAGCCTGAAAACCGGTCGCTATCACTTTGAGCAGTCTTCCACGACTACTCAGCCGGCACGCTAAGACAGGAGTCGTCAAATGATCGATTATGAAGTATTGCGTTTTATCTGGTGGCTGCTGGTTGGCGTTCTGCTGATTGGTTTTGCAGTCACTGACGGTTTCGACATGGGGGTGGGCATGCTCACCCGTTTCCTCGGTCGTAACGACACCGAGCGTCGAATTATGATTAACTCCATTGCACCACACTGGGACGGTAACCAGGTTTGGCTGATCACCGCGGGCGGCGCACTCTTTGCTGCCTGGCCGATGGTCTATGCCGCTGCGTTCTCCGGCTTCTATGTGGCGATGATCCTCGTGCTGGCGTCTTTGTTCTTCCGTCCGGTCGGTTTTGACTACCGCTCCAAGATTGAAGAAACCCGCTGGCGTAACATGTGGGACTGGGGCATCTTCATT  >  NC_000913/772871‑773416

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: