Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,602,509 C→T W117* (TGG→TAG)  yjjP ← putative succinate exporter YjjP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,602,5090CT97.4% 209.5 / ‑1.7 76W117* (TGG→TAG) yjjPputative succinate exporter YjjP
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (30/44);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

CGGCCATCGCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCG  >  NC_000913/4602264‑4602753
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
cggcCATCGCCATGGTGGGGGAATTGCTGAAAGTGGGGAGTTGCAAAAGCCATCCGGAAATGGTGGTGGCGGCGAAAGCGGTGAGGCAAAAGTTGATCTGTGGAGGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTAACGAAACAAGCGCAAGAAATGCCAACCGTTAAGGATACCCGCTATCtt                                                                                                                                                                                                                                                 >  2:133870/1‑251 (MQ=255)
       cGCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAAc                                                                                                                                                                                                                                          <  1:135068/251‑1 (MQ=255)
             ggtggGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAAt                                                                                                                                                                                                                                     <  1:283676/250‑1 (MQ=255)
                   ggaaTTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGaa                                                                                                                                                                                                                              <  1:252485/248‑1 (MQ=255)
                      aTTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAtt                                                                                                                                                                                                                            <  2:172835/250‑1 (MQ=255)
                            gAAAGTGGGGAGTTGCAAAAGCACTCCGGAAATGGTGGCGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGCGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCACCCATTCAGGCTACCAGCTATCTTGGGTAACGTAATGGATGATTTTgactc                                                                                                                                                                                                                      >  2:232530/1‑246 (MQ=255)
                              aaGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGGGATGACGGAACCATCCCAGCCACCGTTATTGAGTTTACAGAAAAAGGCGCAAGAAAGGCCAACCATTAAGGCTAACAGCTATCTTGGGTAACGTAAAGGCTGAATTTGGCTaaa                                                                                                                                                                                                                    >  2:391510/1‑250 (MQ=255)
                                     gAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGCCAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTtct                                                                                                                                                                                                            >  1:156729/1‑251 (MQ=255)
                                       gTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAGAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCCAAGACGGTGATGACGGCACTATCACAGCCACCGTTATTGCGTTTCAAGACACAGGCGCAAGCAAGGCCAAACATTAAGTCAACCAACTATCTTGGGTAACGTAATGGCTGAATTTGGCAAACTCGTTTctct                                                                                                                                                                                                          >  2:473252/1‑251 (MQ=255)
                                         tGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTa                                                                                                                                                                                                         <  1:169188/250‑1 (MQ=255)
                                          gCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCCGCCACCGTTATTGAGTTTACAGAAACAGGAGCAAGAAAGGCCAACCATTAAGGCTACCAGATATATTGGGTAACGTAATGGCTGAATTTGGATAAATCGTTtcttatcg                                                                                                                                                                                                       >  2:4968/1‑246 (MQ=255)
                                          gCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTGCCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTAc                                                                                                                                                                                                        <  1:68640/250‑1 (MQ=255)
                                          gCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACg                                                                                                                                                                                                       <  1:395731/251‑1 (MQ=255)
                                            aaaaGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGATACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGcc                                                                                                                                                                                                     <  1:300602/251‑1 (MQ=255)
                                             aaaGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCt                                                                                                                                                                                                    >  1:217699/1‑251 (MQ=255)
                                             aaaGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGAGGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCt                                                                                                                                                                                                    <  2:276911/251‑1 (MQ=255)
                                              aaGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCtt                                                                                                                                                                                                   >  1:493790/1‑251 (MQ=255)
                                                  aaTCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGAAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACAAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGt                                                                                                                                                                                                >  2:198814/1‑250 (MQ=255)
                                                        gAAATGGTGGTGGCTGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGTAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCcag                                                                                                                                                                                         <  2:288191/251‑1 (MQ=255)
                                                          aaTGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGGTTCTCTACGCCTTTGTAATCcagca                                                                                                                                                                                       <  1:432030/251‑1 (MQ=255)
                                                          aaTGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGGGGATGAAGATGACGTGGTGCCAGCAGCTGGCGGATATACATCGCGGGCGTACTGGCAAAGAAGGGGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACTCCTTTGTAATCcagca                                                                                                                                                                                       <  2:306321/251‑1 (MQ=255)
                                                              gtggtggCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatg                                                                                                                                                                                   <  1:479443/251‑1 (MQ=255)
                                                                ggtggCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAgatgat                                                                                                                                                                                 <  1:536414/251‑1 (MQ=255)
                                                                        cgAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCa                                                                                                                                                                         <  1:27239/251‑1 (MQ=255)
                                                                        cgAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACAGCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCa                                                                                                                                                                         <  1:40824/251‑1 (MQ=255)
                                                                         gAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGTAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCaa                                                                                                                                                                        <  1:528263/251‑1 (MQ=255)
                                                                               ggTAAGGCAAAAGTTGATATGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAAGCCAGCAGATGATGCTCCGCAAGaaaca                                                                                                                                                                  <  2:295002/251‑4 (MQ=255)
                                                                                 gaaGGCAAAAGTTGATCTGTGGATGAAGAAGGCGGTGTGGCAGCAGCTGGGGGATATACAACGCGGTCGTACTGGCAAAGAAAGTGATGACGGCACCAATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCACAAGAATcac                                                                                                                                                                 <  2:109445/250‑1 (MQ=255)
                                                                                  aaGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCGTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATcaca                                                                                                                                                                <  1:251183/250‑1 (MQ=255)
                                                                                    ggCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGAAAAGAAGGGGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtg                                                                                                                                                             >  1:299812/1‑251 (MQ=255)
                                                                                      cAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAATGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAAtgtg                                                                                                                                                           <  2:493790/251‑1 (MQ=255)
                                                                                           gTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGGGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGcgggg                                                                                                                                                      >  2:404804/1‑247 (MQ=255)
                                                                                             ggatttgtgGATGAAGATGAGGTTGTGCCAGCAGCTGGCGGATATACAGCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTAGTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCATCTAGCTTGGTTAACGTAATGGCTGAATTTGGCTAAATCGTTTCCCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGGATCACAATGTGCCGGACt                                                                                                                                                    <  2:514564/246‑1 (MQ=255)
                                                                                                 cTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGTTGCTCAGCAAGAATCACAATGTGCTGGACTTc                                                                                                                                                  >  2:161747/1‑249 (MQ=255)
                                                                                                             aTGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGc                                                                                                                                     >  1:270285/1‑250 (MQ=255)
                                                                                                               gACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGAGGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCata                                                                                                                                  <  1:133870/251‑1 (MQ=255)
                                                                                                                  gTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTa                                                                                                                               <  1:272631/251‑1 (MQ=255)
                                                                                                                        ccAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGCGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCCACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCcg                                                                                                                         >  1:510950/1‑251 (MQ=255)
                                                                                                                           gcagcTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCcgcg                                                                                                                       >  1:234130/1‑250 (MQ=255)
                                                                                                                            cagcTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGAt                                                                                                                     <  1:479358/251‑1 (MQ=255)
                                                                                                                                     gATATACATAGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCAACCGTTCTTGAGTTTACAGCAACAGGCGCCAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTACTCCAGCAGCTGATGCTCCGCAAGAATCACAATGTGCTGGACTTCCGTCCCCACATGCATATAAATGcccagttcgtgctttt                                                                                                            >  2:544102/1‑237 (MQ=255)
                                                                                                                                          acctCGCGGTCGTACTGGCACAGAAGGTGATGAAGGTACCCACCCAGCCACCGTTATTGAGTGTACAGAAACAGGCTCAAGAAAGGCCAACCATTAAGGCTACCAGCTATTTTGGGTAACGTAACGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTACTCCAGCAGATGATGCTCTGCAAGAATCATAATGTGCTGTACTTCAGGCACCGCATGCATAGTAATGCCGCGATCGTGAttatcacg                                                                                                        <  2:29863/247‑7 (MQ=255)
                                                                                                                                          aCATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACgt                                                                                                       <  1:285769/251‑1 (MQ=255)
                                                                                                                                          aCATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGACAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGt                                                                                                                                                                                                >  2:394610/1‑162 (MQ=255)
                                                                                                                                          aCATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGACAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCGTTGt                                                                                                                                                                                                <  1:394610/162‑1 (MQ=255)
                                                                                                                                              cccgGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGTCTTCAGTCACCACATGCAGATTAATGCCGCGATCGTGATTTTTACgtgt                                                                                                     <  2:99080/247‑1 (MQ=255)
                                                                                                                                               gcgGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCTCCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAAGGCCGCGATCGTGATTTTTACGTGTCGa                                                                                                  <  1:404267/251‑1 (MQ=255)
                                                                                                                                                cgGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTACTGCCGCGATCGTGTTTTTTACGGGTCgct                                                                                                 >  2:505053/1‑249 (MQ=255)
                                                                                                                                                cgGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGAt                                                                                                 <  2:270285/251‑1 (MQ=255)
                                                                                                                                                    cGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTc                                                                                              <  1:447810/250‑1 (MQ=255)
                                                                                                                                                       aCTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCGTGGGTAACGTAATGGCTTAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATTATGCTCCGCAAGAATCACAATGTGCTTGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGc                                                                                          <  1:527547/251‑1 (MQ=255)
                                                                                                                                                          ggCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAAGTGGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCAtt                                                                                       <  1:513778/251‑1 (MQ=255)
                                                                                                                                                            cAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGc                                                                                     >  1:21141/1‑251 (MQ=255)
                                                                                                                                                              atgaagCTGATGACGGCACCATCCCAGCCACCGTTATTGAGGTTACAGAAACAGGCGCAAGAAAGGCCAGCCACGAATGCTACCGGCTATCAGGGGGAACGGAAGGGCTTAATTGGGGTAAATCGTATCCCTACGCCTGTGTAATCCAGCTGATGATGGTCTGCAAGAATCACAATGTGCTTGACTTCAGTCACCATTTGCATGCAAAGGCCTCGATCGTGATTGTGACGTGTAGATGTCatgccttgacg                                                                                   <  2:517240/249‑11 (MQ=255)
                                                                                                                                                              aagaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTc                                                                                                                                              >  1:413696/1‑192 (MQ=255)
                                                                                                                                                              aagaagGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTc                                                                                                                                              <  2:413696/192‑1 (MQ=255)
                                                                                                                                                                           cGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGtt                                                                      >  1:391352/1‑251 (MQ=255)
                                                                                                                                                                            ggCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTg                                                                     >  1:147871/1‑251 (MQ=255)
                                                                                                                                                                                     cccAGCCACAGTTATTGGGTTTACAGAAACAGGCGCAAGAACGGCCAACCATTAAGTCTCCAAGCTATATTGGGTAACGTAATGGCTGGATTTGGCTAAATCGTTTCTCTACGCCTTTGTAAACCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGACCCCCCATGAATATTAATGCCGCAATCGTGATTTTTACGTGTCCATGTCCGGCATTGTCCACCTTTAACAGGTGTcgacact                                                             >  2:14322/1‑245 (MQ=255)
                                                                                                                                                                                         gCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATgg                                                         <  1:232530/251‑1 (MQ=255)
                                                                                                                                                                                          ccACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTAACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGc                                                        <  1:4968/251‑1 (MQ=255)
                                                                                                                                                                                              cGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTc                                                    <  2:147871/251‑1 (MQ=255)
                                                                                                                                                                                                 tATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATTCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATgg                                                         <  2:222179/243‑1 (MQ=255)
                                                                                                                                                                                                 tATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATTCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATgg                                                         >  1:222179/1‑243 (MQ=255)
                                                                                                                                                                                                 tATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGa                                                  >  1:332496/1‑250 (MQ=255)
                                                                                                                                                                                                  aTTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAag                                                <  1:391510/251‑1 (MQ=255)
                                                                                                                                                                                                    tGAGTGTACAGAACCGGGCGCAAGAAAGGCCAACCAGTAAGTGTACCAGCTATCTTTGGTAACGCAATGGCTGAATTTGGCTAAAACGTTTCTCTACGCCTTTTTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGTACTTCAGTCACCATATGCAAATTAATGCCGCGAGCGTTATTTTTACGTGTCGATGTCAGGCATTGCCCAACTTTAATAGTTGTCGGCACTATGGCGTTCGAaga                                               <  2:480180/250‑1 (MQ=255)
                                                                                                                                                                                                            ccgAAACAGGCGCAAGAAAGGCCAACCATTAATGCTACCAGCTATCTTGGGTAAAGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTTGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCAATGCCCATCTTTAATAGTTGTCAGCACTATTGCGTTCGAAGAGATAGAACt                                      <  2:49384/249‑1 (MQ=255)
                                                                                                                                                                                                                        cAAGAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTc                          >  1:397437/1‑251 (MQ=255)
                                                                                                                                                                                                                           gAAAGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTTTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCag                       <  2:155554/251‑2 (MQ=255)
                                                                                                                                                                                                                             aaGGCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTc                     >  1:508126/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCCAACCATTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTccc                   >  1:408372/1‑250 (MQ=255)
                                                                                                                                                                                                                                  cAACCATTAAGGCTACCAGCTATCTTGGGTTACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCCCCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGGCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGt                <  2:508126/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaTTAAGGCTACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGTACTTCAGTCACCACATGCATATTAATGCCGCGAGCGTGATTTTTACGTGGCGATGTCAGGCATTGCCCATCTTTATTAGTTGTCAGCCCTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTgacg            <  2:156729/250‑4 (MQ=255)
                                                                                                                                                                                                                                               tACCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTc   >  1:219819/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aCCAGCTATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGCTCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCg  >  1:337539/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
CGGCCATCGCAATGGTGGGGGTATTGCTGAAAGTGGGGAGTTGCAAAAGCAATCCGGAAATGGTGGTGGCGGCGAAAGCGGTAAGGCAAAAGTTGATCTGTGGATGAAGATGACGTTGTGCCAGCAGCTGGCGGATATACATCGCGGTCGTACTGGCAAAGAAGGTGATGACGGCACCATCCCAGCCACCGTTATTGAGTTTACAGAAACAGGCGCAAGAAAGGCCAACCATTAAGGCTACCAGCCATCTTGGGTAACGTAATGGCTGAATTTGGCTAAATCGTTTCTCTACGCCTTTGTAATCCAGCAGATGATGCTCCGCAAGAATCACAATGTGCTGGACTTCAGTCACCACATGCATATTAATGCCGCGATCGTGATTTTTACGTGTCGATGTCAGGCATTGCCCATCTTTAATAGTTGTCAGCACTATGGCGTTCGAAGAGATAGAACTTTCGACGCTGTCCATTCCCAGTGCCCGACCCAGTCG  >  NC_000913/4602264‑4602753

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: