Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,560,455 +G intergenic (‑2/+1) glpR ← / ← glpR pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,560,4551.G98.2% 178.7 / ‑4.1 56intergenic (‑2/+1)glpR/glpRpseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon/pseudogene, DNA‑binding transcriptional repressor;regulator; Energy metabolism, carbon: Anaerobic respiration; repressor of the glp operon
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (27/28);  total (27/29)
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.99e-01

GATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTG  >  NC_000913/3560316‑3560587
                                                                                                                                             |                                                                                                                                     
gATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGcgca                                                                                                                                     >  1:194976/1‑140 (MQ=255)
     ccAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                <  1:396979/142‑1 (MQ=255)
     aCAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATc                                                                                                                                <  1:145261/143‑1 (MQ=255)
       aGCGTCGAGCCCTGGGGGATGTGTGCCGCCACTTGGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGGGGCCTGGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAg                                                                                                                              <  1:225285/143‑1 (MQ=255)
         cGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGt                                                                                                                            >  1:9460/1‑143 (MQ=255)
           tCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGttt                                                                                                                          <  1:337497/143‑1 (MQ=255)
             gggccggtGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCGTTTTCGTCGGTCGGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                        <  2:986974/136‑1 (MQ=255)
             gAGCCTTTGGGGTTTTGCTCCGCCCCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                        <  2:689054/143‑1 (MQ=255)
             gAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                                                                                                        <  1:1103346/143‑1 (MQ=255)
                  aTTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCCCCATGATGGCGCAGGATCAGGTTTTGCTCCg                                                                                                                   >  1:712277/1‑143 (MQ=255)
                          tttGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGCTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                                                                                           <  1:595938/143‑1 (MQ=255)
                            tGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCTCCGCACCTCCATGATGGCGAAGGAACAGTATTTACTCCGCCAGCTCAtt                                                                                                         <  1:854990/143‑1 (MQ=255)
                               tCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAg                                                                                                      >  2:915680/1‑143 (MQ=255)
                                 cGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTACGCCAGCTCATTGAggc                                                                                                    >  2:923535/1‑142 (MQ=255)
                                  gCCACTTTGCGGGCGATGCGCTCTTTTCCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTc                                                                                                   <  1:759629/143‑1 (MQ=255)
                                   ccACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg                                                                                                  >  1:185179/1‑143 (MQ=255)
                                     aCTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgcg                                                                                                >  1:996290/1‑143 (MQ=255)
                                      tctcgtgggcggggcgCTCCGTTTCTTCGGTCTGGGTGGCCTTGCGGTCGGGCCACGGCGTGTTCACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                                                                                              <  2:683874/137‑1 (MQ=255)
                                         tGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTGTTTGTCCGCCAGCTCATTGAGGTCGCGGCGa                                                                                           >  1:902594/1‑143 (MQ=255)
                                            gggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATa                                                                                        <  1:1122581/143‑1 (MQ=255)
                                             ggCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTCAACCGAACTGGAACGCAGCGCCGCACCGCCAGGATGGCGCAGGAGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAg                                                                                       >  2:1112237/1‑143 (MQ=255)
                                              gCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATCGAGGTCGCGGCGAATCGt                                                                                      >  2:99013/1‑143 (MQ=255)
                                                  tGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGTTCGTGCCACGGCGTGTTAACCGAACTGGATGGCAGCGCCGCACTGCTATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGc                                                                                  <  2:322598/143‑1 (MQ=255)
                                                    cgcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCGACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                <  1:486858/143‑1 (MQ=255)
                                                    cgcTCGTTGTCTTCGGTCTGGGTGGCGTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGAGGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCgg                                                                                <  2:33621/143‑1 (MQ=255)
                                                     gcTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCGGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg                                                                               >  1:1078536/1‑143 (MQ=255)
                                                        ctTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg                                                                            >  2:893896/1‑143 (MQ=255)
                                                           tttCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGACTAGTCTGCGGGCTGCCg                                                                         >  2:337608/1‑143 (MQ=255)
                                                           tgtcttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg                                                                         <  1:528023/141‑1 (MQ=255)
                                                               ttcGGTCTGGGTGGCCTTGCGGTCGTGCCGCGGCGTGTGAGCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                     <  1:517885/143‑1 (MQ=255)
                                                               ttcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCCTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGgaga                                                                     >  1:472447/1‑143 (MQ=255)
                                                                tcGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGCTCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGGGAGTTGCCGTCGTGCTGACggggga                                                                    >  1:782824/1‑139 (MQ=255)
                                                                 cGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGaaa                                                                   <  2:986089/143‑1 (MQ=255)
                                                                         gTGGCCCTGCGAGCGCGCCCCGGCGGGTTGGCCGAGCTGGGAGGCAGCGCCGCAACGCCGTGATGGAGCAGTTGCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGCATAGTCCGCGGGCTGACGGAGAAATGCTCTAc                                                           <  2:76624/143‑1 (MQ=255)
                                                                          gcgCCGTGCGATCGGGCCACGGCGTGTTATCCGAACGGGGAGGCCGCGCCGCACCGCCATGATGGCGCAGGAGCAGGTTTTGCCCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTAcc                                                          <  2:53319/141‑1 (MQ=255)
                                                                             ccgtGCGATCGTGCCACGGCGTGTTAACCGAGCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGc                                                       <  2:1100815/140‑1 (MQ=255)
                                                                                    aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGCTGGGATCAGCGCGGGCTGCCGGAGAAAGGCTCTACCAGCTCTTCga                                                >  2:1128499/1‑142 (MQ=255)
                                                                                        cgCCACGGCGCGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACt                                            <  1:468635/142‑1 (MQ=255)
                                                                                         gCCACGGCGTGTTAGCGGCACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                           <  2:784641/143‑1 (MQ=255)
                                                                                          ccACGGCGTGTTAACCGAACCGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGa                                          <  1:144922/143‑1 (MQ=255)
                                                                                           cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGAc                                         <  1:1082303/143‑1 (MQ=255)
                                                                                            aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACa                                        >  1:553511/1‑143 (MQ=255)
                                                                                                   gtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCGCATTGAAACAACGCATGAGTGTGTTTGTGCTGCCGGAGCAATGCGCTACCAGCTCTCCGGTCCTGACATAACCct                                 >  1:538623/1‑143 (MQ=255)
                                                                                                        ccGAACTGGAAGGCCGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGGAATGATCTAGCAGCTCTGCGGTACTGACATAACCCTCCTGt                            >  2:221688/1‑143 (MQ=255)
                                                                                                        ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGt                            <  1:962081/143‑1 (MQ=255)
                                                                                                         cGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCGTTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGACATGCTCTACCAGCTCTCCGGTACTGACGTAACCGTGCTGtc                           >  2:46346/1‑142 (MQ=255)
                                                                                                           aaCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGGCGGAGAAATGCTCTACCAGCTCTTCGGTACTGGCATAACCCTGCTGTTTa                         >  1:1051621/1‑143 (MQ=255)
                                                                                                              tGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTCTGCTCCGCCGGCTCCTTGAGGTCGCGGCGCATAGTCTGCGGGCTGCCGGAGAAATGCTCTACCAGCTCTACGGTACTGACATAACCCTGCTGTTTAAcc                      >  1:288268/1‑143 (MQ=255)
                                                                                                               ggAAGGCAGTGCCGCCCCGCCAGGATGGCGCAGGATCAGGTTGTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCa                     <  1:1112237/143‑1 (MQ=255)
                                                                                                                 aaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                           >  1:129316/1‑119 (MQ=255)
                                                                                                                 aaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTg                                           <  2:129316/119‑1 (MQ=255)
                                                                                                                      aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGCGAAATGCTCTACCAGCGCTTCGGTTCTGCCATACCCCTGCTGTTCTACCAGCTCCat              >  1:20224/1‑143 (MQ=255)
                                                                                                                       gcgcCGCCCCGCCAGGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGata             <  1:719929/143‑1 (MQ=255)
                                                                                                                         gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGataat           >  1:101182/1‑143 (MQ=255)
                                                                                                                            gcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCCGCGACTAGTCGGCGGGCTGACGAACAAATGCTCTACCAGCTCTTCGGTACCGACATAACCCTCCTGTGTAACCCGTTCGAAAAGAcc        >  2:289497/1‑143 (MQ=255)
                                                                                                                                cGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACCGCGCAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTtg    >  2:632721/1‑143 (MQ=255)
                                                                                                                                  ccATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGCGATAGGCCTTCGGTTTAACCAGTTCGATAATACCGTtgtg  >  1:486009/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                     
GATAAACAGCGTCGAGCCATTGGGGATTTGCTCCGCCACTTTGCGGGCGATGCGCTCTTTTTCTTCGGTCTGGGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACGGAGAAATGCTCTACCAGCTCTTCGGTACTGACATAACCCTGCTGTTTAACCAGTTCGATAATACCGTTGTG  >  NC_000913/3560316‑3560587

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