Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 151,522 G→T 13.8% A26A (GCC→GCA yadK ← putative fimbrial protein YadK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913151,5220GT13.8% 131.9 / 14.1 58A26A (GCC→GCAyadKputative fimbrial protein YadK
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base T (3/5);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

GTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAA  >  NC_000913/151380‑151664
                                                                                                                                              |                                                                                                                                              
gTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATg                                                                                                                                                >  1:753975/1‑143 (MQ=255)
  ttaTAACTGTCAGCAGTAACATTACCTGAAAAATATCCCACCCCATCGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCGTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGc                                                                                                                                              <  2:318111/141‑1 (MQ=255)
   gATACCTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCt                                                                                                                                             <  2:202326/143‑1 (MQ=255)
        ccttCAGGAGTAACACTACCTTAAAAATATCCCCCCCTAACGACGCCAACATCAATATTACCCTGGTTACTGATTCCGCTTTTTCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                        <  1:789277/141‑1 (MQ=255)
        cctgcAGGAGTAACATTACCTGAAATATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                        <  2:152826/139‑1 (MQ=255)
        cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                        <  2:69320/143‑1 (MQ=255)
        cTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCCATTGCTGGACa                                                                                                                                        <  1:591174/143‑1 (MQ=255)
             aGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGa                                                                                                                                   <  2:549128/143‑1 (MQ=255)
              ggAGTTACAGCACCTGAAAAATATCCCACCCCAGCGACGCCAAGAGCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                                                                                                  <  2:251258/143‑1 (MQ=255)
                 gTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAACGGCTGGACAGGCGATGCa                                                                                                                               >  1:39991/1‑143 (MQ=255)
                  tcaCATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaa                                                                                                                              <  2:236247/141‑1 (MQ=255)
                     cATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCTGCAATGGCTGGACAGGCGATGCAAAAc                                                                                                                           >  2:422153/1‑143 (MQ=255)
                      attttcTGAAAAATAACCAACCCCAACGAGGCCAAGATCAAGATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCCATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACa                                                                                                                          <  1:350387/138‑1 (MQ=255)
                            tGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTATGCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTGTTTTACATTGGCCACAAGGTCGATATCCTGTCCAGCAACGGCTGGACAGGCGATGCAAAACAGTAATg                                                                                                                    >  2:637433/1‑143 (MQ=255)
                              aaaaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGac                                                                                                                  <  2:889843/143‑1 (MQ=255)
                                aaaTATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGCTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGCTGCAAAACAGTAATGacag                                                                                                                >  2:806962/1‑143 (MQ=255)
                                          ccccAACTACGCCAAGATCAATATTACCCTGTTTACTGATTCCGCTTTTGCTGGTGTTGTTTTTTATTTTGTCAACAAGGTCAATAGCCTGGCGAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatt                                                                                                      <  1:621386/143‑1 (MQ=255)
                                           cccAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                                                                                                     >  1:745132/1‑143 (MQ=255)
                                              aaCGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCCATGCAAAACAGTAATGACAGCAACAGAATTATtc                                                                                                  >  1:692691/1‑143 (MQ=255)
                                                 gacgCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                                                                                               >  2:415864/1‑143 (MQ=255)
                                                      cAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc                                                                                          >  2:161498/1‑143 (MQ=255)
                                                          aTCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGCCAGAAACAGAATTATTCTCTTCATCAGCt                                                                                      >  1:21494/1‑143 (MQ=255)
                                                             aaTATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACCGACACCGCATTATTCTCTTCATCAGCttgc                                                                                   >  1:281005/1‑141 (MQ=255)
                                                                    cccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTATTCATTAGCTTCCGTTGAGt                                                                            >  2:880831/1‑143 (MQ=255)
                                                                     ccTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTg                                                                           >  1:309361/1‑143 (MQ=255)
                                                                          ttACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGCCAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATg                                                                      <  1:702435/143‑1 (MQ=255)
                                                                            aCTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCGAGCAATGGCTGGACAGGCGATGCATAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                                                                    <  2:444121/143‑1 (MQ=255)
                                                                            aCTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGATCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAACCAGTCATGCCAGAAACAGCATTATTCTCTTCATCAGCTTACGTTGAGTCGGCTGCa                                                                    >  1:781471/1‑143 (MQ=255)
                                                                            aCTGATTCCGCTTTTGCAGGTGCTGCTGTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCa                                                                    <  2:644761/143‑1 (MQ=255)
                                                                              tGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCACTATCCTGTCCAGCAATGGCTGGACAGGCGGTGCAAAACAGTAATGACAGAAACAAAATTATTCTCTTCATCAGCTTACGTTGGGTGGGATGCata                                                                  >  1:819385/1‑143 (MQ=255)
                                                                                aTTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGGGGGATGCagata                                                                >  2:178790/1‑143 (MQ=255)
                                                                                     gCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATCTATAGtt                                                           >  1:623924/1‑143 (MQ=255)
                                                                                      cTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTc                                                          <  1:149401/143‑1 (MQ=255)
                                                                                       ttttGCAGGTGCGGTTGTCTACGTTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc                                                         <  1:26344/143‑1 (MQ=255)
                                                                                       ctttGCAGGTGCTGCGTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTcc                                                         <  1:361316/142‑1 (MQ=255)
                                                                                        tttGCAGGTGCTGTTTTTTACATTGGCAACACGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGACGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCACCCGCTACCGCTGAGTGGGAGCCATATCTAGTTCCt                                                        >  2:132552/1‑143 (MQ=255)
                                                                                            cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGCAt                                                    >  2:815750/1‑143 (MQ=255)
                                                                                            cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt                                                    <  2:694658/143‑1 (MQ=255)
                                                                                            cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAt                                                    >  1:532980/1‑143 (MQ=255)
                                                                                                 gcggTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTa                                               <  2:667581/140‑1 (MQ=255)
                                                                                                  cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTaa                                              <  1:247472/143‑1 (MQ=255)
                                                                                                     ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAACACAGTAATGACAGAAACAGAATTATTCTCTTCATGAGCTTCCGCTGAGTACGCTGCAAATATAGTTCCTGAATCCGTAAGTc                                           >  2:630735/1‑143 (MQ=255)
                                                                                                     ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTc                                           >  1:822668/1‑143 (MQ=255)
                                                                                                            cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAACTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttagt                                    <  2:447944/143‑1 (MQ=255)
                                                                                                              ttGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTgttcgtta                                  >  1:154906/1‑143 (MQ=255)
                                                                                                              ttGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATGATTCTCTTCATCAGCTTACGTAGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGtgtgacagtga                                  >  2:735100/1‑136 (MQ=255)
                                                                                                                ggCAACAAGGTCAATTTCCTGCCCAGCAATGGCGTGACAGGCGATCCAAAACAGTAATGACAGAAACAGAATTATTCTCTGCATCAGCTTACGTTGAGTTGGATGCATATATAGTTCCTGAATCTGTAAGTCTCTTAGTtata                                <  2:714815/143‑1 (MQ=255)
                                                                                                                 gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATTTCGTTCCTGAATCTGTAAGTCTGTTAGTTATAg                               <  1:89890/143‑1 (MQ=255)
                                                                                                                     caaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTa                           <  1:454732/143‑1 (MQ=255)
                                                                                                                      aaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTAc                          >  1:205428/1‑143 (MQ=255)
                                                                                                                         gTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAAGCTGTAAGTCTGTTAGGTATAGGTTACTGt                       >  2:174220/1‑143 (MQ=255)
                                                                                                                            aaTATCCTGTCCAGCAAAGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTaaa                    <  1:457762/143‑1 (MQ=255)
                                                                                                                              tcgccTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAAtg                  >  2:205060/4‑143 (MQ=255)
                                                                                                                                 ccTGTCCTGCGATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGACTCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTgg               <  2:109230/143‑1 (MQ=255)
                                                                                                                                 ccTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTgg               >  1:770590/1‑143 (MQ=255)
                                                                                                                                  cTGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGCTTTAGGTTACTGTCAATTTGGc              >  1:142238/1‑143 (MQ=255)
                                                                                                                                   tGTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCt             >  1:477345/1‑143 (MQ=255)
                                                                                                                                          cAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGGGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAg      >  1:69440/1‑143 (MQ=255)
                                                                                                                                             tGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaa   <  1:316727/143‑1 (MQ=255)
                                                                                                                                              ggCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaaa  >  2:656476/1‑143 (MQ=255)
                                                                                                                                              ggCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGaaaa  <  1:488561/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GTTGATAACTTTCAGGAGTAACATTACCTGAAAAATATCCCACCCCAACGACGCCAAGATCAATATTACCCTGGTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAGTGGGATGCATATATAGTTCCTGAATCTGTAAGTCTGTTAGTTATAGGTTACTGTAAATGTGGCTGGTGCAGAAAA  >  NC_000913/151380‑151664

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: