Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,969,103 2 bp→GC coding (232‑233/1263 nt) wecC → UDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,969,1030CG99.3% 401.9 / NA 148R78A (CGC→GCC) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (79/68);  minor base T (0/1);  total (79/69)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,969,1040GC98.6% 374.1 / NA 143R78A (CGC→GCC) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base C (73/67);  minor base A (2/0);  total (75/68)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

CGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCG  >  NC_000913/3969035‑3969171
                                                                    ||                                                                   
tgAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggct                                                               >  2:553365/2‑76 (MQ=255)
ggAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggct                                                               <  2:5324035/75‑1 (MQ=255)
 gAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggctg                                                              >  2:2443621/1‑76 (MQ=255)
 gAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggctg                                                              >  1:2671335/1‑76 (MQ=255)
  ggtgtAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGa                                                             <  2:4329810/75‑1 (MQ=255)
    tgtAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                           <  1:2045513/76‑1 (MQ=255)
      tagtaAAACCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGc                                                         <  1:2794505/76‑1 (MQ=255)
      tagtaAAACCGGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGc                                                         <  1:179493/76‑1 (MQ=255)
      tagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGc                                                         <  1:5450704/76‑1 (MQ=255)
       agtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                        >  1:1375061/1‑76 (MQ=255)
       agtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                        >  2:2001420/1‑76 (MQ=255)
         taAAACCTGCCGGAGAAGGCGGTGTTTTACGAGGGAGCAGGACGCCGGTTGAAGCGGATGCCTGGCTGATTGCTGt                                                      <  2:1035919/76‑1 (MQ=255)
          aaaacCTGCCGTAGAGGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTa                                                     <  2:2157772/76‑1 (MQ=255)
            aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                   >  1:3256200/1‑76 (MQ=255)
            aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                   >  2:2288772/1‑76 (MQ=255)
             aaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAccc                                                  >  1:5470265/1‑76 (MQ=255)
              ccTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCa                                                 <  2:1315850/75‑1 (MQ=255)
              aCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCa                                                 <  1:553365/76‑1 (MQ=255)
               cTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                <  1:5451150/76‑1 (MQ=255)
               cTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                >  2:5336222/1‑76 (MQ=255)
                tGCCGTGGAGGGCGTTTTTTTACGAGGGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                               <  1:3470525/76‑1 (MQ=255)
                tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                               >  1:803433/1‑76 (MQ=255)
                tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                               >  1:4629425/1‑76 (MQ=255)
                tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGACGCCTGGCTGGCTGATGTAGCCACg                                               >  1:2357829/1‑76 (MQ=255)
                 gCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGc                                              <  2:3257876/76‑1 (MQ=255)
                 gCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGc                                              <  1:476922/76‑1 (MQ=255)
                 gCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGc                                              >  1:5575484/1‑76 (MQ=255)
                  ccGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGcc                                             >  1:3420599/1‑76 (MQ=255)
                   cGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCg                                            >  2:4143753/1‑76 (MQ=255)
                    gTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGt                                           >  2:3238406/1‑76 (MQ=255)
                      aGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGctt                                         >  2:1165153/1‑76 (MQ=255)
                       gAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTa                                        <  2:3965360/76‑1 (MQ=255)
                        aaGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGGCGTTTaa                                       >  2:338555/1‑76 (MQ=255)
                         aGGCGGTTTTTTACGAGCGAGCACGACGCCCGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAg                                      >  2:2458665/1‑76 (MQ=255)
                          ggCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                     <  2:3059087/76‑1 (MQ=255)
                          ggCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAgg                                     <  1:108980/76‑1 (MQ=255)
                           gCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                    <  2:428528/76‑1 (MQ=255)
                           gCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                    >  1:4489204/1‑76 (MQ=255)
                           gCGGTTTTTGACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCGGGCTGATTGCTGTACCCACGCCGGTGAAggg                                    <  2:471446/76‑1 (MQ=255)
                            gggTTTTTTACGAGCGGGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGc                                   <  2:3752341/75‑1 (MQ=255)
                            cGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGc                                   >  1:3886330/1‑76 (MQ=255)
                             ggTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                  >  1:3437174/1‑76 (MQ=255)
                             ggTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                  >  1:4641953/1‑76 (MQ=255)
                               ttttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGAt                                <  1:2870789/76‑1 (MQ=255)
                               ttttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGAt                                <  2:5696443/76‑1 (MQ=255)
                                tttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGCTTGCTGTACCCACGCCGTTTAATGGCGATc                               >  1:4977345/1‑76 (MQ=255)
                                tttttACGAGCGAGCACGACGCCAGCTGAAGAGGATGACCGCACGATTCCCGCACCCACGCCGTCTAAGcgcgccc                               >  2:336295/1‑73 (MQ=255)
                                 tttgACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                              <  2:568519/76‑1 (MQ=255)
                                   ttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATg                            >  1:2828138/1‑76 (MQ=255)
                                    tACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGa                           >  2:2917693/1‑76 (MQ=255)
                                    tACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGa                           <  1:1348197/76‑1 (MQ=255)
                                    gaagagcgaggAGGGCGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGa                           <  2:5579527/73‑1 (MQ=255)
                                     aCGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAg                          >  2:2709203/1‑76 (MQ=255)
                                     aCGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAg                          <  1:1111474/76‑1 (MQ=255)
                                     aCGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAg                          <  2:4236554/76‑1 (MQ=255)
                                      cgagcgagCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGc                         <  1:5336223/76‑1 (MQ=255)
                                       gagcgagcACGACGCCAGTAGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGcc                        <  2:3351301/76‑1 (MQ=255)
                                        cgcgggcACGCCGCCATTTGAAGCGTATTGCTGGCTGCTTGCTGTCCCCACGCCTTTTACGGCCGATCATGAGCCa                       <  1:4739976/75‑1 (MQ=255)
                                        agcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                       >  2:5258890/1‑76 (MQ=255)
                                        agcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                       >  2:2678208/1‑76 (MQ=255)
                                         gcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAg                      <  1:4396605/76‑1 (MQ=255)
                                         gcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAg                      >  1:652473/1‑76 (MQ=255)
                                          cgagcACGACGCCAGTTGAAGCGGATGCCTGGCCGATTGCGGTACCCACGCCGTTCAAGGGCGATCACGAGCCAGa                     >  2:1080326/1‑76 (MQ=255)
                                           gagcAGGACGCCAGTTGAAGCGGATGCCTGGCTGATTGATGGAGCCACGCCGTTTAAGGGCGCTCATGAGCGAggt                    >  1:1186429/1‑74 (MQ=255)
                                           gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                    <  1:1570429/76‑1 (MQ=255)
                                           gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                    <  1:550377/76‑1 (MQ=255)
                                           gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                    >  1:835252/1‑76 (MQ=255)
                                           gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                    >  1:1847321/1‑76 (MQ=255)
                                           gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                    <  2:540152/76‑1 (MQ=255)
                                           gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                    >  1:63599/1‑76 (MQ=255)
                                           gagcACGAAGCCAGTTGAAGCGCATGCCAGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGAGCATGAGGCAGat                    >  2:3516766/1‑76 (MQ=255)
                                            agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGata                   <  1:2001420/76‑1 (MQ=255)
                                             gcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGatat                  <  2:2562011/76‑1 (MQ=255)
                                              cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATg                 <  1:3712417/76‑1 (MQ=255)
                                              cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGGTTAAGGGACATCATGAGCCAGATATg                 >  2:1190103/1‑76 (MQ=255)
                                               acgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGa                <  1:2288772/76‑1 (MQ=255)
                                               acgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGa                >  2:1659336/1‑76 (MQ=255)
                                                cgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc               >  1:4164764/1‑76 (MQ=255)
                                                 gacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATg                 <  2:1312223/73‑1 (MQ=255)
                                                  acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCt             <  1:195279/76‑1 (MQ=255)
                                                  acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCt             <  1:1078457/76‑1 (MQ=255)
                                                   cgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTa            <  1:5776344/76‑1 (MQ=255)
                                                   cgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTa            >  1:2541797/1‑76 (MQ=255)
                                                    gCCAGTTGAAGCGGATGCCTGTCTGATTGCTTTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc           <  2:5246081/76‑1 (MQ=255)
                                                    gCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc           <  1:2315924/76‑1 (MQ=255)
                                                    gCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc           >  2:5443545/1‑76 (MQ=255)
                                                    gCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc           <  1:3024900/76‑1 (MQ=255)
                                                     ccAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACg          >  1:3128147/1‑76 (MQ=255)
                                                       aGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGtt        >  1:213747/1‑76 (MQ=255)
                                                       aGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGtt        <  2:4727195/76‑1 (MQ=255)
                                                        gTTGAAGCGGATGCCTGGCTGGTTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg       <  2:5736131/76‑1 (MQ=255)
                                                         ttGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTCCGTTgc      <  1:4092513/76‑2 (MQ=255)
                                                         ttGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGa      <  2:4402559/76‑1 (MQ=255)
                                                          tGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGaa     >  1:3809387/1‑76 (MQ=255)
                                                          tGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGaa     >  1:161422/1‑76 (MQ=255)
                                                           gAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGCTGact    >  2:583586/1‑74 (MQ=255)
                                                            aaGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATc   >  1:3221493/1‑76 (MQ=255)
                                                            aaGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATc   <  1:3378638/76‑1 (MQ=255)
                                                            aaGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATc   >  2:2122443/1‑76 (MQ=255)
                                                             aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGATCTAAGTTGAAtag  >  1:1206142/1‑74 (MQ=255)
                                                                    ||                                                                   
CGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCG  >  NC_000913/3969035‑3969171

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: