Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,068,520 G→C 7.0% R388P (CGC→CCC)  PP_5322 → metal ion transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,068,5200GC 7.0% 180.2 / 2.7 57R388P (CGC→CCC) PP_5322metal ion transporter
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base C (2/2);  total (26/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.59e-01

ATCGGCTACCTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGCATCATCATCAAGAAGATGAAAGGCCCCAAGA  >  NC_002947/6068375‑6068666
                                                                                                                                                 |                                                                                                                                                  
aTCGGCTACCTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCct                                                                                                                                                 >  6:264296/1‑149 (MQ=255)
  cGGCTACCTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCctgc                                                                                                                                               >  4:171075/1‑149 (MQ=255)
   ggCTACCTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCctgct                                                                                                                                              <  5:293073/149‑1 (MQ=255)
   ggCTACCTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCctgct                                                                                                                                              >  4:279347/1‑149 (MQ=255)
   ggCTACCTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCctgc                                                                                                                                               <  3:234162/148‑1 (MQ=255)
     cTACCTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTgg                                                                                                                                            <  3:171075/149‑1 (MQ=255)
        ccTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTa                                                                                                                                         >  8:77652/1‑149 (MQ=255)
            aCCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCg                                                                                                                                      <  3:279347/148‑1 (MQ=255)
             ccATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGAc                                                                                                                                    <  6:44848/149‑1 (MQ=255)
             ccATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGAc                                                                                                                                    <  4:59023/149‑1 (MQ=255)
                        gTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACAt                                                                                                                         >  1:214058/1‑149 (MQ=255)
                         tGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGATGGGTATCGACCTTGACCACATc                                                                                                                        >  1:100774/1‑149 (MQ=255)
                             ggAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGATACCCAGCAGGCGCTCGACCTTGAAca                                                                                                                            >  2:166762/1‑141 (MQ=38)
                             ggAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGATACCCAGCAGGCGCTCGACCTTGAAca                                                                                                                            <  1:166762/141‑1 (MQ=38)
                             ggAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAgg                                                                                                                    <  8:70652/149‑1 (MQ=255)
                              gAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGc                                                                                                                   <  6:194150/149‑1 (MQ=255)
                                  tGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGa                                                                                                                >  1:11717/1‑148 (MQ=255)
                                     tGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTg                                                                                                                                <  2:96067/129‑1 (MQ=255)
                                     tGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTg                                                                                                                                >  1:96067/1‑129 (MQ=255)
                                          ggCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGc                                                                                                       >  7:238001/1‑149 (MQ=255)
                                             gACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCgg                                                                                                    <  2:214058/149‑1 (MQ=255)
                                                  ccAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGATACCCAGCAGGCGCTCGACCTTGaac                                                                                                                             >  7:121034/1‑117 (MQ=37)
                                                  ccAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGATACCCAGCAGGCGCTCGACCTTGaac                                                                                                                             <  8:121034/119‑3 (MQ=37)
                                                    aGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGCGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATc                                                                                             <  5:216688/149‑1 (MQ=255)
                                                       aCGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTAc                                                                                          <  7:77652/149‑1 (MQ=255)
                                                         gAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACga                                                                                        >  5:20738/1‑149 (MQ=255)
                                                         gAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCCGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACga                                                                                        <  5:264296/149‑1 (MQ=255)
                                                           gCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGa                                                                                                                                     >  7:111429/1‑102 (MQ=255)
                                                           gCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGa                                                                                                                                     <  8:111429/102‑1 (MQ=255)
                                                            cACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGAc                                                                                     >  6:306973/1‑149 (MQ=255)
                                                             aCCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACg                                                                                    <  2:11717/149‑1 (MQ=255)
                                                               cGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCt                                                                                  >  2:35196/1‑149 (MQ=255)
                                                                  aaGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGTTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCaa                                                                               <  2:100774/149‑1 (MQ=255)
                                                                       cGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGcgcg                                                                          >  4:32036/1‑149 (MQ=255)
                                                                           cGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGTGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGcc                                                                      <  7:280266/149‑1 (MQ=255)
                                                                             tGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCt                                                                     >  6:107540/1‑148 (MQ=255)
                                                                               gCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTgaa                                                                  <  3:32036/149‑1 (MQ=255)
                                                                                       gCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAagaaga                                                             <  5:306973/146‑1 (MQ=255)
                                                                                        ccTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGaa                                                         <  2:2265/149‑1 (MQ=255)
                                                                                         cTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAg                                                        >  2:172900/1‑149 (MQ=255)
                                                                                              aGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTg                                                   <  7:142529/149‑1 (MQ=255)
                                                                                               gCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTgg                                                  <  8:238001/149‑1 (MQ=255)
                                                                                                            cgggTGCGTGGCGAGACCCCACTGTTCAATCTCGAGCGCCTGCTTGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTAGAATGGCt                                     <  8:33310/147‑1 (MQ=255)
                                                                                                               gTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCt                                                                     <  1:2226/114‑1 (MQ=255)
                                                                                                               gTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCt                                                                     >  2:2226/1‑114 (MQ=255)
                                                                                                                   gTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGGGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGcatc                              <  1:141918/149‑1 (MQ=255)
                                                                                                                     ggCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTgaagaa                                                               <  1:135574/114‑1 (MQ=255)
                                                                                                                     ggCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTgaagaa                                                               >  2:135574/1‑114 (MQ=255)
                                                                                                                      gCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGcatcatc                           <  6:20738/149‑1 (MQ=255)
                                                                                                                           acCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGCatcatcatc                        >  3:124846/1‑147 (MQ=255)
                                                                                                                             cccACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGCATCATCATCaag                     >  5:212888/1‑148 (MQ=255)
                                                                                                                             cccACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGCATCATCATCaag                     >  4:87651/1‑148 (MQ=255)
                                                                                                                             cccACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACAACGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGCATCATCATCaag                     >  1:226294/1‑148 (MQ=255)
                                                                                                                                       aaGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAg                                                >  4:180664/1‑111 (MQ=255)
                                                                                                                                       aaGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAg                                                <  3:180664/111‑1 (MQ=255)
                                                                                                                                          gTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGaa                                                         <  1:72388/99‑1 (MQ=255)
                                                                                                                                          gTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGaa                                                         >  2:72388/1‑99 (MQ=255)
                                                                                                                                               gcgcCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGCATCATCATCAAGAAGATGAAAGGcccc      >  4:115781/1‑145 (MQ=255)
                                                                                                                                               gcgcCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGCATCATCATCAAGAAGATGAAAGGCCCCAAGa  <  1:172900/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATCGGCTACCTGACCATGGAAGACGTGCTGGAAGTGCTGGTAGGCGACATCCAGGACGAGCACCGCAAGACCGAGCGTGGCATTCTTGCCTACCAGCCGGGCAAGCTGCTGGTGCGTGGCGACACCCCACTGTTCAAGGTCGAGCGCCTGCTGGGTATCGACCTTGACCACATCGAGGCCGAAACCCTGGCCGGGTTGATCTACGAGACGCTCAAGCGCGTGCCTGAAGAAGAGGAAGTGCTGGAGGTCGAAGGGCTGCGCATCATCATCAAGAAGATGAAAGGCCCCAAGA  >  NC_002947/6068375‑6068666

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: