Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9010TA41.1% 69.3 / 26.9 56Y44F (TAT→TTT) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (24/9);  new base A (0/23);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 9.18e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGG  >  NC_000913/4336756‑4337050
                                                                                                                                                  |                                                                                                                                                      
gagGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATg                                                                                                                                                   >  1:689805/1‑151 (MQ=255)
gagGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATg                                                                                                                                                   <  2:717954/151‑1 (MQ=255)
   gACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGc                                                                                                                                                >  1:669635/1‑150 (MQ=255)
     cGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAAt                                                                                                                                             >  1:902783/1‑151 (MQ=255)
               aTAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATACATAAATGGCAATTCCGcc                                                                                                                                       >  2:796884/1‑147 (MQ=255)
                            gggCCAAGGCAGCGGCGGGTATAAGCGTAAAAACCACCGGGCTTTGGGTCGAGAAACTACTTTTGGGATTTTCCCACCGAGCCCCCGAGTGCACTGATAATCGTCACCAACCACCCAAAAGGGGCAATCCCGCCAGTAGAGGCCAGGTTTg                                                                                                                      <  2:595435/151‑1 (MQ=255)
                                cAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGACCCAGTGCACCGATAATCGTCACCAAACATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAgg                                                                                                                  >  2:142852/1‑151 (MQ=255)
                                 aaagagggggcgggggAAAGCGAAAGACCCCCCAGGATTGGGGGCGAGAAACAACATTTGGGCGTATCCCACCGAGAGCCCCAGTGCACCGAAAATCGTCCCCACCCATCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAg                                                                                                                   <  2:791416/141‑1 (MQ=255)
                                  aggcagcGGGGGGCGTAGGCATAAAACCCACCAGGCTTGGGTTCGAGGACCGACTTTTTGGCTTATCCCACCGAGACCCCCAGTGCACCGATACTCGTCCCCACCCATCAATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTa                                                                                                                <  2:399371/149‑1 (MQ=255)
                                     cagcGGGGGGCGTAAGCATAAGAACCCCCAGGACTTGGGTCGAGGAACGACATTTGGGCGTATCCCATCGAGACCCCCAGTGCCCCGATAATCGTCACCACCCATCCAAAAAGGAAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTaa                                                                                                               <  1:935775/149‑1 (MQ=255)
                                                 taagcgtaagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                                                                                 >  1:1059663/1‑151 (MQ=255)
                                                  aagctaaagACCCCCCAGGATTGGGGTGGAGGAACGAAATTTTGGCGTATACCAGCGAGACCCCCAGTGCACAGATATTCGCCACAAACCACCAAAAAATGACAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                                                                                 <  2:505172/144‑1 (MQ=255)
                                                  aagcgtaagACCCACCAGGATTGGGGTCGAGGAACGACATTTGGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCCCCACCCATCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTg                                                                                                 <  2:365286/150‑1 (MQ=255)
                                                     ctaaaaacccaccgGGACTTGGGGCGAGAAACTACTTTTGGGCGTATACCATAGAGAGCCCCAGTGCCCCTATATCCGTCCCAACCCATCAAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAcc                                                                                              <  2:372037/142‑1 (MQ=255)
                                                         agACCCACCAGGACTTGGGTCGAGGAACGACTTTTTGGCGTATCCCATCGAGACCCCCAGTCCACCGATAATCGTCACCACCCACCAAAAAATGACAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCata                                                                                         <  2:322662/151‑1 (MQ=255)
                                                            accaccAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataat                                                                                       >  2:736688/1‑150 (MQ=255)
                                                                  aGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTcccc                                                                                >  1:370291/1‑151 (MQ=255)
                                                                  aGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGAAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTcccc                                                                                >  1:1020604/1‑151 (MQ=255)
                                                                   ggACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCCCCGATAATCGTCACCAACCACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCg                                                                               <  1:541279/151‑1 (MQ=255)
                                                                   ggACTGGGGTCGAGGAACGACATTTTGGCGTATCCCATCGAGGCCCCCAGTGCCCCGATAACCGTCACAAACCACCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCg                                                                               <  1:712627/151‑1 (MQ=255)
                                                                         gggTCGAGGAACGACATTTGGGCGTATACCAGCGAGACCCCCAGTGCACCGATAATCGCCCCCAACCATCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCa                                                                         <  2:999418/151‑1 (MQ=255)
                                                                             tCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATc                                                                       >  2:441601/1‑150 (MQ=255)
                                                                                 ggAACGACATTTTGGCGTATCCCATCGAGACCCCCGGTCCACCGATAACCGCCCCAAACCATCCAAAAATGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGt                                                                  <  2:23953/151‑1 (MQ=255)
                                                                                   aaCGACATTTTGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCTTCACCAACCACCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGa                                                                <  1:1119249/151‑1 (MQ=255)
                                                                                       aCATTTTGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCGCCACCACCCACCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACggg                                                            <  1:976951/151‑1 (MQ=255)
                                                                                          tttgGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCACCCATCCATAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGa                                                          <  1:1035564/150‑1 (MQ=255)
                                                                                           tttGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAtt                                                        >  2:729825/1‑151 (MQ=255)
                                                                                           tttGGCGAATACCAACGAGCGCCCCAGTGCACCGATAACCGCCACCACCCACCCATAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                         <  1:502903/150‑1 (MQ=255)
                                                                                            ttGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCGTCCCCAACCACCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAtt                                                        <  1:872405/150‑1 (MQ=255)
                                                                                              ggCGTATCCCACCGAGACCCCCAGTGCCCCGATAATCGTCACCACCCACCCATAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTaa                                                      <  1:767906/150‑1 (MQ=255)
                                                                                                 gTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGccc                                                  >  1:113369/1‑151 (MQ=255)
                                                                                                  taacctaTCGAGCGCCCAGGGCCACAGATAGGCGCCGCCACCCATCAAAAAACGCCACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGcc                                                   <  2:349236/143‑1 (MQ=255)
                                                                                                   ataCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCAc                                                >  1:1063061/1‑151 (MQ=255)
                                                                                                   ataCCAGGGAGAGCCCCAGGGCACCGATAATCGTCACCACCCACCCATAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGccc                                                  <  2:264302/149‑1 (MQ=255)
                                                                                                      ccATCGAGACCCCCAGTGCACCGATAATCGTCACCACCCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACt                                               <  1:829173/149‑1 (MQ=255)
                                                                                                        aTCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgt                                           >  1:644988/1‑151 (MQ=255)
                                                                                                          cGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                          >  1:1115064/1‑150 (MQ=255)
                                                                                                           gagaGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                          >  2:1039716/1‑149 (MQ=255)
                                                                                                            agagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTtgtg                                          <  1:620052/148‑1 (MQ=255)
                                                                                                                ccccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                    >  1:451025/1‑150 (MQ=255)
                                                                                                                ccccAGTGCACAGATAATCGTCACCAACCATCCATAAAGGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcat                                     <  1:529113/149‑1 (MQ=255)
                                                                                                                   cAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatcagc                                 >  2:487185/1‑150 (MQ=255)
                                                                                                                          ccGATAATCGCCACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                         <  1:51378/151‑1 (MQ=255)
                                                                                                                            gATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGAc                       >  2:479454/1‑151 (MQ=255)
                                                                                                                             aTAATCGTCCCCAACCATCCATAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGa                           <  1:459734/146‑1 (MQ=255)
                                                                                                                             aTAACCGCCACCACCCACCAAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGAc                       <  2:513848/150‑1 (MQ=255)
                                                                                                                               aaccgcccccaaccaCCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                    <  2:508192/143‑1 (MQ=255)
                                                                                                                                aTCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                    >  1:1056446/1‑150 (MQ=255)
                                                                                                                                  cgccACCAACACCCAAAAAAGGCCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAt                  <  2:217092/147‑1 (MQ=255)
                                                                                                                                  cGTCACAACCGATCAAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAt                  <  2:800593/150‑1 (MQ=255)
                                                                                                                                    tCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                >  2:442190/1‑150 (MQ=255)
                                                                                                                                     caccacCCACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaaa              <  1:729164/151‑1 (MQ=255)
                                                                                                                                     cACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTaa               >  1:1087049/1‑150 (MQ=255)
                                                                                                                                      accaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAAc             >  1:868815/1‑151 (MQ=255)
                                                                                                                                       ccaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACt            >  1:153841/1‑151 (MQ=255)
                                                                                                                                         aaccaTCAAAAAAGGGCATTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAt                  <  2:689822/143‑1 (MQ=255)
                                                                                                                                           ccatccatAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTg        >  1:22046/1‑151 (MQ=255)
                                                                                                                                            catcaaaAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGc       <  1:54326/151‑1 (MQ=255)
                                                                                                                                                 atAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAgg  >  2:13172/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGG  >  NC_000913/4336756‑4337050

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