Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,753,888 G→T 10.9% A72S (GCG→TCG)  czcA‑II → cation efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,753,8880GT10.9% 146.4 / 9.7 55A72S (GCG→TCG) czcA‑IIcation efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base T (3/3);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

GTGGGTATCCACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTGGCCCGCGAGCAAT  >  NC_002947/2753750‑2754035
                                                                                                                                          |                                                                                                                                                   
gTGGGTATCCACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGCGGTCTCGAtggcaaa                                                                                                                                           <  5:319837/149‑4 (MQ=38)
gTGGGTATCCACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGAtggctgg                                                                                                                                           <  3:90985/149‑1 (MQ=255)
aagggTATCCACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGCGGTCTCGAtggcaaa                                                                                                                                           <  3:183816/147‑4 (MQ=37)
  gggTATCCACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGcc                                                                                                                                         >  7:56084/1‑149 (MQ=255)
   ggTATCCACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCt                                                                                                                                        >  2:117075/1‑149 (MQ=255)
     tATCCACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCTGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGc                                                                                                                                      >  5:328855/1‑149 (MQ=255)
        ccACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAg                                                                                                                                   >  3:192891/1‑149 (MQ=255)
          acaGCTACCAGAAACTGCCGATCGATGCTGTGCCTGACATCACTAACTTCCAGGTGCAGATCAACACCTCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAggg                                                                                                                                 <  8:324198/149‑1 (MQ=255)
              cTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCACCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCa                                                                                                                             <  1:80994/149‑1 (MQ=255)
                 ccAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAg                                                                                                                           >  2:236339/1‑148 (MQ=255)
                  cAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAg                                                                                                                           >  6:238925/1‑147 (MQ=255)
                  cAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATAAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAg                                                                                                                           >  1:122597/1‑147 (MQ=255)
                            cgatcgatGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAggg                                                                                                                                 >  8:52926/1‑131 (MQ=255)
                            cgatcgatGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAggg                                                                                                                                 <  7:52926/131‑1 (MQ=255)
                              atcgatGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCAAGGGCTCAAGCAGACCcgctcgct                                                                                                             >  2:97864/1‑149 (MQ=255)
                                   tGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTcgc                                                                                                        >  5:146833/1‑149 (MQ=255)
                                      tgtgCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGACATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCt                                                                                                     >  5:167515/1‑149 (MQ=255)
                                       gtgCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTc                                                                                                    >  1:50824/1‑149 (MQ=255)
                                          cctgACATCACTAACGTCCATGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCggg                                                                                                 <  5:278881/146‑1 (MQ=255)
                                          ccGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCgg                                                                                                  <  4:184696/148‑1 (MQ=255)
                                            ggACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGCCCGAGCAGCGCATCGCCTTTGCCGTCGCGCCGGCGATGGCTGGCCGGCCGGGGCGGGAGCGGTCCCGCTCGCGGTCGGGCGCGGGcg                                                                                               >  4:278214/1‑148 (MQ=255)
                                            ggACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGcc                                                                                               >  2:5169/1‑149 (MQ=255)
                                                 tCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTgctt                                                                                          <  2:200815/149‑4 (MQ=255)
                                                  cACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCg                                                                                          >  5:348568/1‑148 (MQ=255)
                                                  cACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCg                                                                                          >  1:113942/1‑148 (MQ=255)
                                                   aCTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGc                                                                                                                                      <  2:69794/103‑1 (MQ=255)
                                                   aCTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGc                                                                                                                                      >  1:69794/1‑103 (MQ=255)
                                                    cTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCaa                                                                                       <  8:56084/149‑1 (MQ=255)
                                                    cTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCaa                                                                                       <  4:192891/149‑1 (MQ=255)
                                                       aCGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTc                                                                                    <  1:163596/149‑1 (MQ=255)
                                                       aCGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTc                                                                                    <  5:238925/149‑1 (MQ=255)
                                                                       cAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCGAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCgat                                                                     >  4:32086/1‑148 (MQ=255)
                                                                             cgccgcCCCTGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCa                                                              >  3:8152/1‑149 (MQ=255)
                                                                              gccgccCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCAc                                                             <  3:278214/149‑1 (MQ=255)
                                                                              gccgccCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCAc                                                             <  1:5169/149‑1 (MQ=255)
                                                                                        gCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATAtcttc                                                   <  6:348568/149‑1 (MQ=255)
                                                                                         cTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGGCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCTCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATAtcttct                                                  <  6:349861/149‑1 (MQ=255)
                                                                                               gccgcTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCGCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTcgccc                                            <  8:341182/149‑1 (MQ=255)
                                                                                                         aCCGAGCAGCGCATCACCTTTGACATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGt                                  <  6:167515/149‑1 (MQ=255)
                                                                                                               cagcGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGa                            <  3:168924/149‑1 (MQ=255)
                                                                                                                  cgcATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAgcg                         >  5:113139/1‑149 (MQ=255)
                                                                                                                   gcATCACCTTTCCCATCGAGCCCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAgcgc                        <  6:284252/149‑1 (MQ=255)
                                                                                                                          cTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGTCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAgg                 <  1:97864/149‑1 (MQ=255)
                                                                                                                           ccagCCATCGCGGTCTCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGt                >  6:187783/4‑149 (MQ=255)
                                                                                                                           ccagCCATCGCGGTCTCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGt                >  7:258806/4‑149 (MQ=255)
                                                                                                                            cagCCATCGCGGTCTCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCACATGGTCAACGAGCGCCTGCAGGTg               >  7:352042/3‑149 (MQ=255)
                                                                                                                             agCCATCGCGGTCTCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTgg              <  5:187783/148‑1 (MQ=255)
                                                                                                                              gCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCGGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCGCGCCAACTGGGCAGCGCGCGCCTGCAGGTgga             >  6:266618/1‑148 (MQ=255)
                                                                                                                                cATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGat                                                         <  3:246853/103‑1 (MQ=255)
                                                                                                                                cATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGat                                                         >  4:246853/1‑103 (MQ=255)
                                                                                                                                 aTCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTGGCCcg          <  6:328855/149‑1 (MQ=255)
                                                                                                                                  tCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTGGCCcgc         >  5:1998/1‑149 (MQ=255)
                                                                                                                                    gagaCCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTGGCCCGCGa       <  2:50824/149‑1 (MQ=255)
                                                                                                                                       aCCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTGGCCCGCGAGc     >  6:69604/1‑148 (MQ=255)
                                                                                                                                       aCCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTGGCCCGCGAGc     >  4:63908/1‑148 (MQ=255)
                                                                                                                                         cgcgATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAATGAGCGCCTGCAGGTGGCCCGCGa       <  7:77368/144‑1 (MQ=255)
                                                                                                                                         cgcgATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAATGAGCGCCTGCAGGTGGCCCGCGa       >  8:77368/1‑144 (MQ=255)
                                                                                                                                         cgcgATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACCGTCATTTTCGATGATGGCAACGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTGGCCCGCGAGCAAt  >  5:161438/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GTGGGTATCCACAGCTACCAGAAACTGCCGATCGATGCTGTGCCGGACATCACTAACGTCCAGGTGCAGATCAACACCGCCGCCCCCGGCTATTCGCCGCTGGAGACCGAGCAGCGCATCACCTTTGCCATCGAGACCGCGATGGCTGGCCTGCCAGGGCTCAAGCAGACCCGCTCGCTGTCGCGCTCGGGCCTGTCGCAAGTCACGGTCATTTTCGATGATGGCACCGATATCTTCTTCGCCCGCCAACTGGTCAACGAGCGCCTGCAGGTGGCCCGCGAGCAAT  >  NC_002947/2753750‑2754035

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