Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,449 T→C 7.2% intergenic (+14/‑78) gpt → / → frsA xanthine phosphoribosyltransferase; xanthine‑guanine phosphoribosyltransferase/fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913256,4490TC 7.2% 155.5 / 3.5 56intergenic (+14/‑78)gpt/frsAxanthine phosphoribosyltransferase; xanthine‑guanine phosphoribosyltransferase/fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity
Reads supporting (aligned to +/- strand):  ref base T (28/23);  new base C (2/2);  total (30/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 = 2.84e-01

TCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAG  >  NC_000913/256311‑256596
                                                                                                                                          |                                                                                                                                                   
tgttagCAAAACCGGCGGGGGGTCCGCTGGTTGATGACCATGTAGTTGACATCCCGCAAGATACCTGGATGGAACCGCCGTGGTATAGGGTCGTCGTATGCGTCCCGCCAAGCTCCGGTCGCTAATCTTTGTAACGCCGGGCACTGCCg                                                                                                                                           <  2:223112/144‑1 (MQ=255)
 cctCGCAAACCCGGCTGGTCGTCCGCTGGTTGCTGACTATGTTGTTGATATCGCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGCATTCGTCCCGCCAATCTCCGGTCTCTAATCTTTTCAACGCCTGGCACTGCCgg                                                                                                                                          <  2:438772/147‑1 (MQ=255)
  ttCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGAACGATACCTGGGTTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCggg                                                                                                                                         >  1:623964/1‑149 (MQ=255)
   tCGCAAAACCGGCGGGTCGGCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAAGCTCGGGTCGCGAATCTTTTCAACGCCTGGCACTGCCGGGc                                                                                                                                        <  1:794541/149‑1 (MQ=255)
      cAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCgtt                                                                                                                                     <  2:416263/149‑1 (MQ=255)
        aacccGGCTGGTCGTCCGCTGGTTGATGACTATGGTGTTGATACCCCGCAAGATACCTGGCTTGAACAGCCGTGGGATATGGGCGTCGTATTCGTTCCGCCAATCTCCGGTCGCTAATCTTTACAACGCCCGGCAGTGCCAGGCgttga                                                                                                                                   <  2:438254/146‑2 (MQ=255)
             ggCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCttt                                                                                                                              <  1:87086/149‑1 (MQ=255)
              gCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCCGGCAGTGCCAGGCgttgaaaagat                                                                                                                             <  1:342729/149‑8 (MQ=255)
                 ggTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTaa                                                                                                                          >  1:38350/1‑149 (MQ=255)
                   tcgtcCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACCGCAATCAACCCTTGTCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCGTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACt                                                                                                                        >  1:750969/1‑149 (MQ=255)
                     gtcCGCTGGTTGATGACTACGTTGTTGGGATCCAGCAAGCTGCCTGGATTGGGCAGCCGTGGGATATGGGCGTCGTATTCGTCTCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCTGGCGTTGTTCTTTTTAGCTTc                                                                                                                      <  2:662058/149‑1 (MQ=255)
                      accGCTTGTGGATGACTATGTTGTTGATACCGCGCACGATACCTGGATTGAACGGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCGTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCa                                                                                                                     <  2:436822/148‑1 (MQ=255)
                           tGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCgg                                                                                                                <  1:457374/149‑1 (MQ=255)
                                   gaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCTTTGTTCTTTTTAACTTCAGGCGGGTTaacag                                                                                                        >  2:702542/1‑145 (MQ=255)
                                   gaCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACaag                                                                                                        >  1:51642/1‑148 (MQ=255)
                                    aCTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTAAGGCGGGTTACAATa                                                                                                       >  2:519265/1‑149 (MQ=255)
                                      tATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAgc                                                                                                     >  2:305596/1‑148 (MQ=255)
                                      tATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGt                                                                                                     >  1:644761/1‑149 (MQ=255)
                                      tATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGt                                                                                                     >  2:743970/1‑149 (MQ=255)
                                      tATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGGTACAATAGt                                                                                                     >  1:213269/1‑149 (MQ=255)
                                          ttgttgATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGGTACGATAGCTTcc                                                                                                 >  2:457296/1‑149 (MQ=255)
                                           tgttgATAGCTCGCAAGATACCTGGTTTGTGCTGCTGTTGGATATGGGCGTTGTATTCGTCCCGCTAATTTCCGGTTGCTAATCTTTTCAACGCCTGGCACGGCGGGGCGTTGTTCTTTTTATCTTCAGGCGGGTTACAATAGTTTCCa                                                                                                <  1:458572/149‑1 (MQ=255)
                                               gATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaa                                                                                            >  1:581055/1‑149 (MQ=255)
                                                agatCCCGCAAGTTACCTTGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCCACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCagtaag                                                                                           <  1:174055/147‑1 (MQ=255)
                                                       gCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTg                                                                                    <  1:631312/149‑1 (MQ=255)
                                                               cgtGGACTGAGCAGCCGTGGGATATGGGCGGCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTACCTTCGGGCGGGGGACAATAGTTTCCAGTAAGTATTCTGGAGGCTGc                                                                            <  2:646329/147‑1 (MQ=255)
                                                               ccTGGATTTGACGGCCGTGGGATATCGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGc                                                                            <  2:184268/149‑1 (MQ=255)
                                                                  ggATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATc                                                                         >  1:300347/1‑149 (MQ=255)
                                                                     ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTCCAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                      >  2:456764/1‑149 (MQ=255)
                                                                     ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTGTTCTGGAGGCTGCATCCAt                                                                      >  2:217443/1‑149 (MQ=255)
                                                                     ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                      >  1:251430/1‑149 (MQ=255)
                                                                     ttGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCAt                                                                      >  1:176751/1‑149 (MQ=255)
                                                                       gAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAACCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGa                                                                    <  1:148794/149‑1 (MQ=255)
                                                                           aGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGTAGGCTGCATCCATGAcaca                                                                >  1:67170/1‑149 (MQ=255)
                                                                           aGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGAcaca                                                                >  1:769662/1‑149 (MQ=255)
                                                                                      atGGGCGTCGCATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGa                                                     <  2:63586/149‑1 (MQ=255)
                                                                                                     gTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTGAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCa                                      >  1:65865/1‑149 (MQ=255)
                                                                                                        ccGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAAc                                   <  1:649070/149‑1 (MQ=255)
                                                                                                          gCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAAccc                                 >  2:434714/1‑149 (MQ=255)
                                                                                                            aaaTCTCCGGTCGCTAATCTTGTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCg                               <  1:145516/148‑1 (MQ=255)
                                                                                                                   cGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTaaa                        <  2:581055/149‑1 (MQ=255)
                                                                                                                   cGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTaaa                        <  1:47038/149‑1 (MQ=255)
                                                                                                                      tCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACAt                     <  1:706956/149‑1 (MQ=255)
                                                                                                                        gCTAATCTTTTCAACGCATAATTTTGTCGGGCGGTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTGTAAACATcc                   <  2:763504/149‑1 (MQ=255)
                                                                                                                         cTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCACACAACGGTTAGACTGTGTTCAAACCCCGCTTTAAACATCCt                  <  1:194179/149‑1 (MQ=255)
                                                                                                                           aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCCGTTCAAACCCCGCCTTAAACATCCCGa                >  2:795616/1‑149 (MQ=255)
                                                                                                                           aaTCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTGTTTAACTTCAGGCGGGTTACAATAGATTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGa                >  1:520091/1‑149 (MQ=255)
                                                                                                                            aaaagaacAACGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGTTGCCGCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGaa               >  1:786570/8‑149 (MQ=255)
                                                                                                                             tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGaaa              >  1:653471/1‑149 (MQ=255)
                                                                                                                             tCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCGGaaa              >  1:435028/1‑149 (MQ=255)
                                                                                                                                   cAACGGCTGGCACTTCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGa        <  1:727058/149‑1 (MQ=255)
                                                                                                                                   cAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGa        <  1:98536/149‑1 (MQ=255)
                                                                                                                                     aCGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGCCg      >  2:283706/1‑149 (MQ=255)
                                                                                                                                      cGCCCGGCAGTGCCAGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTCTAAACATCCTGAAACCCCCACGc     >  2:315766/1‑149 (MQ=255)
                                                                                                                                         cTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAg  >  1:274134/1‑149 (MQ=255)
                                                                                                                                         cTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAg  >  1:743790/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
TCTTCGCAAAACCGGCTGGTCGTCCGCTGGTTGATGACTATGTTGTTGATATCCCGCAAGATACCTGGATTGAACAGCCGTGGGATATGGGCGTCGTATTCGTCCCGCCAATCTCCGGTCGCTAATCTTTTCAACGCCTGGCACTGCCGGGCGTTGTTCTTTTTAACTTCAGGCGGGTTACAATAGTTTCCAGTAAGTATTCTGGAGGCTGCATCCATGACACAGGCAAACCTGAGCGAAACCCTGTTCAAACCCCGCTTTAAACATCCTGAAACCTCGACGCTAG  >  NC_000913/256311‑256596

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