Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,150,316 G→A E219K (GAG→AAG)  yghA → NADP(+)‑dependent aldehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,150,3160GA100.0% 242.2 / NA 75E219K (GAG→AAG) yghANADP(+)‑dependent aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/38);  total (37/38)

GGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGT  >  NC_000913/3150078‑3150559
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
ggTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGAGTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGATGAAAGGTATTc                                                                                                                                                                                                                                         <  2:245343/251‑1 (MQ=255)
   tGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAACGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGCCTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACCGGTCGCGAGGAAAGGGATTCggg                                                                                                                                                                                                                                      >  2:13211/1‑251 (MQ=255)
      cATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCACCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCCTTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGAATTCGGGTg                                                                                                                                                                                                                                    >  1:532865/1‑250 (MQ=255)
         tCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAat                                                                                                                                                                                                                                 >  1:445922/1‑250 (MQ=255)
         gccGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCGAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGAGTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAat                                                                                                                                                                                                                                 <  2:68954/249‑1 (MQ=255)
          ccGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCCTTAACGTTTTCGCGCTGTTCTGGATAACCCAGGAAGCGATCCCACTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACAAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAAATACAGCCGAGGCTTGGCAAACCAGGTCGCGAAGAAAGGGATTCGGGAGAatat                                                                                                                                                                                                                               >  2:413143/1‑251 (MQ=255)
            ggATATTGCAGACCTCACCAGCGAACCGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTg                                                                                                                                                                                                                             >  2:388560/1‑251 (MQ=255)
            ggATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTg                                                                                                                                                                                                                             >  1:168264/1‑251 (MQ=255)
              atatTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTc                                                                                                                                                                                                                           >  1:188870/1‑251 (MQ=255)
                   gCAGCCCTCACCAGCGAACAGTTTCAAAAGACCGTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGcc                                                                                                                                                                                                                      <  2:474494/251‑1 (MQ=255)
                   gCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGcc                                                                                                                                                                                                                      <  1:195363/251‑1 (MQ=255)
                            aCCAGCGAACAGTTTAAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGACAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAACGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGAt                                                                                                                                                                                                             >  2:297668/1‑251 (MQ=255)
                                    aCAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTACCCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGAc                                                                                                                                                                                                       >  1:342739/1‑249 (MQ=255)
                                    aCAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACGTCGTCAAGCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAg                                                                                                                                                                                                     <  2:222524/251‑1 (MQ=255)
                                      aGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGc                                                                                                                                                                                                    >  1:474576/1‑250 (MQ=255)
                                          tCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGc                                                                                                                                                                                               >  2:142685/1‑251 (MQ=255)
                                            aaaaGACCTTTGCCATTATCGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATAATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTACGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCa                                                                                                                                                                                              >  1:361972/1‑250 (MQ=255)
                                             aaaGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCaa                                                                                                                                                                                             >  2:549807/1‑250 (MQ=255)
                                                  ccTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTTGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGATAGGTATTCGGGGGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTcc                                                                                                                                                                                       <  2:532865/251‑1 (MQ=255)
                                                     gtGCCATTAACGTTGTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCTCTTCGGCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcgg                                                                                                                                                                                     <  2:573239/249‑1 (MQ=255)
                                                      tGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcggcg                                                                                                                                                                                   >  2:445367/1‑251 (MQ=255)
                                                       gCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCcggcgg                                                                                                                                                                                  <  2:62306/251‑1 (MQ=255)
                                                              aCGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACg                                                                                                                                                                           >  2:489776/1‑251 (MQ=255)
                                                                           ggtCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGAt                                                                                                                                                                <  1:269431/247‑1 (MQ=255)
                                                                               tGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGc                                                                                                                                                           >  2:558189/1‑250 (MQ=255)
                                                                                ggCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAAc                                                                                                                                                                            <  2:50347/232‑1 (MQ=255)
                                                                                ggCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAAc                                                                                                                                                                            >  1:50347/1‑232 (MQ=255)
                                                                                 gCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTAACATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACGATGCGGCTACGAAGGCGTCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATGTCCGGCGGACAAACGCAGGATAAGATCCCGCAGt                                                                                                                                                        <  2:517858/251‑1 (MQ=255)
                                                                                       ccAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTGCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTc                                                                                                                                                  <  2:120646/251‑1 (MQ=255)
                                                                                          ggAAGCGATCCCCCTGCTACCGAAAGGTGCAAGGATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACGACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGc                                                                                                                                               <  2:484066/251‑1 (MQ=255)
                                                                                           gAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGCTCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCCGGATAAGATCCCGCAGTTTGGTCAGc                                                                                                                                               >  2:444991/1‑250 (MQ=255)
                                                                                                gatcccaatGCTACCGAAAGGTGCAAGTCTCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGACGCGAGGGAAGGTATTCGGGAGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGCCAAACTCAGGACAAGATCCCGCAGTTTGGTCAGCCAACGc                                                                                                                                         >  2:388644/9‑251 (MQ=255)
                                                                                                gATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGc                                                                                                                                         >  1:62325/1‑251 (MQ=255)
                                                                                                   cccccTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCg                                                                                                                                       <  2:198979/250‑1 (MQ=255)
                                                                                                     cccTGCTACCTAAAGGTGCAATTATCATCACCACTTCGTCAATCCAGGCATTCCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGAta                                                                                                                                    <  2:8205/251‑2 (MQ=255)
                                                                                                               gAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCgg                                                                                                                          <  1:549807/251‑1 (MQ=255)
                                                                                                                aaaGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCggg                                                                                                                         >  1:251417/1‑251 (MQ=255)
                                                                                                                         aGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGACCTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGGCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGTCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTAAGCAAACGCCGGTGCAACGTGAGGGGCAACcggc                                                                                                                >  2:210438/1‑251 (MQ=255)
                                                                                                                              catcaCCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGaa                                                                                                            <  1:46880/250‑1 (MQ=255)
                                                                                                                              catcaCCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGaa                                                                                                            <  1:445367/250‑1 (MQ=255)
                                                                                                                                   ccacTTCGTCAAGCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTGTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGTCGGAACTGGcc                                                                                                      <  2:240582/251‑1 (MQ=255)
                                                                                                                                         cGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGt                                                                                                 <  1:444991/250‑1 (MQ=255)
                                                                                                                                             aaTCCAGGCATCCCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGCGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCCGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAAATGGCCCCTGTATATg                                                                                            >  2:483811/1‑251 (MQ=255)
                                                                                                                                                ccAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGttt                                                                                         >  1:454586/1‑251 (MQ=255)
                                                                                                                                                  aGGCATGCCAGCCAATTCCGCATTTACTTGACTATTCGGCGACGAAGGCGGCGATTCTGAACTCCAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATAGGGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCGGTGTGGTCAGCAAACGCCGATGAAACGTGCTGGGCAACCCGCGTAACTGGTTGCTGGATATGTTTAt                                                                                       <  2:361972/251‑1 (MQ=255)
                                                                                                                                                     cATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTg                                                                                    <  1:80100/251‑1 (MQ=255)
                                                                                                                                                      aTACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTgg                                                                                   <  2:144875/251‑1 (MQ=255)
                                                                                                                                                      aTACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTgg                                                                                   >  1:325459/1‑251 (MQ=255)
                                                                                                                                                            gCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTc                                                                             <  2:342739/251‑1 (MQ=255)
                                                                                                                                                                 gTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGGTGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAgg                                                                          >  1:405234/1‑249 (MQ=255)
                                                                                                                                                                    cGCATTTACTGGACTATGCGGCGACTATGGCGGCGATTCTGAACTACAGCCGTGTCTTGGCTAAACAGTTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAATCGCAGGATAAGATCCCGCTGTTTGGTCAGCAAACGCCGATTAAACGTGCGGGTCAAGCGGCGTAACTGGCCCCTGTATAGGTGTATCTGGCAAGTCAGTAGTCg                                                                     <  2:352391/251‑1 (MQ=255)
                                                                                                                                                                      cATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCACGTCAGGAGTCg                                                                     >  1:199369/1‑249 (MQ=255)
                                                                                                                                                                      cATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAg                                                                   <  1:90139/251‑1 (MQ=255)
                                                                                                                                                                           aCTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAACCGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGCGCCGAGCTACg                                                              >  2:517462/1‑251 (MQ=255)
                                                                                                                                                                             tGGACTATTCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGGCGCGAAGAAAGGAATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCCAAATTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTAGTCCGCAAACGACGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCACTGAATATGTTTATCTGGCAAGGCAGGAGTCGAGCTACGTc                                                            >  2:479289/1‑251 (MQ=255)
                                                                                                                                                                              ggACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACg                                                              <  2:199369/248‑1 (MQ=255)
                                                                                                                                                                                aCTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCAcc                                                         <  1:297668/251‑1 (MQ=255)
                                                                                                                                                                                aCTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCAcc                                                         <  1:575926/251‑1 (MQ=255)
                                                                                                                                                                                       ggCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGa                                                    >  1:228699/1‑249 (MQ=255)
                                                                                                                                                                                       ggCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGa                                                    >  1:279989/1‑249 (MQ=255)
                                                                                                                                                                                       ggCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAg                                                  <  2:437854/251‑1 (MQ=255)
                                                                                                                                                                                        gCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGt                                                 <  1:558189/251‑1 (MQ=255)
                                                                                                                                                                                                 gcggcgATGCGGAACGCCAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATGCGGGTGAATATTGTCGCGCCAGGCCCTATCGGGACAGCACTGCAAATTGCCGGCGGACAAACGCAGGATAAGAGCCCGGAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGTCGTAACTGGCCCCTGTATATGTGTATCTGGGAAGTCAGTAGTCGAGCTACGTCACCGCAGAAGTGCACGGCgt                                        <  1:522676/251‑1 (MQ=255)
                                                                                                                                                                                                   ggcgATTCTGAACTACAGCCGAGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGACGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGACACCGCAGAAGAGCACGGCgtgt                                      >  1:393121/1‑251 (MQ=255)
                                                                                                                                                                                                     cgATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTgc                                    >  1:166081/1‑251 (MQ=255)
                                                                                                                                                                                                     cgATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTgc                                    >  1:242291/1‑251 (MQ=255)
                                                                                                                                                                                                            gAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGTATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTgcggcggcg                             <  1:388644/251‑1 (MQ=255)
                                                                                                                                                                                                               atACAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGc                          <  2:168264/250‑1 (MQ=255)
                                                                                                                                                                                                                 aCAGCCGTGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGTGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATTTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCa                         <  1:413143/250‑1 (MQ=255)
                                                                                                                                                                                                                        tGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGt                                                                                                 <  2:490606/171‑1 (MQ=255)
                                                                                                                                                                                                                        tGGCTTGGCAAAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGt                                                                                                 >  1:490606/1‑171 (MQ=255)
                                                                                                                                                                                                                            ttGGCACAACAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCCCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGCCAAACGCAGGCTAAGATCCCGCAGTTTGGTCAGCAAAAGCCGATGAACCGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGCTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCAAGGCGTGTGCGGCGGCGAGCATTTAGGTTaaa              >  2:278035/1‑250 (MQ=255)
                                                                                                                                                                                                                                   aaCAGGTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATATGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGATTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGCGCATTTAGGTTAAAAAAATGcc      >  1:219528/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cAGGTCGCGAATAAAGGTAGTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCTGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCATCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAATAAAATGCCCg    <  2:448315/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggTCGCGAAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGAAGAAGTGCACGGCGTGTGCGTCGGCGAGCATTTAGGTTAAAAAAATGCCCGGt  <  1:186322/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
GGTTGCCATTCCGGATATTGCAGACCTCACCAGCGAACAGTTTCAAAAGACCTTTGCCATTAACGTTTTCGCGCTGTTCTGGCTAACCCAGGAAGCGATCCCCCTGCTACCGAAAGGTGCAAGTATCATCACCACTTCGTCAATCCAGGCATACCAGCCAAGTCCGCATTTACTGGACTATGCGGCTACGAAGGCGGCGATTCTGAACTACAGCCGTGGCTTGGCAAAACAGGTCGCGGAGAAAGGTATTCGGGTGAATATTGTCGCGCCAGGCCCGATCTGGACAGCACTGCAAATTTCCGGCGGACAAACGCAGGATAAGATCCCGCAGTTTGGTCAGCAAACGCCGATGAAACGTGCGGGGCAACCGGCGGAACTGGCCCCTGTATATGTTTATCTGGCAAGTCAGGAGTCGAGCTACGTCACCGCAGAAGTGCACGGCGTGTGCGGCGGCGAGCATTTAGGTTAAAAAAATGCCCGGT  >  NC_000913/3150078‑3150559

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