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,1030CG100.0% 324.1 / NA 121R78A (CGC→GCC) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (67/54);  total (67/54)
*NC_0009133,969,1040GC99.1% 288.4 / ‑6.2 113R78A (CGC→GCC) wecCUDP‑N‑acetyl‑D‑mannosaminuronic acid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (59/53);  minor base A (1/0);  total (60/53)
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.74e-01

Only 100 of 122 total aligned reads displayed.

GCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCG  >  NC_000913/3969034‑3969179
                                                                     ||                                                                           
gCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggc                                                                        >  2:4117272/1‑76 (MQ=255)
 ggAGTGTAGAAAAACCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggct                                                                       <  1:1374291/75‑1 (MQ=255)
  gAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggctg                                                                      >  2:3330756/1‑76 (MQ=255)
  gAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggctg                                                                      >  1:4041187/1‑76 (MQ=255)
  gAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggctg                                                                      >  1:1717027/1‑76 (MQ=255)
  gAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGcctggctg                                                                      >  1:2663188/1‑76 (MQ=255)
   aGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGa                                                                     >  2:871010/1‑76 (MQ=255)
    gtgtAGTAAAAACTGCCCTAGAAGGCGGCTTTTTACGACCGAGCCCCACGCCCGTTCCAGCGGATGCCTGGCTGAt                                                                    >  2:4191409/1‑76 (MQ=255)
     tgtAGTAAAACCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                   <  1:4794521/76‑1 (MQ=255)
     ggtagaaAAACCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                   <  1:118836/75‑1 (MQ=255)
     ggtagaaAAACCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGAtt                                                                   <  2:2781518/75‑1 (MQ=255)
       tagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGc                                                                 >  2:3428624/1‑76 (MQ=255)
       tagtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGc                                                                 <  1:2271795/76‑1 (MQ=255)
       tagtaAAAACTGCCGTAGAAGGCGGTGTTTTAAGAGCGAGCACGACGACAGTTGAAGCGGATGCCTGGCTGATTGc                                                                 >  1:1824461/1‑76 (MQ=255)
        agtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                                >  2:2892103/1‑76 (MQ=255)
        agtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCt                                                                >  1:2964209/1‑76 (MQ=255)
         gtaAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTg                                                               >  2:3769657/1‑76 (MQ=255)
          aaaaacCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGt                                                              <  1:1427065/75‑1 (MQ=255)
           aaaacCTGCCCGAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTa                                                             <  2:4722714/76‑1 (MQ=255)
             aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                           <  1:1950392/76‑1 (MQ=255)
             aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                           >  2:3413025/1‑76 (MQ=255)
             aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                           >  1:305846/1‑76 (MQ=255)
             aaaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAcc                                                           >  1:1523853/1‑76 (MQ=255)
              accTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAccc                                                          <  2:3231444/74‑1 (MQ=255)
              aaCTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTAccc                                                          >  2:2194196/1‑76 (MQ=255)
               aCTGCCGTAGAAGGCGGTTTTTTACGAGCGCGCACGACGCCCGTTGACGCGGAAGCCTGGCTGATTGCTGTACCCa                                                         <  2:3902153/76‑1 (MQ=255)
                cTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                        <  1:2048447/76‑1 (MQ=255)
                cTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCAc                                                        >  1:4614897/1‑76 (MQ=255)
                 tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                       >  2:3550335/1‑76 (MQ=255)
                 tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                       >  1:1153287/1‑76 (MQ=255)
                 tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                       <  1:3632134/76‑1 (MQ=255)
                 tGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACg                                                       >  1:1442411/1‑76 (MQ=255)
                   ccGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGcc                                                     <  1:780992/76‑1 (MQ=255)
                    cGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCg                                                    <  1:4371878/76‑1 (MQ=255)
                    cGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCg                                                    <  2:3646916/76‑1 (MQ=255)
                    cGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCg                                                    >  1:184288/1‑76 (MQ=255)
                    cGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCg                                                    <  2:3572358/76‑1 (MQ=255)
                     gTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGt                                                   >  2:2898197/1‑76 (MQ=255)
                     gTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGt                                                   >  1:1025357/1‑76 (MQ=255)
                     gTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGt                                                   >  1:1072648/1‑76 (MQ=255)
                      tAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGtt                                                  >  2:2643802/1‑76 (MQ=255)
                       aGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGttt                                                 >  1:3652187/1‑76 (MQ=255)
                       aGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGttt                                                 <  2:454777/76‑1 (MQ=255)
                         gaGGCGGTGTTTTACGAGCGAGCGCGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTaa                                               <  2:3431232/75‑1 (MQ=255)
                         aaGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTaa                                               >  1:1998876/1‑76 (MQ=255)
                          aGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAg                                              <  1:4369826/76‑1 (MQ=255)
                          aGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGCTTTAg                                              >  2:1745993/1‑76 (MQ=255)
                           ggCGGTTTTTTATGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTTATTGCTGTACCCTCGCCGTTTAAgg                                             <  2:3465652/76‑1 (MQ=255)
                            gCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                            <  1:2751846/76‑1 (MQ=255)
                            gCGGTTTTTGACGAGCGCGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAggg                                            <  1:2900240/76‑1 (MQ=255)
                             cGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGc                                           <  1:1032977/76‑1 (MQ=255)
                             cGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGc                                           >  1:816785/1‑76 (MQ=255)
                             cGGTTTTTGACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTCTAACGGc                                           >  2:280256/1‑76 (MQ=255)
                              ggTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCg                                          >  1:785236/1‑76 (MQ=255)
                                  ttttACGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCa                                      <  2:3334029/76‑1 (MQ=255)
                                     ttggagcggggAGGACGGCAGGTGGAGGGGATGCCTGGGTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGa                                   <  2:327599/73‑1 (MQ=255)
                                      aCGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAg                                  <  1:4126167/76‑1 (MQ=255)
                                      aCGAGCGAGCACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAg                                  >  2:3961549/1‑76 (MQ=255)
                                        gagcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGcc                                <  1:4438203/76‑1 (MQ=255)
                                        gagcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGcc                                >  1:4193391/1‑76 (MQ=255)
                                        gagcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGcc                                <  2:816785/76‑1 (MQ=255)
                                        gagcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGcc                                >  1:2352170/1‑76 (MQ=255)
                                        gagcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGCTTAAGGGCGATCATGAGcc                                >  2:897350/1‑76 (MQ=255)
                                         agcgagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCa                               <  1:4120590/76‑1 (MQ=255)
                                            gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                            >  2:2409803/1‑76 (MQ=255)
                                            gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                            >  1:4238743/1‑76 (MQ=255)
                                            gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                            >  2:840154/1‑76 (MQ=255)
                                            gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                            <  2:3593269/76‑1 (MQ=255)
                                            gagcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGat                            >  1:4640468/1‑76 (MQ=255)
                                            gagcACGACGCCAGCTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGAGAGATCATGAGCCAGaa                            >  2:755584/1‑75 (MQ=255)
                                             agcACGCCGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAtata                           <  1:617633/76‑1 (MQ=255)
                                             agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGata                           <  1:78173/76‑1 (MQ=255)
                                             agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGata                           <  2:4087419/76‑1 (MQ=255)
                                             agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGata                           >  2:2384231/1‑76 (MQ=255)
                                             agcACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGata                           >  2:2615513/1‑76 (MQ=255)
                                               cACGACGCGAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATg                         <  2:3552266/76‑1 (MQ=255)
                                               cACGACGCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATg                         >  2:3253029/1‑76 (MQ=255)
                                                acgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGa                        <  2:1153286/76‑1 (MQ=255)
                                                 cgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc                       <  2:3411941/76‑1 (MQ=255)
                                                 cgacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAc                       <  2:3877891/76‑1 (MQ=255)
                                                  gacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAcc                      >  1:3255053/1‑76 (MQ=255)
                                                  gacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAcc                      >  2:4347594/1‑76 (MQ=255)
                                                  gacgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGAcc                      >  2:3814893/1‑76 (MQ=255)
                                                   acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCt                     <  1:2283747/76‑1 (MQ=255)
                                                   acgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCt                     >  1:2576051/1‑76 (MQ=255)
                                                    cgCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTa                    <  1:3865297/76‑1 (MQ=255)
                                                     gCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATGATGAGCCAGATATGACCTAc                   <  2:1705841/76‑1 (MQ=255)
                                                     gCCAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTAc                   <  1:731339/76‑1 (MQ=255)
                                                      cctgTTGAAGAGGATGCCGGGCTGATTGCTGGACCCACGCCGCTTAAGGGCGATCATCAGCCAGATATCACCTACg                  >  2:1320615/4‑76 (MQ=255)
                                                       cAGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGt                 <  1:3609951/76‑1 (MQ=255)
                                                        aGTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGtt                >  2:3512407/1‑76 (MQ=255)
                                                         gTTGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTg               >  2:3487649/1‑76 (MQ=255)
                                                           tGAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGaa             >  2:3692873/1‑76 (MQ=255)
                                                            gAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAt            <  2:4376377/76‑1 (MQ=255)
                                                            gAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAt            <  2:4238743/76‑1 (MQ=255)
                                                            gAAGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAAt            >  1:4392666/1‑76 (MQ=255)
                                                             aaGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATc           >  1:3167812/1‑76 (MQ=255)
                                                             aaGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATc           >  2:1772286/1‑76 (MQ=255)
                                                              aGCGGATGCCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCg          <  1:1781850/76‑1 (MQ=255)
                                                                      gctgggtgATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTgctcg  <  1:328693/76‑1 (MQ=255)
                                                                     ||                                                                           
GCGAGTGTAGTAAAAACTGCCGTAGAAGGCGGTTTTTTACGAGCGAGCACGACGCCAGTTGAAGCGGATCGCTGGCTGATTGCTGTACCCACGCCGTTTAAGGGCGATCATGAGCCAGATATGACCTACGTTGAATCGGCTGCTCG  >  NC_000913/3969034‑3969179

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