Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,843,542 G→A 7.2% R418C (CGC→TGC)  gacS ← sensor protein GacS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,843,5420GA 7.2% 157.3 / 2.3 55R418C (CGC→TGC) gacSsensor protein GacS
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base A (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

AGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGG  >  NC_002947/1843395‑1843690
                                                                                                                                                   |                                                                                                                                                    
aGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCgg                                                                                                                                                     <  4:70252/149‑1 (MQ=255)
  ttGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCggg                                                                                                                                                    <  2:304283/148‑1 (MQ=255)
     tCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGaa                                                                                                                                                >  2:150439/1‑149 (MQ=255)
         ccTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTg                                                                                                                                            <  1:150439/149‑1 (MQ=255)
             gCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATgg                                                                                                                                        <  4:239162/149‑1 (MQ=255)
              cTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGc                                                                                                                                       <  7:253677/149‑1 (MQ=255)
                                   cgcACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGa                                                                                                                   <  1:166124/148‑1 (MQ=255)
                                   cgcACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAg                                                                                                                  <  2:65916/149‑1 (MQ=255)
                                    gcACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGt                                                                                                                                                   <  8:281504/115‑1 (MQ=255)
                                    gcACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGt                                                                                                                                                   >  7:281504/1‑115 (MQ=255)
                                        aTCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCt                                                                                                             >  6:320947/1‑149 (MQ=255)
                                         tCCTGAGGGTACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGAGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTg                                                                                                            <  3:266126/149‑1 (MQ=255)
                                                     aGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCaa                                                                                                >  2:64952/1‑149 (MQ=255)
                                                      ggCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGATGATCTGCTTGAGCCGCAAc                                                                                               <  6:110558/149‑1 (MQ=255)
                                                       gCCGATGCCGGTGTCCTGCACGCTTATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACg                                                                                              >  7:49979/1‑149 (MQ=255)
                                                         cGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGa                                                                                                                   >  2:9251/1‑126 (MQ=255)
                                                         cGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTCGATGGCGTTGCTGACCAGGTTGGTGa                                                                                                                   <  1:9251/126‑1 (MQ=255)
                                                               cGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGAATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGcc                                                                                       >  7:326519/1‑148 (MQ=255)
                                                                 gtgtCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGt                                                                                                                         <  6:277835/112‑1 (MQ=255)
                                                                 gtgtCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGt                                                                                                                         >  5:277835/1‑112 (MQ=255)
                                                                 gaagacTGCGCGCAGTTGCGCATCAGCGTGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTAc                                                                                    >  8:146907/6‑149 (MQ=255)
                                                                      cTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGACTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGc                                                                                >  8:15764/1‑148 (MQ=255)
                                                                                aTGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgt                                                                       <  5:320947/147‑1 (MQ=255)
                                                                                  gcgcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGCCCTTGATGGCGTTGCTGACCAGGGTGGTGAGGATCTGCTTGAGCCGCCACGGGTCGCCTACCAGCGAAGCAGgtgtgt                                                                   >  5:274622/1‑149 (MQ=255)
                                                                                    gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAg                                                                         >  5:19242/1‑141 (MQ=255)
                                                                                    gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAg                                                                         <  6:19242/141‑1 (MQ=255)
                                                                                      aaCTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTcgcg                                                               >  3:353045/1‑149 (MQ=255)
                                                                                            gcgcAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCgg                                                              >  3:191238/1‑144 (MQ=255)
                                                                                                  tCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCAGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgt                                                                       <  3:208919/129‑1 (MQ=255)
                                                                                                  tCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCAGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgt                                                                       >  4:204278/1‑129 (MQ=255)
                                                                                                  tCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCAGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgt                                                                       >  4:208919/1‑129 (MQ=255)
                                                                                                  tCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCAGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgt                                                                       <  3:204278/129‑1 (MQ=255)
                                                                                                    ttCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCt                                                   >  3:315730/1‑147 (MQ=255)
                                                                                                      cGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTc                                                  <  1:64952/146‑1 (MQ=255)
                                                                                                        tGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCaa                                             >  8:53832/1‑149 (MQ=255)
                                                                                                         gCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCgg                                                              <  8:89592/131‑1 (MQ=255)
                                                                                                         gCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCgg                                                              >  7:89592/1‑131 (MQ=255)
                                                                                                            cGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGACAGTCCCGGAAGAACACCCTCTGCAAATc                                          >  2:142292/1‑148 (MQ=255)
                                                                                                             gTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCa                                        <  7:65351/149‑1 (MQ=255)
                                                                                                                  ccAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctgc                                   >  1:234755/1‑149 (MQ=255)
                                                                                                                  ccAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctgc                                   >  3:229339/1‑149 (MQ=255)
                                                                                                                   cAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACtt                                                                                                                                             <  6:56171/42‑1 (MQ=255)
                                                                                                                   cAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctgct                                  >  1:361709/1‑149 (MQ=255)
                                                                                                                     gCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTg                                >  2:131653/1‑149 (MQ=255)
                                                                                                                       aTGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGc                              <  8:49979/149‑1 (MQ=255)
                                                                                                                         ggCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGt                            >  2:153684/1‑149 (MQ=255)
                                                                                                                            gcgcACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTggg                         <  4:191238/149‑1 (MQ=255)
                                                                                                                             cgcACGACGATGGGGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGggc                        <  7:146907/149‑1 (MQ=255)
                                                                                                                                  gacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATc                                                        >  5:332905/1‑112 (MQ=255)
                                                                                                                                  gacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATc                                                        <  6:332905/112‑1 (MQ=255)
                                                                                                                                  gacgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCg                   <  1:153684/149‑1 (MQ=255)
                                                                                                                                     gaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAg                <  8:326519/149‑1 (MQ=255)
                                                                                                                                         gTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTTCTTGGCGTGGGCGGCCGGAGCGa             <  7:53832/148‑1 (MQ=255)
                                                                                                                                          tGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGGAGGTGTGTCGCTGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTggg                         <  7:76022/135‑1 (MQ=255)
                                                                                                                                          tGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAg            >  7:139663/1‑148 (MQ=255)
                                                                                                                                          tGCATTCGCGGGTGAACTTGATGGAGTTGCTGACCAGGTTGGTGAGGATCAGCTTGAGCCTCAACGGGTCGCCTACAAGCGAAGCAGGTGTGTCGCGGTAGATCAGCCTCAGCAATTCCAGCTGCTTGGCGTgtg                         >  8:76022/1‑133 (MQ=255)
                                                                                                                                            ccTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGAt         <  2:234755/149‑1 (MQ=255)
                                                                                                                                                   gggTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCTCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAgg  <  7:15764/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTGCGCACATCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGGCGGCCGGAGCGAGGATGGTCAGG  >  NC_002947/1843395‑1843690

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: