Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,164,719 C→T 7.4% intergenic (+69/+48) PP_3665 → / ← PP_3666 AraC family transcriptional regulator/MHS family metabolite MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,164,7190CT 7.4% 157.1 / 2.3 54intergenic (+69/+48)PP_3665/PP_3666AraC family transcriptional regulator/MHS family metabolite MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/26)
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.99e-01

CTGCGCATTTTTCCAAGAGTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGCAAGTATCCTGCCGGGGCCAGCTTG  >  NC_002947/4164571‑4164866
                                                                                                                                                     |                                                                                                                                                   
cggcgcATTTTTCCAAGAGTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTc                                                                                                                                                     <  8:280033/147‑1 (MQ=255)
   cgcATTTTTCCAAGAGTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTcc                                                                                                                                                    <  1:346853/147‑1 (MQ=255)
   cgcATTTTTCCAAGAGTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGc                                                                                                                                                  >  1:38720/1‑149 (MQ=255)
       tttttCCAAGAGTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCcaca                                                                                                                                              <  2:213914/149‑1 (MQ=255)
              aaGAGTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGtt                                                                                                                                       <  4:255505/149‑1 (MQ=255)
                 agTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTa                                                                                                                                    <  1:87271/149‑1 (MQ=255)
                  gTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAg                                                                                                                                   <  6:305268/149‑1 (MQ=255)
                   ttATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAg                                                                                                                                   <  5:343198/148‑1 (MQ=255)
                   ttATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGAGCCTACACGCACATCGGTGGTCAGCCACCGGCCGTTGTAg                                                                                                                                   >  6:343198/1‑148 (MQ=255)
                         ccGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGt                                                                                                                             >  8:74756/1‑148 (MQ=255)
                         ccGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGt                                                                                                                             >  5:177330/1‑148 (MQ=255)
                           ggTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTgag                                                                                                                          <  7:74756/149‑1 (MQ=255)
                            gTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTgag                                                                                                                          >  7:127433/1‑148 (MQ=255)
                             ttcttcGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGTCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGTCGCGTGagaa                                                                                                                        <  6:177330/149‑1 (MQ=255)
                             ttcttcGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGagaa                                                                                                                        <  2:38720/149‑1 (MQ=255)
                                ttcGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTTCGCCACAGGCCGTTGTAGGCGCGTGagaag                                                                                                                       >  7:140943/1‑147 (MQ=255)
                                ttcGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTTCGCCACAGGCCGTTGTAGGCGCGTGagaag                                                                                                                       >  7:136235/1‑147 (MQ=255)
                                         tCGCCTTTAGGGGAGCGTCGGCAGGCCTCTCTCCACTGAGGTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCTCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTg                                                                                                              <  7:308983/147‑1 (MQ=255)
                                                tAGGGGAGCGGCGGCAGGCCTCTCTCAACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAAc                                                                                                     >  1:9754/1‑149 (MQ=255)
                                                     gAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAgg                                                                                                                                            <  5:62801/105‑1 (MQ=255)
                                                     gAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAgg                                                                                                                                            >  6:62801/1‑105 (MQ=255)
                                                     gAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATc                                                                                                >  4:32395/1‑149 (MQ=255)
                                                      aGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACAGCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCt                                                                                               >  7:189585/1‑149 (MQ=255)
                                                          gcggcAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTc                                                                                             >  4:236808/1‑147 (MQ=255)
                                                                 gCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCcgc                                                                                    >  2:2026/1‑149 (MQ=255)
                                                                   ctctctCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGc                                                                                       <  3:338388/144‑1 (MQ=255)
                                                                   ctctctCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGc                                                                                       >  4:338388/1‑144 (MQ=255)
                                                                        tcCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTgc                                                                                <  7:291297/146‑1 (MQ=255)
                                                                          cACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGt                                                                           <  3:32395/149‑1 (MQ=255)
                                                                               aGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTccc                                                                      <  2:173495/149‑1 (MQ=255)
                                                                               aGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTccc                                                                      <  3:236808/149‑1 (MQ=255)
                                                                                        aGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACAt                                                               >  8:325970/1‑147 (MQ=255)
                                                                                             ttGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATg                                                        >  2:257432/1‑149 (MQ=255)
                                                                                                gCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGc                                                      <  8:127433/148‑1 (MQ=255)
                                                                                                      cTGGCCTCTTCGCGGGCAAGCCCGCGCTCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAg                                                >  1:197823/1‑149 (MQ=255)
                                                                                                      cTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAgc                                               >  4:2712/1‑149 (MQ=255)
                                                                                                         gCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAgcgc                                             >  7:72314/1‑148 (MQ=255)
                                                                                                         gCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAgcg                                              <  8:189585/147‑1 (MQ=255)
                                                                                                          tgtaggcgcgGGCTTGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGcca                                           >  1:75907/7‑149 (MQ=255)
                                                                                                           ctctTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGccac                                          <  1:2026/149‑1 (MQ=255)
                                                                                                                cgcgGGCAAGCCCGCGCCTACACGCACATCGGTGGTTCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGccacca                                        <  8:140943/146‑1 (MQ=255)
                                                                                                                cgcgGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGt                                      <  1:257432/148‑1 (MQ=255)
                                                                                                                cgcgGGCAAGCACGCGCCTACACGCACATCGGTGGTTCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGccacca                                        <  8:136235/146‑1 (MQ=255)
                                                                                                                   gggCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCCTGGCCACCAGCGCCACCAGtgctg                                  >  6:257935/1‑149 (MQ=255)
                                                                                                                     gCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGaca                                                                  >  6:74793/1‑115 (MQ=255)
                                                                                                                     gCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGaca                                                                  <  5:74793/115‑1 (MQ=255)
                                                                                                                     gCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTgcc                                >  3:87586/1‑149 (MQ=255)
                                                                                                                          cccGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAg                           >  5:307051/1‑149 (MQ=255)
                                                                                                                          cccGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAg                           >  7:47245/1‑149 (MQ=255)
                                                                                                                             gcgcCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGc                          >  2:167398/1‑147 (MQ=255)
                                                                                                                             gcgcCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGCaa                        <  2:75907/149‑1 (MQ=255)
                                                                                                                                  cTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGCAAGTAt                    >  4:225322/1‑149 (MQ=255)
                                                                                                                                     cacGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGCAAGTATc                   <  3:2712/147‑1 (MQ=255)
                                                                                                                                      acGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGCAAGTATCCTg                >  4:213165/1‑149 (MQ=255)
                                                                                                                                            aTCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGCAAGTATCCTGCCgggg          <  4:87586/149‑1 (MQ=255)
                                                                                                                                                    tCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGCAAGTATCCTGCCGGGGCCAGCTtg  <  3:213165/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTGCGCATTTTTCCAAGAGTTATAGCCGGTTCTTCGGGCTGTCGCCTTTAGGGGAGCGGCGGCAGGCCTCTCTCCACTGAGTTGGGTCAGGTATTGGCATTGCTGGCCTCTTCGCGGGCAAGCCCGCGCCTACACGCACATCGGTGGTCCGCCACAGGCCGTTGTAGGCGCGTGAGAAGAAGCCAGTGCAGGTCAACTAATCTTCAAGGGCCGCGTGCGCCGTCTCCCGACACATCAGCATGGCCACCAGCGCCACCAGTGCTGCCGCCAGCAAGTATCCTGCCGGGGCCAGCTTG  >  NC_002947/4164571‑4164866

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: