Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,104,342 T→A 8.7% L276Q (CTG→CAG)  PP_2714 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,104,3420TA 8.7% 159.7 / 5.1 55L276Q (CTG→CAG) PP_2714sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (21/29);  new base A (2/3);  total (23/32)
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.37e-01

CTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACC  >  NC_002947/3104193‑3104490
                                                                                                                                                     |                                                                                                                                                    
cTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCt                                                                                                                                                      <  1:764524/150‑1 (MQ=255)
 ttGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTg                                                                                                                                                     >  2:237101/1‑150 (MQ=255)
  tGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTgg                                                                                                                                                    <  2:484149/150‑1 (MQ=255)
   gCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTggg                                                                                                                                                   <  2:219751/150‑1 (MQ=255)
    cTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGTTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGc                                                                                                                                                  <  1:1044704/150‑1 (MQ=255)
     cgcagcAGATCGTCGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCt                                                                                                                                                 >  1:366657/5‑150 (MQ=255)
       gcgcAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCttt                                                                                                                                               >  1:171058/1‑150 (MQ=255)
       cagcAGATCGTCGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCttt                                                                                                                                               <  2:366657/148‑1 (MQ=255)
          cAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCtt                                                                                                                                            <  1:91834/150‑1 (MQ=255)
          cAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCtt                                                                                                                                            <  1:654272/150‑1 (MQ=255)
          cAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCAGGGCTTTCtg                                                                                                                                            <  2:736599/150‑2 (MQ=255)
             aTCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTtgtg                                                                                                                                         <  2:547639/150‑1 (MQ=255)
                   gaTCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCggg                                                                                                                                   <  1:994591/150‑1 (MQ=255)
                          tgcGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATc                                                                                                                            <  1:237101/150‑1 (MQ=255)
                            cggcgTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTTGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATTCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCgg                                                                                                                          <  1:1123563/150‑1 (MQ=255)
                            cggcgTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCgg                                                                                                                          <  1:301906/150‑1 (MQ=255)
                              gcgTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCt                                                                                                                        <  1:195553/150‑1 (MQ=255)
                                  gAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGaca                                                                                                                    <  2:351087/150‑1 (MQ=255)
                                      aTCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACggg                                                                                                                >  2:654413/1‑150 (MQ=255)
                                         gCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCTAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGCCTGACACGggtgg                                                                                                             <  2:171058/150‑1 (MQ=255)
                                              aTTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAAc                                                                                                        >  1:293767/1‑150 (MQ=255)
                                               ttGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCAGGGCTTTCtg                                                                                                                                            <  2:802550/113‑2 (MQ=255)
                                               ttGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCAGGGCTTTCtg                                                                                                                                            >  1:802550/1‑112 (MQ=255)
                                                 gCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAg                                                                                                     <  1:654413/150‑1 (MQ=255)
                                                 gCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAg                                                                                                     <  2:665277/150‑1 (MQ=255)
                                                     aTGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACTGGTGGTCAACCAGTTGc                                                                                                 >  2:236092/1‑150 (MQ=255)
                                                           cccACGAAATGCGTACGCAGCTGGCCGCATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCt                                                                                           <  2:578855/150‑1 (MQ=255)
                                                                gAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCtt                                                                                          >  2:597677/1‑146 (MQ=255)
                                                                gAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCtt                                                                                          <  1:597677/146‑1 (MQ=255)
                                                                gAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGcc                                                                                      >  1:298209/1‑150 (MQ=255)
                                                                                  gCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTg                                                                    >  2:281728/1‑150 (MQ=255)
                                                                                      tATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCa                                                                                                      <  2:177020/112‑1 (MQ=255)
                                                                                      tATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCa                                                                                                      >  1:177020/1‑112 (MQ=255)
                                                                                      tATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGc                                                                >  1:682954/1‑150 (MQ=255)
                                                                                        ttGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAg                                                              >  2:545079/1‑150 (MQ=255)
                                                                                               ggCACACGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATGGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGGCGAGCCAAGACCGGGGCAGCGAGCCa                                                       >  2:313170/2‑150 (MQ=255)
                                                                                                 cACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAg                                                      <  1:56808/149‑1 (MQ=255)
                                                                                                 cACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAg                                                      >  2:56808/1‑149 (MQ=255)
                                                                                                  aCGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTg                                                                                                  <  1:250964/104‑1 (MQ=255)
                                                                                                  aCGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTg                                                                                                  >  2:250964/1‑104 (MQ=255)
                                                                                                             gcgcTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACTGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGAc                                         <  1:236092/150‑1 (MQ=255)
                                                                                                                        cAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGgt                              >  1:246203/1‑150 (MQ=255)
                                                                                                                        cAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGCTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACACGGATCGACCTGGCCGAGgt                              >  1:246243/1‑150 (MQ=255)
                                                                                                                         aGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGgtc                             <  2:293767/150‑1 (MQ=255)
                                                                                                                               gatgCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCAcc                       <  2:682954/150‑1 (MQ=255)
                                                                                                                                        cacaaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTgg                                                                   <  2:422068/93‑1 (MQ=255)
                                                                                                                                        cacaaGAAAGCCCAGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTgg                                                                   >  1:422068/5‑97 (MQ=255)
                                                                                                                                         gcCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGACTCTCTGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTgg             <  2:885467/150‑1 (MQ=255)
                                                                                                                                          ccAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAAc                  <  2:660924/144‑1 (MQ=255)
                                                                                                                                          ccAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAAc                  >  1:660924/1‑144 (MQ=255)
                                                                                                                                               aaGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAAc       <  1:501197/150‑1 (MQ=255)
                                                                                                                                                  cccTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCa                                                               <  2:523847/91‑1 (MQ=255)
                                                                                                                                                  cccTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCa                                                               >  1:523847/1‑91 (MQ=255)
                                                                                                                                                    cTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTgg             >  1:750955/1‑139 (MQ=255)
                                                                                                                                                    cTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTgg             <  2:750955/139‑1 (MQ=255)
                                                                                                                                                    cTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGAcc  >  1:705018/1‑150 (MQ=255)
                                                                                                                                                     tGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAg                                <  2:586332/119‑1 (MQ=255)
                                                                                                                                                     tGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAg                                >  1:586332/1‑119 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTGCTGGCGCAGATCGACGATCTGCTGCGGCGTGAACATCGCTTCATTGCCGATGCTGCCCACGAAATGCGTACGCCTCTGGCCGTATTGCGCCTGCACGCGCAGAACGCGCTGGCCGCCAGCAATGATGCCGAGCGCCAGAAAGCCCTGGGCTTTCTTGTGGGCGGGGTCGATCGGCTGACACGGGTGGTCAACCAGTTGCTGACCTTGGCCCGCGTCGAGCCAAGACTGGGGCAGCGAGCCAGCACCCGGATCGACCTGGCCGAGGTCGTCACCGAAACCCTGGCTGAACTGACC  >  NC_002947/3104193‑3104490

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: ‑