Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,073,161:1 +G 7.1% coding (346/1308 nt) mdtA ← multidrug transporter membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,073,1611.G 7.1% 223.2 / 10.7 56coding (346/1308 nt)mdtAmultidrug transporter membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base . (30/22);  new base G (2/2);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

TGCTGTCTTCAGCGTACAGGCCTTTATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAGGGGGAAGTCGCCTATCCGCACCGGCACGAC  >  NC_002947/4073014‑4073302
                                                                                                                                                    |                                                                                                                                             
tGCTGTCTTCAGCGTACAGGCCTTTATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCt                                                                                                                                             <  1:7307/150‑1 (MQ=255)
     tCTTCAGCGTACAGGCCTTTATAGCGGGCCAGTTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGAcg                                                                                                                                        >  1:608712/1‑150 (MQ=255)
       ttCAGCGTACAGGCCTTTATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCTAAGGTGCCAACACGGAAGAGCAGGGGTATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCCGGTCGCCAGCCTTGAcgcg                                                                                                                                      >  1:461571/1‑150 (MQ=255)
           gCGTACAGGCCTTTATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTgg                                                                                                                                  >  1:553326/1‑150 (MQ=255)
              tACAGGCCTTTATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGcct                                                                                                                               >  1:285547/1‑150 (MQ=255)
                 aGGCCTTTATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGccttcc                                                                                                                            >  1:1016686/1‑150 (MQ=255)
                      tttATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCCAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGaa                                                                                                                       >  1:1243522/1‑150 (MQ=255)
                      tttATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGaa                                                                                                                       >  2:121208/1‑150 (MQ=255)
                      tttATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGaa                                                                                                                       <  1:149792/150‑1 (MQ=255)
                      tctataGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGaa                                                                                                                       >  1:141104/3‑150 (MQ=255)
                         ataGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGATGCGCTGGCCTTCCTTGAAGTg                                                                                                                    <  1:1268267/150‑1 (MQ=255)
                                           cTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCgcc                                                                                                  >  1:108441/1‑150 (MQ=255)
                                                   tcttcAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCAcgc                                                                                          <  2:608712/150‑1 (MQ=255)
                                                   tcttaAACTGTGCCTGGTTCTGCGCCAAGGCGCCTTCGGCCTGCGGCAGGTCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCAcgc                                                                                          <  2:1016686/150‑1 (MQ=255)
                                                       cAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTCCGCCAGCCACGCGGCt                                                                                      <  2:461571/150‑1 (MQ=255)
                                                                cTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACAtt                                                                             >  2:278422/1‑150 (MQ=255)
                                                                 tGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGGTTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTg                                                                            >  1:12778/1‑150 (MQ=255)
                                                                    ttCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGAcc                                                                         <  2:841669/150‑1 (MQ=255)
                                                                    ttCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGGCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGAcc                                                                         >  2:1155124/1‑150 (MQ=255)
                                                                     tCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGAGCGGGTCGATTTCGGCCAGCAGGTCGACAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCg                                                                        <  1:1387081/150‑1 (MQ=255)
                                                                       tGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTa                                                                      >  1:269748/1‑150 (MQ=255)
                                                                          gcCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAg                                                                                                                      <  1:926805/99‑1 (MQ=255)
                                                                          gcCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAg                                                                                                                      >  2:926805/1‑99 (MQ=255)
                                                                            cAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGc                                                                                                                                                <  2:1282556/72‑1 (MQ=255)
                                                                            cAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGc                                                                                                                                                >  1:1282556/1‑72 (MQ=255)
                                                                                 tGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCgg                                                            <  2:1243522/150‑1 (MQ=255)
                                                                                   ccTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTg                                                                                                                         >  1:60750/1‑87 (MQ=255)
                                                                                   ccTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTg                                                                                                                         <  2:60750/87‑1 (MQ=255)
                                                                                   ccTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTa                                                          >  2:1219476/1‑150 (MQ=255)
                                                                                     ttCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAAc                                                        >  1:82502/1‑150 (MQ=255)
                                                                                      tCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAAcc                                                       <  1:121208/150‑1 (MQ=255)
                                                                                        ggCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGt                                                     <  2:315875/150‑1 (MQ=255)
                                                                                        ggCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGt                                                     >  1:315875/1‑150 (MQ=255)
                                                                                         gCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTg                                                                            <  1:651584/126‑1 (MQ=255)
                                                                                         gCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTg                                                                            >  2:651584/1‑126 (MQ=255)
                                                                                          ccTGCTGCAGGGCGATGCGATAAGGGCGCGGCTAGATTTCGGAGAGCAGCTCGCCCGCGTTGACGCTCTGGCCGTCCGGGAACGGGATGTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGc                                                   >  2:1230933/1‑150 (MQ=255)
                                                                                                    ggCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCttt                                         <  2:553326/150‑1 (MQ=255)
                                                                                                        atgcgatAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTgtag                                     >  1:990622/1‑150 (MQ=255)
                                                                                                         tgcgatAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTgtagt                                    <  1:1219476/150‑1 (MQ=255)
                                                                                                             atAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTagag                                >  1:1266039/1‑150 (MQ=255)
                                                                                                               aaGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAggg                              <  1:1301316/150‑1 (MQ=255)
                                                                                                                  gaactggtgaaaatccacttcaaggaagGCCAGCGCGTCAAGGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAgg                               <  1:172481‑M1/118‑1 (MQ=255)
                                                                                                                  gaactggtgaaaatccacttcaaggaagGCCAGCGCGTCAAGGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAgg                               >  2:172481‑M1/29‑146 (MQ=255)
                                                                                                                           cGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACTGTGCCCAAGGCTTTGTAGTAGAGGGGGAAGTCGCCTa                  <  1:1155124/150‑1 (MQ=255)
                                                                                                                              tttCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAGGGGGAAGTCGCCTATcc               >  2:630093/1‑150 (MQ=255)
                                                                                                                                  ggCCAGCAGGTCGCCAGCCTTGACGCGNTGGCCTTCCTTGNAGTGGANTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCACGGCTTTGTAGTAGAGGGGGAAGTCGCCTATCCGCAc           >  1:1139595/1‑150 (MQ=255)
                                                                                                                                    ttcaaggaagGCCAGCGCGTCAAGGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAg                                              <  1:355632‑M1/103‑1 (MQ=255)
                                                                                                                                    ttcaaggaagGCCAGCGCGTCAAGGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAg                                              >  2:355632‑M1/11‑113 (MQ=255)
                                                                                                                                      agcagGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAGGGGGAAGTCGCCTATCCGCACCGGc       >  1:765371/1‑150 (MQ=255)
                                                                                                                                           gTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCATTTTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAGGGGGAAGTCGCCTATCCGCACCGGCacgac  >  1:100037/1‑150 (MQ=255)
                                                                                                                                           gTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAGGGGGGAGTCGCCTATCCGCACCGGCacgac  >  2:157242/1‑150 (MQ=255)
                                                                                                                                              GCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTaga                                 >  1:732814‑M1/1‑116 (MQ=255)
                                                                                                                                              GCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTaga                                 <  2:732814‑M1/116‑1 (MQ=255)
                                                                                                                                               CCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTagag                                <  2:1082317‑M1/116‑1 (MQ=255)
                                                                                                                                               CCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTagag                                >  1:1082317‑M1/1‑116 (MQ=255)
                                                                                                                                                catccTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTg                                                                   <  1:878611/77‑1 (MQ=255)
                                                                                                                                                CAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTg                                                                   <  2:878611‑M1/80‑1 (MQ=255)
                                                                                                                                                   CCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTg                                                    <  2:381135‑M1/92‑1 (MQ=255)
                                                                                                                                                   CCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTg                                                    >  1:381135‑M1/1‑92 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TGCTGTCTTCAGCGTACAGGCCTTTATAGCGGGCCAGGTCGACCTGAGCGTTCTTCAACTGTGCCTGGTTCTGCGCCAAGGTGCCTTCGGCCTGCTGCAGGGCGATGCGATAAGGGCGCGGGTCGATTTCGGCCAGCAGGTCGCCAGCCTTGACGCGCTGGCCTTCCTTGAAGTGGATTTTCACCAGTTCGCCAGCCACGCGGCTGCGCACATTGACCGTATTGGTCGCGGTAACCGTGCCCAAGGCTTTGTAGTAGAGGGGGAAGTCGCCTATCCGCACCGGCACGAC  >  NC_002947/4073014‑4073302

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