Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 133,209 A→T 6.9% *211L (TAG→TTG)  dsbA → thiol:disulfide interchange protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947133,2090AT 6.9% 176.3 / 3.1 58*211L (TAG→TTG) dsbAthiol:disulfide interchange protein
Reads supporting (aligned to +/- strand):  ref base A (25/29);  new base T (2/2);  total (27/31)
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.90e-01

AGGAATTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCA  >  NC_002947/133062‑133355
                                                                                                                                                   |                                                                                                                                                  
aGGAATTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAg                                                                                                                                                   >  2:124963/1‑149 (MQ=255)
    aTTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAggcgg                                                                                                                                               <  1:124963/149‑1 (MQ=255)
          gaagaaGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTGCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCAc                                                                                                                                         <  5:189985/149‑1 (MQ=255)
            agaagTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCAccc                                                                                                                                       <  6:268830/149‑1 (MQ=255)
             gaagTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCa                                                                                                                                      >  3:111536/1‑149 (MQ=255)
                gTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGc                                                                                                                                   <  8:104240/149‑1 (MQ=255)
                gTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGc                                                                                                                                   <  1:30863/149‑1 (MQ=255)
                gTATCAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGc                                                                                                                                   <  2:245221/149‑1 (MQ=255)
                     agaTCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCg                                                                                                                               <  2:175623/148‑1 (MQ=255)
                      gaTCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCt                                                                                                                             >  7:245441/1‑149 (MQ=255)
                       aTCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCtt                                                                                                                            <  6:93912/149‑1 (MQ=255)
                          aCCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCtt                                                                                                                            >  2:284128/1‑146 (MQ=255)
                           ccGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCaga                                                                                                                        <  8:128726/149‑1 (MQ=255)
                                  aCCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCCGCCTACTTGgcagcgc                                                                                                                                            <  4:140022/122‑2 (MQ=255)
                                  aCCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCCGCCTACTTGgcagcgc                                                                                                                                            >  3:140022/1‑121 (MQ=255)
                                  aCCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCCGCCTACTTGgcagcgc                                                                                                                                            >  3:134841/1‑121 (MQ=255)
                                  aCCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCCGCCTACTTGgcagcgc                                                                                                                                            <  4:134841/122‑2 (MQ=255)
                                    cGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGt                                                                                                               >  3:254238/1‑149 (MQ=255)
                                     gAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTg                                                                                                              <  1:178800/149‑1 (MQ=255)
                                         cTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCAt                                                                                                          >  8:253163/1‑149 (MQ=255)
                                            gtcgtcAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTgg                                                                                                       <  4:254238/149‑1 (MQ=255)
                                                 cAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTg                                                                                                   >  8:123050/1‑148 (MQ=255)
                                                   aCGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCAc                                                                                                <  7:253163/149‑1 (MQ=255)
                                                        aaGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGcc                                                                                           >  7:66576/1‑149 (MQ=255)
                                                        aaGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGcc                                                                                           <  8:66576/149‑1 (MQ=255)
                                                         aGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCg                                                                                          >  8:157104/1‑149 (MQ=255)
                                                            aCCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCaga                                                                                                                        >  5:358217/1‑116 (MQ=255)
                                                             ccGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGTCCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAgg                                                                                      <  1:284128/149‑1 (MQ=255)
                                                             ccGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAgg                                                                                      <  2:333747/149‑1 (MQ=255)
                                                                 ttCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCaa                                                                                  >  1:144367/1‑149 (MQ=255)
                                                                  tCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAgacggc                                                                                  >  5:33224/1‑143 (MQ=255)
                                                                   cGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCaa                                                                                  <  7:289850/147‑1 (MQ=255)
                                                                      ccTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACAAGc                                                                             <  7:157104/149‑1 (MQ=255)
                                                                                      ggACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGcccc                                                             >  6:305339/1‑149 (MQ=255)
                                                                                        aCCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCgg                                                           <  4:111536/149‑1 (MQ=255)
                                                                                            gAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCt                                                       >  2:227630/1‑149 (MQ=255)
                                                                                             aaGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTg                                                      >  7:335351/1‑149 (MQ=255)
                                                                                                gCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGccc                                                   >  5:178849/1‑149 (MQ=255)
                                                                                                        aaTGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCAc                                                                        <  1:68827/120‑1 (MQ=255)
                                                                                                        aaTGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCAc                                                                        >  2:68827/1‑120 (MQ=255)
                                                                                                              gCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCAcc                                                                       <  6:196501/115‑1 (MQ=255)
                                                                                                              gCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCAcc                                                                       >  5:196501/1‑115 (MQ=255)
                                                                                                                  aCCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCgg                                 <  2:144367/149‑1 (MQ=255)
                                                                                                                    cAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGa                                                 >  5:165440/1‑131 (MQ=255)
                                                                                                                    cAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGa                                                 <  6:165440/131‑1 (MQ=255)
                                                                                                                          aTCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCCTCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGc                                   <  7:43867/139‑1 (MQ=255)
                                                                                                                          aTCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGc                                   >  8:43867/1‑139 (MQ=255)
                                                                                                                              aCAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCa                          >  4:32360/1‑144 (MQ=255)
                                                                                                                              aCAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCa                          <  3:32360/144‑1 (MQ=255)
                                                                                                                                aaGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCt                                                       <  8:176870/113‑1 (MQ=255)
                                                                                                                                aaGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCt                                                       >  7:176870/1‑113 (MQ=255)
                                                                                                                                aaGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCaa                   <  8:335351/149‑1 (MQ=255)
                                                                                                                                         gcTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCAACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGt          <  5:305339/149‑1 (MQ=255)
                                                                                                                                             ccAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGgctgc                             <  4:242957/126‑1 (MQ=255)
                                                                                                                                             ccAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGgctgc                             >  3:242957/1‑126 (MQ=255)
                                                                                                                                             ccAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTggg       >  5:147325/1‑148 (MQ=255)
                                                                                                                                             ccAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTggg       >  6:104783/1‑148 (MQ=255)
                                                                                                                                                 gTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCa  <  6:178849/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AGGAATTGGCGAAGAAGTATGAGATCACCGGCGTACCGAGCCTGGTCGTCAACGGCAAGTACCGCTTCGACCTGGGGACTGCCGGCGGACCGGAAGGCGTGCTGAATGTCGCCGACCAGCTGATCGACAAGGAGCGCGCTGCCAAGTAGGCGGCGTCACCCATGCCCCGCTTCAGAGCTACGCGTGGCATTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCA  >  NC_002947/133062‑133355

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