Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,311,616 C→T 6.4% intergenic (‑82/+30) apbE ← / ← ompC putative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,311,6160CT 6.4% 929.6 / 35.5 289intergenic (‑82/+30)apbE/ompCputative thiamine‑synthetic flavin transferase lipoprotein/outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (140/130);  new base T (12/7);  total (152/137)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-02

Only 100 of 292 total aligned reads displayed.

TTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGTACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCG  >  NC_000913/2311346‑2311886
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                              
    ttAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGTACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCCGCAGTGCCACCCGGGTAAAGCTTATTTCCCTACGTGAGTATCCTCATGCGAACGGTCGCAAGATTACACCAAAAAACTGAGTTTTTCAGCTGAACACAATGAAAAAAGGGCCCgagagccc                                                                                                                                                                                                                                                                        >  4:465945/1‑268 (MQ=255)
                 ggcgttggggttaaaTGCGCGGGCTGCTGGGACNNNAGTAGGGACCCGTCGTTTTGCCATCAAGAACGGTAAGTTCAGTAGGTTGGGTTTTGGGGGGTTGTGGTTGTTGATCGCAACCAAGAAAGAAGAGGGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATCCGTGATGATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATc                                                                                                                                                                                                                                                       <  3:265692/266‑1 (MQ=255)
                       gggtgtaaaTGCCCGGGGTGCTGGGAAGCATGAAGGTACCCAAAGTTTTGTCGTCAAGAGCGGTAAGTACAAAAGCGTGGGTTTTGGGGGGTGGTAGTTTTGAGTCGCAACCGACAAAGAAGAGCGCGGCAGCCAGCAGTACCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGATCGCAAGAGGACACCAAAAAACTGAGTGTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATc                                                                                                                                                                                                                                                       <  4:192720/266‑1 (MQ=255)
                                              gCACGCCAGAAGGTACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGAcc                                                                                                                                                                                                                              <  2:644498/275‑1 (MQ=255)
                                                               gtaaGTTTTGCCTTCCAGTACAGTAGATTCGGTAGCTGGGTGTTTGGCGGGTTGTTCTGTTTTCTCGCAACCAACAAAGAAGAGCGCGGAAACCAGCAATGCGTGGCGGGTAAAGCTTATTTCCACACGGGATTACCCTCAGGCGAAGGGTCGCAAGAGTCGACCAAAGAACTAAGTTTGGACGCTGAAAACAAACAAAAAAGGGCCCGGAGGCCCTTTGTTCGATATCAATAGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGCTGTTAt                                                                                                                                                                                                             <  4:613913/272‑1 (MQ=255)
                                                                                                             ggCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagca                                                                                                                                                                >  2:547655/1‑275 (MQ=255)
                                                                                                                                tCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAgg                                                                                                                                                                                                                                                                           >  6:1517788/1‑149 (MQ=255)
                                                                                                                                tCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTAGCAGTGTTGAGGCCAGCATAGCGAGTGAACGTGTTGTCGTCCAGACGGTTGATAGTGTAGGCAAc                                                                                                                                             >  4:302208/1‑275 (MQ=255)
                                                                                                                                       caacaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCTTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCtttt                                                                                                                                                                                                                                                                    <  6:1042950/149‑2 (MQ=255)
                                                                                                                                         acaaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGc                                                                                                                                                                                                                                                                          <  5:436522/141‑1 (MQ=255)
                                                                                                                                          caaAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTc                                                                                                                                                                                                                                                                 <  7:1608386/149‑1 (MQ=255)
                                                                                                                                           aaaGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCg                                                                                                                                                                                                                                                                >  7:904747/1‑149 (MQ=255)
                                                                                                                                             agaagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGAGCCCGCAGGCCCTTTGTTCGat                                                                                                                                                                                                                                                              >  5:1303340/1‑149 (MQ=255)
                                                                                                                                              gaagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGata                                                                                                                                                                                                                                                             <  6:835313/149‑1 (MQ=255)
                                                                                                                                               aagaGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAAAGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGatat                                                                                                                                                                                                                                                            <  7:1234891/149‑1 (MQ=255)
                                                                                                                                                agagCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTATTCGATATc                                                                                                                                                                                                                                                           >  7:1422531/1‑149 (MQ=255)
                                                                                                                                                  agCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGNCCGCAGGCCCTTTGTTCGATATCNATCCNGATTAGAACTGGTAAACCAGAccc                                                                                                                                                                                                                             >  4:701817/1‑177 (MQ=255)
                                                                                                                                                  agCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttttt                                                                                                                           >  2:492008/1‑275 (MQ=255)
                                                                                                                                                  agCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGttttt                                                                                                                           >  2:153343/1‑275 (MQ=255)
                                                                                                                                                    cgcgGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATc                                                                                                                                                                                                                                                       <  5:596757/149‑1 (MQ=255)
                                                                                                                                                      cgGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCga                                                                                                                                                                                                                                                     <  7:914224/149‑1 (MQ=255)
                                                                                                                                                       ggCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACGCCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCgag                                                                                                                                                                                                                                                    <  5:56228/149‑1 (MQ=255)
                                                                                                                                                       ggCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGa                                                                                                                      <  4:547651/275‑1 (MQ=255)
                                                                                                                                                        gCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCgaga                                                                                                                                                                                                                                                   >  8:815146/1‑149 (MQ=255)
                                                                                                                                                         cagccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGaa                                                                                                                     >  2:202025/1‑274 (MQ=255)
                                                                                                                                                           gccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTa                                                                                                                                                                                                                                                >  5:793191/1‑149 (MQ=255)
                                                                                                                                                            ccagcAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAg                                                                                                                                                                                                                                               >  6:713120/1‑149 (MQ=255)
                                                                                                                                                               gcagTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAAc                                                                                                                                                                                                                                            >  5:268260/1‑149 (MQ=255)
                                                                                                                                                                     ccACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACt                                                                                                                                                                                                                                           <  6:591665/144‑1 (MQ=255)
                                                                                                                                                                     ccACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTaa                                                                                                                                                                                                                                      <  5:757659/149‑1 (MQ=255)
                                                                                                                                                                       aCGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAAc                                                                                                                                                                                                                                    >  8:1348120/1‑149 (MQ=255)
                                                                                                                                                                         gcgGGTAAAGCTTATTTCCATACGTGATTATCCGCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTttgttg                                                                                                                       <  2:83380/256‑1 (MQ=255)
                                                                                                                                                                           gggTAAAGCTTATTTCCATACGGGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacat                                                                                                 <  1:270331/276‑1 (MQ=255)
                                                                                                                                                                             gTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGAcc                                                                                                                                                                                                                              >  8:75285/1‑149 (MQ=255)
                                                                                                                                                                              tAAAGCTTATTTCCATACGTGATTAGCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACcaacatcaac                                                                                             <  1:15198/277‑1 (MQ=255)
                                                                                                                                                                                  gCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCagag                                                                                                                                                                                                                         >  7:1379903/1‑149 (MQ=255)
                                                                                                                                                                                   cTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGc                                                                                                                                                                                                                        >  5:1140552/1‑149 (MQ=255)
                                                                                                                                                                                    ttATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCt                                                                                                                                                                                                                       >  7:1033892/1‑149 (MQ=255)
                                                                                                                                                                                       tttCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAACAGGGCCCGCATGCCCTTTGTTCGATATCATTCGAGATTAGAACTGGTAACCCAGACCCAGCGCTACGATGTTATAAGTTTTGATGCCAGCGTCACGAGTGAACTTGTTGTCGTCCAGCAGGTTGTTTTTGTAGTCAACTTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACACACATCAACTTAtt                                                                                        >  4:341576/1‑273 (MQ=255)
                                                                                                                                                                                        ttCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGa                                                                                                                                                                                                                   <  5:907424/149‑1 (MQ=255)
                                                                                                                                                                                          ccATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGTCCCAGAGCTACGATg                                                                                                                                                                                                                 >  6:1345017/1‑149 (MQ=255)
                                                                                                                                                                                          ccATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATg                                                                                                                                                                                                                 <  5:1600097/149‑1 (MQ=255)
                                                                                                                                                                                            aTACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCTGCGGGCCCttttt                                                                                                                                                                                                                                                                   >  7:780227/1‑95 (MQ=255)
                                                                                                                                                                                             tACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATc                                                                                                                                                                                                                                                           >  2:339753/1‑104 (MQ=255)
                                                                                                                                                                                             tACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTa                                                                                                                                                                                                              >  8:1592506/1‑149 (MQ=255)
                                                                                                                                                                                                gTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCa                                                                                                                                                                                                           <  5:1035924/149‑1 (MQ=255)
                                                                                                                                                                                                 tGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAGGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTa                                                                                                                                                                                                              <  5:190142/145‑1 (MQ=255)
                                                                                                                                                                                                  gATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAgt                                                                                                                                                                                                         <  8:1432550/149‑1 (MQ=255)
                                                                                                                                                                                                        ccTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGAt                                                                                                                                                                                                   <  8:1650639/149‑1 (MQ=255)
                                                                                                                                                                                                         cTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAAc                                                                                                                                                                                                                                    <  7:1212815/115‑1 (MQ=255)
                                                                                                                                                                                                          tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                                                                                                                                 >  7:1542674/1‑149 (MQ=255)
                                                                                                                                                                                                          tCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGc                                                                                                                                                                                                 >  6:198397/1‑149 (MQ=255)
                                                                                                                                                                                                            aTGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGAt                                                                                                                                                                                                                                                  >  2:699538/1‑98 (MQ=255)
                                                                                                                                                                                                                     ttCGCAAGAGTACACCAAAAAACTGAGTTTGTAAGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGAccc                                                        <  2:416104/274‑1 (MQ=255)
                                                                                                                                                                                                                     gTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGa                                                                                                                                                                                      <  6:526228/149‑1 (MQ=255)
                                                                                                                                                                                                                      tCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGatat                                                                                                                                                                                                                                                            <  1:241536/78‑1 (MQ=255)
                                                                                                                                                                                                                        gCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGaa                                                                                                                                                                                                                                             >  7:188850/1‑91 (MQ=255)
                                                                                                                                                                                                                        gCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGaa                                                                                                                                                                                                                                             <  8:188850/91‑1 (MQ=255)
                                                                                                                                                                                                                        gCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTg                                                                                                                                                                                   >  6:787710/1‑149 (MQ=255)
                                                                                                                                                                                                                        gCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTg                                                                                                                                                                                   >  6:872160/1‑149 (MQ=255)
                                                                                                                                                                                                                         cAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGt                                                                                                                                                                                    >  7:166440/1‑147 (MQ=255)
                                                                                                                                                                                                                            gagTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACt                                                                                                                                                                               <  7:815146/149‑1 (MQ=255)
                                                                                                                                                                                                                              gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGt                                                                                                                                                                                                                                        <  3:55389/90‑1 (MQ=255)
                                                                                                                                                                                                                              gTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTGGATTTTg                                                                                                                                                     >  2:738003/1‑173 (MQ=255)
                                                                                                                                                                                                                               tACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCa                                                                                                                                               >  4:738630/1‑178 (MQ=255)
                                                                                                                                                                                                                                acacCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGtt                                                                                                                                                                           <  8:1265073/149‑1 (MQ=255)
                                                                                                                                                                                                                                acacCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTAcc                                             >  2:52478/1‑275 (MQ=255)
                                                                                                                                                                                                                                  acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGatat                                                                                                                                                                                                                                                            <  5:1430486/66‑1 (MQ=255)
                                                                                                                                                                                                                                  acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                                                                                                                                                         >  7:721672/1‑149 (MQ=255)
                                                                                                                                                                                                                                  acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                                                                                                                                                         >  6:30034/1‑149 (MQ=255)
                                                                                                                                                                                                                                  acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                                                                                                                                                         >  8:350333/1‑149 (MQ=255)
                                                                                                                                                                                                                                  acCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgt                                                                                                                                                                         >  6:646369/1‑149 (MQ=255)
                                                                                                                                                                                                                                    cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTaaa                                                                                                                                                                                                                                     >  7:563557/1‑87 (MQ=255)
                                                                                                                                                                                                                                    cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTaaa                                                                                                                                                                                                                                     <  8:563557/87‑1 (MQ=255)
                                                                                                                                                                                                                                    cAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTgtcg                                                                                                                                                                       >  6:1639059/1‑149 (MQ=255)
                                                                                                                                                                                                                                       aaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTa                                                                                                                                                                                                                                       >  4:678726/1‑82 (MQ=255)
                                                                                                                                                                                                                                       aaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTcc                                                                                                                                                                    <  5:1517788/149‑1 (MQ=255)
                                                                                                                                                                                                                                       aaaaCTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGAc                                      <  4:153341/275‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAAcc                                                                                                                                                                                                                                   <  2:678731/82‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagca                                                                                                                                                                >  5:652833/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tGAGTTTGTACGCTGAAAACCATGAAAAATGGGCCCTCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGGGAACTGGTTGTCGTCcagcag                                                                                                                                                               >  6:1210171/1‑149 (MQ=255)
                                                                                                                                                                                                                                            tGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCcagcag                                                                                                                                                               >  7:1107526/1‑149 (MQ=255)
                                                                                                                                                                                                                                               gTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGtt                                                                                                                                                            >  7:997623/1‑149 (MQ=255)
                                                                                                                                                                                                                                               gTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGtt                                                                                                                                                            >  8:937258/1‑149 (MQ=255)
                                                                                                                                                                                                                                               gTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGtt                                                                                                                                                            >  5:1470418/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     aCGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGAGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGAtttt                                                                                                                                                      <  6:1470418/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       gCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGt                                                                                                                                                    <  8:702369/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGt                                                                                                                                           >  6:1420843/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aaTGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGt                                                                                                                                           >  5:311793/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTTTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAg                                                                                                                                         >  5:1678663/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   gAACAAAGGGCCTGCGGGCCCTTTTTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAgg                                                                                                                                        >  6:397602/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                                                                                                                                             <  7:452586/111‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aaaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTgg                                                                                                                                                                             >  8:452586/1‑111 (MQ=255)
                                                                                                                                                                                                                                                                     acaaaGGGCCTGCGGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg                                                                                                                                      >  7:233454/3‑149 (MQ=255)
                                                                                                                                                                                                                                                                     aaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACg                                                                                                                                            <  6:800455/143‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaaaaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTg                                                                                                                                      >  7:1173729/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       aaaGGGCCTGCGGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGa                                                                                                                                    >  6:565554/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        aaGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGAc                                                                                                                                   <  8:1107526/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                           ggCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCg  <  2:302212/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                             cccGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGa                                                                                                                                                                                      >  8:1196760/1‑93 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                              
TTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGTACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGCGAACGGTCGCAAGAGTACACCAAAAAACTGAGTTTGTACGCTGAAAACAATGAAAAAAGGGCCCGCAGGCCCTTTGTTCGATATCAATCGAGATTAGAACTGGTAAACCAGACCCAGAGCTACGATGTTATCAGTGTTGATGCCAGCGTCACGAGTGAACTGGTTGTCGTCCAGCAGGTTGATTTTGTAGTCAACGTAGGTGGACATGTTTTTGTTGAAGTAGTAGGTAGCACCAACATCAACATATTTCAGGATATCTTCGTCGTCGTAGCCACGACCCAGGTTTTTACCTTTAGACTGCAGGTAAGCCAGGGACGGACGCAGACCGAAGTCG  >  NC_000913/2311346‑2311886

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