Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,685,197 T→C 16.3% K464K (AAA→AAG fumC ← fumarate hydratase (fumarase C),aerobic Class II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,685,1970TC16.3% 173.8 / 25.6 79K464K (AAA→AAGfumCfumarate hydratase (fumarase C),aerobic Class II
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base C (7/6);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCT  >  NC_000913/1685100‑1685296
                                                                                                 |                                                                                                   
aTGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCttt                                                                                                   <  2:856034/100‑1 (MQ=255)
  gCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCa                                                                                                 <  2:828147/100‑1 (MQ=255)
   ccGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCAt                                                                                                >  2:1286016/1‑100 (MQ=255)
   ccGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCAt                                                                                                <  2:242886/100‑1 (MQ=255)
   ccGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCAt                                                                                                <  1:1286010/100‑1 (MQ=255)
     gACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATAc                                                                                              <  2:586309/100‑1 (MQ=255)
      aCAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACt                                                                                             >  1:537688/1‑100 (MQ=255)
        aaTGTGCAGCACCGTTATAAACCGCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGc                                                                                           <  1:1631569/100‑1 (MQ=255)
         aTGTGCAGCACCGTTATAAACGTCAGGCGCAGCCGCTTCGCTTAATCATTCCACGGCTGCAACTGTATGTTGCCGATTAACGCCCGGCTTTCATACTGcc                                                                                          >  1:816787/1‑100 (MQ=255)
         aTGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGcc                                                                                          <  1:218218/100‑1 (MQ=255)
          tgtgCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCg                                                                                         <  2:368521/100‑1 (MQ=255)
           gtgCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGa                                                                                        <  2:39500/100‑1 (MQ=255)
            tgCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGAc                                                                                       >  2:671145/1‑100 (MQ=255)
            tgCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGAc                                                                                       >  1:443024/1‑100 (MQ=255)
               agcaCCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCAt                                                                                    >  2:509574/1‑100 (MQ=255)
               agcaCCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCAt                                                                                    >  1:246636/1‑100 (MQ=255)
                  aCCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTg                                                                                 >  2:2383259/1‑100 (MQ=255)
                   ccGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGt                                                                                >  2:1071215/1‑100 (MQ=255)
                       tataAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTg                                                                            <  1:733002/100‑1 (MQ=255)
                       tataAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTg                                                                            >  2:1204976/1‑100 (MQ=255)
                        ataAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTgg                                                                           <  1:2288897/100‑1 (MQ=255)
                        ataAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTgg                                                                           >  2:1781498/1‑100 (MQ=255)
                         taAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTgg                                                                           >  1:1920226/1‑99 (MQ=255)
                         taAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGc                                                                          <  2:2211772/100‑1 (MQ=255)
                         taAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGc                                                                          <  1:633219/100‑1 (MQ=255)
                           aaCCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCg                                                                        <  1:1715336/100‑1 (MQ=255)
                            aCCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGt                                                                       >  1:1258852/1‑100 (MQ=255)
                            aCCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGt                                                                       <  1:1755879/100‑1 (MQ=255)
                               tCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTAcc                                                                    <  2:1912574/100‑1 (MQ=255)
                                      cAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGt                                                             <  2:1778941/100‑1 (MQ=255)
                                        gCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTc                                                            >  1:648712/1‑99 (MQ=255)
                                              tCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCg                                                     <  2:1909165/100‑1 (MQ=255)
                                                  ttGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCtt                                                 <  1:2567305/100‑1 (MQ=255)
                                                   tGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCCAACTCGGCTTc                                                >  1:722574/1‑100 (MQ=255)
                                                    gttCATTCCACGGCTGCACCTGTATTTTGCATATTAACGCCCGGCTTTCATACTTCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCgg                                                    <  1:2193182/93‑1 (MQ=255)
                                                    gATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTc                                                <  1:558107/99‑1 (MQ=255)
                                                    gATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCg                                               >  2:1335639/1‑100 (MQ=255)
                                                     aTCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGGCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGc                                              <  1:1405537/100‑1 (MQ=255)
                                                     aTCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGc                                              >  2:274571/1‑100 (MQ=255)
                                                     aTCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGc                                              <  1:1856263/100‑1 (MQ=255)
                                                      tCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCt                                             >  2:451368/1‑100 (MQ=255)
                                                      tCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCt                                             >  1:845402/1‑100 (MQ=255)
                                                      tCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCt                                             <  1:1161369/100‑1 (MQ=255)
                                                       cATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTa                                            >  2:1541685/1‑100 (MQ=255)
                                                         ttCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAg                                          <  1:2344746/100‑1 (MQ=255)
                                                         ttCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAg                                          >  1:519803/1‑100 (MQ=255)
                                                           ccACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGa                                         <  1:305801/99‑1 (MQ=255)
                                                           ccACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGAt                                        >  2:1636332/1‑100 (MQ=255)
                                                           ccACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGAt                                        >  2:988440/1‑100 (MQ=255)
                                                            cACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATa                                       >  2:1638101/1‑100 (MQ=255)
                                                              cGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATAcc                                     >  1:1045090/1‑100 (MQ=255)
                                                              cGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATAcc                                     <  2:919257/100‑1 (MQ=255)
                                                                gCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATAcccc                                   >  2:1431106/1‑100 (MQ=255)
                                                                  tGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATAcccc                                   <  1:1571464/98‑1 (MQ=255)
                                                                    cACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAgcg                               <  1:1517491/100‑1 (MQ=255)
                                                                          tATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAggg                         <  1:294735/100‑1 (MQ=255)
                                                                          tATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAggg                         <  1:1675715/100‑1 (MQ=255)
                                                                             gTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCg                      >  1:206581/1‑100 (MQ=255)
                                                                              ttGCAGATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGc                     <  1:1647488/100‑1 (MQ=255)
                                                                               tGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCa                    <  1:2254455/100‑1 (MQ=255)
                                                                                 cAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGc                  >  1:927904/1‑100 (MQ=255)
                                                                                 cAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGc                  >  2:2018882/1‑100 (MQ=255)
                                                                                 cAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGc                  >  2:2501537/1‑100 (MQ=255)
                                                                                  aGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCt                 <  2:479240/100‑1 (MQ=255)
                                                                                  aGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCt                 >  1:2139789/1‑100 (MQ=255)
                                                                                   gATTAACGCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCtt                >  2:1245484/1‑100 (MQ=255)
                                                                                    aTTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCttt               <  2:537690/100‑1 (MQ=255)
                                                                                    aTTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCttt               <  2:782022/100‑1 (MQ=255)
                                                                                     ttAACGCCCGGCCTTCGTACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCtttt              >  2:777657/1‑100 (MQ=255)
                                                                                      tAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTa             <  1:2304377/100‑1 (MQ=255)
                                                                                       aaCGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTaa            >  1:2153587/1‑100 (MQ=255)
                                                                                       aaCGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTaa            <  1:2501518/100‑1 (MQ=255)
                                                                                         cGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAgg          <  1:154881/100‑1 (MQ=255)
                                                                                         cGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAgg          >  2:1027598/1‑100 (MQ=255)
                                                                                         cGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGCGCCGCAGCTTTTAAgg          >  1:1511101/1‑100 (MQ=255)
                                                                                          gCCCGGCTTTCGTACTGCCGACCTTCTGTTTGTGCTTTTCCCAGCTTTCAAACTCGGCTTCGCTAAGTTACCCCAGCGCCAGGTCTGCATCtggtagtg          <  1:1822587/99‑8 (MQ=255)
                                                                                          gCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGt         >  1:768260/1‑100 (MQ=255)
                                                                                          gCCCGGCCTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGt         >  1:157698/1‑100 (MQ=255)
                                                                                              ggCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGc     >  2:865742/1‑100 (MQ=255)
                                                                                                cTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGccc   >  2:278039/1‑100 (MQ=255)
                                                                                                 tttCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCct  <  2:510877/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
ATGCCGACAATGTGCAGCACCGTTATAAACCTCAGGCGCAGCCGCTTCGTTTGATCATTCCACGGCTGCACCTGTATGTTGCAGATTAACGCCCGGCTTTCATACTGCCGACCATCTGTTCTGGCCGTACCCAGCTGTCAAACTCGGCTTCGCTAAGATACCCCAGCGCAAGGGCCGCAGCTTTTAAGGTCAGCCCT  >  NC_000913/1685100‑1685296

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