Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,744,424 A→G 11.6% L83S (TTA→TCA)  apt ← adenine phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,744,4240AG11.6% 137.0 / 9.8 59L83S (TTA→TCA) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base G (5/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

TCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTT  >  CP000730/1744289‑1744564
                                                                                                                                       |                                                                                                                                            
tCAATCGTACCACCAGTAGCGAATAAATCATCCGTCCCTTACACGCGTTGACGTACCTTAATTGCATATTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATCACTTCACGCGGTAATTTc                                                                                                                                         >  1:181292/1‑141 (MQ=255)
 cAATCGTACCACCATTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                                                                                                                                        <  1:53978/141‑1 (MQ=255)
 cAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTTTCGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                                                                                                                                        <  1:369953/141‑1 (MQ=255)
 cAACCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTcc                                                                                                                                        <  1:54044/141‑1 (MQ=255)
    tCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAAGTAATTTCCctt                                                                                                                                     >  1:660136/1‑141 (MQ=255)
          caccaGGAGCTAATAAAACATCTGTAATGAACACACGTTTACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTc                                                                                                                               <  2:452679/141‑1 (MQ=255)
           accaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCAGGCGGTGATTTCCCTTCTTTTCt                                                                                                                              >  2:392205/1‑141 (MQ=255)
              aGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAc                                                                                                                           <  2:355433/141‑1 (MQ=255)
               gTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTTATAACGAATGACTTCACGAGGTAATTTTCCTTCTTTTTTAAct                                                                                                                          >  1:322574/1‑140 (MQ=255)
                        aaaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAg                                                                                                                 <  2:185503/141‑1 (MQ=255)
                         aaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGc                                                                                                                >  2:646447/1‑141 (MQ=255)
                             atcTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAAt                                                                                                            >  2:707571/1‑141 (MQ=255)
                                 ttAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATcccc                                                                                                        <  2:423567/140‑1 (MQ=255)
                                 gTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATcccc                                                                                                        <  1:209562/141‑1 (MQ=255)
                                  tAATTACCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCa                                                                                                       <  1:343746/141‑1 (MQ=255)
                                  tAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCa                                                                                                       <  1:762557/141‑1 (MQ=255)
                                      tAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGa                                                                                                   <  1:38570/141‑1 (MQ=255)
                                      tAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGa                                                                                                   >  1:628580/1‑141 (MQ=255)
                                        acacacGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAAt                                                                                                 >  2:209615/1‑141 (MQ=255)
                                            acGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGc                                                                                             <  1:401790/141‑1 (MQ=255)
                                            acGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGc                                                                                             <  1:294022/141‑1 (MQ=255)
                                                tGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCACTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACa                                                                                         >  1:412464/1‑141 (MQ=255)
                                                tGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAAGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACa                                                                                         <  2:660136/141‑1 (MQ=255)
                                                      ggTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATAACTTCACGAGGTAATTTCCCTTCTTTTCTACCAGGTGCAAAGCCAATCCCCATTGAATCAGCCACAGGACac                                                                                   >  2:577620/1‑140 (MQ=255)
                                                         ttAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATACTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCCTTGAATAAGCTACAGGACAGCCa                                                                                >  2:3176/1‑141 (MQ=255)
                                                             ttGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACGGGACAGCCAATGa                                                                            >  1:244787/1‑141 (MQ=255)
                                                              tGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGAt                                                                           <  2:182616/141‑1 (MQ=255)
                                                                cATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaa                                                                         >  1:645982/1‑141 (MQ=255)
                                                                   cTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGc                                                                      <  2:47342/141‑1 (MQ=255)
                                                                   cTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGCCAGCCAATAATAAAGc                                                                      >  2:236761/1‑141 (MQ=255)
                                                                     ttGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCa                                                                    <  2:446324/141‑1 (MQ=255)
                                                                      tgtgCATTGTTAACACATTTGTACCATACTCTAGGTCATACTCATAACGGATGACTTCACGAGTTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGGTAAAGcacc                                                                   >  2:529979/1‑138 (MQ=255)
                                                                        tacATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcgc                                                                 >  1:349824/3‑141 (MQ=255)
                                                                        tacATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATACCCATTGAATAAGCTACAGGACAGCCAATGATAAAACCAcgc                                                                 >  2:687850/3‑141 (MQ=255)
                                                                         gCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcgcg                                                                <  1:667276/141‑1 (MQ=255)
                                                                          cATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACgcgc                                                               >  2:422925/1‑141 (MQ=255)
                                                                              gTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAACGCCACCCGCTTCa                                                           >  2:481786/1‑141 (MQ=255)
                                                                                tAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCGCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTAAAGGACAGACACTGATAAAGCTCCGCGCTCCAgg                                                         >  1:643957/1‑141 (MQ=255)
                                                                                tAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAgg                                                         <  2:758861/141‑1 (MQ=255)
                                                                                 aaaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGt                                                        <  1:254398/141‑1 (MQ=255)
                                                                                        ttGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGCCAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACaa                                                 >  1:133135/1‑141 (MQ=255)
                                                                                         tGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAAc                                                >  1:410874/1‑141 (MQ=255)
                                                                                         tGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAAc                                                >  2:560047/1‑141 (MQ=255)
                                                                                            cccATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGat                                             <  1:356389/140‑1 (MQ=255)
                                                                                               atatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATc                                          <  1:445282/141‑1 (MQ=255)
                                                                                                 atTCTAGGCCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCCTTGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCaa                                        <  2:315026/141‑1 (MQ=255)
                                                                                                   tCTAGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAAAGATAAAGCAAAGCGCTTCAGCTCCTACACCgcgacacaaa                                      >  2:489003/1‑132 (MQ=255)
                                                                                                      aGGTCATACTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtct                                   >  2:22749/1‑141 (MQ=255)
                                                                                                      aGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCACTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtcc                                   >  2:214877/1‑140 (MQ=255)
                                                                                                         tCATCCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTg                                <  2:677238/141‑1 (MQ=255)
                                                                                                             aCTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCtt                            >  2:142451/1‑141 (MQ=255)
                                                                                                              cTCATAACGAATGACTTCACGAGGTGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCttt                           <  1:646447/141‑1 (MQ=255)
                                                                                                                    acaaaTGACTTCACGATGTAATTTCCCCTCTTTTCTAACAGGTGCAAATCCAATCCCCATTGAATAAGCTACAGTACAGCCAATGATAAACCCACTCGCTTCACGTCCTACAACGATATCAACATCTCTGTCTTTTGCTTa                     <  2:619722/138‑1 (MQ=255)
                                                                                                                    aCGAATGGCTTCCCGCTGTAATTGCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTa                     <  2:619705/141‑1 (MQ=255)
                                                                                                                        aTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                 >  1:437700/1‑141 (MQ=255)
                                                                                                                               cacaataaaaTTACCCATTTTATCTAACTGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAttt          <  2:526242/133‑1 (MQ=255)
                                                                                                                                cctaGGTAATTTCACTTCTTTTCTACCCGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtttt         <  2:744597/138‑1 (MQ=255)
                                                                                                                                aCGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAAtttt         <  1:209615/141‑1 (MQ=255)
                                                                                                                                 cGAGGTAATTTCCCTTCTTTTCTCCCATGACCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTa        >  1:757760/1‑141 (MQ=255)
                                                                                                                                   aGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATc      >  1:691305/1‑141 (MQ=255)
                                                                                                                                       aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGGATAAGCTACAGGACACCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGtt  >  2:576426/1‑141 (MQ=255)
                                                                                                                                       aaTTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGtt  <  1:116649/141‑1 (MQ=255)
                                                                                                                                       aaTTTCACTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGtt  <  1:542219/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
TCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTT  >  CP000730/1744289‑1744564

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