Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,509,471 G→A intergenic (‑45/‑159) glk ← / → yfeO glucokinase/putative ion channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,509,4710GA98.4% 181.4 / NA 62intergenic (‑45/‑159)glk/yfeOglucokinase/putative ion channel protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (34/27);  minor base C (0/1);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGAC  >  NC_000913/2509332‑2509610
                                                                                                                                           |                                                                                                                                           
gccTTAGCCCGCGCGTTTGCACCACTGGCAATACCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATaa                                                                                                                                          <  2:573421/141‑1 (MQ=255)
gTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATaa                                                                                                                                          >  1:611004/1‑143 (MQ=255)
 tCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATaaa                                                                                                                                         >  1:481833/1‑143 (MQ=255)
      gCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTTCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCt                                                                                                                                    <  2:581587/143‑1 (MQ=255)
       ccTGCGAGAATTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGAATACTTGGTCATTCTTCAACTGCTCCGCTACAGTCAAAATAATTATTTCTCAAACCATAAATGccgg                                                                                                                                   >  2:627350/1‑141 (MQ=255)
        cTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCt                                                                                                                                    <  2:433741/141‑1 (MQ=255)
        cTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCt                                                                                                                                    >  1:433741/1‑141 (MQ=255)
         tGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTggg                                                                                                                                 <  1:943122/143‑1 (MQ=255)
           cGAGATTGCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCACTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTggggg                                                                                                                               <  1:932607/143‑1 (MQ=255)
            gagaTTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACATGGGGGc                                                                                                                              >  1:427469/1‑143 (MQ=255)
                tttCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATaa                                                                                                                          >  2:837915/1‑143 (MQ=255)
                 ttCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATaac                                                                                                                         <  1:713945/143‑1 (MQ=255)
                  tCACCACTGGCAATATCACACCGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATaaca                                                                                                                        <  2:539017/143‑1 (MQ=255)
                               tatCACACAGAGCAAGACGTGCGTTGAAGCCCCCCAAATCACTGACTAGTGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGaaa                                                                                                           >  1:835752/1‑143 (MQ=255)
                                    cacaGAGCAAGGCGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAaggcag                                                                                                      <  1:221786/143‑1 (MQ=255)
                                    cacaGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAaggcaa                                                                                                      >  2:1040569/1‑142 (MQ=255)
                                    cacaGAGCAAGACGTGCGTTGGTGCCGCACACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCTGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACATGGGGGCATAACAACAACGCCGGaacgacac                                                                                                      >  2:450877/1‑137 (MQ=255)
                                       agagCAAGACGTGCGTTGGTGCCGCCCACATCACCGACGAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATACTTCTTTCCCACACAATAAACACCTCGGCGCATACCAACACCCCCGGAAAGCCAGCCt                                                                                                   >  2:252223/1‑143 (MQ=255)
                                         agCAAGACGTGCCTTGGTGCCCTCCACATCACCGACTAATGCATACTGTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTcc                                                                                                 <  2:899522/143‑1 (MQ=255)
                                             aGACGTGCGTTGCTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGt                                                                                             <  2:712418/143‑1 (MQ=255)
                                                      ttGGTGCCGCCCGCATCAGCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGa                                                                                    <  2:1044123/143‑1 (MQ=255)
                                                         gTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGTGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCt                                                                                 <  2:125795/143‑1 (MQ=255)
                                                                ccACATCACCGACTAATGCATAGTTTGTCTGTTTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATAACGATCTGGGTcac                                                                          >  2:647993/1‑143 (MQ=255)
                                                                      cACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAAtt                                                                    >  2:436330/1‑143 (MQ=255)
                                                                          ggcTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATGCGTTCGCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACtt                                                                <  2:84624/141‑1 (MQ=255)
                                                                          gACTAATGCATACTTTGTCCTTCTAGGACGACGCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACtt                                                                <  1:946188/143‑1 (MQ=255)
                                                                          gACTAATGCATACTTTGTCATTCTTCAACTGCTCCTCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCTTCTATTCATTCGTGATCGATCTGGGTCACACAATTACtt                                                                <  2:587064/143‑1 (MQ=255)
                                                                          gACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACtt                                                                <  2:848091/143‑1 (MQ=255)
                                                                               aTGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATAGCCGGTGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATGTGGGTCACGGCATTACGATATCg                                                           >  2:231144/1‑143 (MQ=255)
                                                                               aTGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCg                                                           >  1:739520/1‑143 (MQ=255)
                                                                                 gCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGtt                                                         <  1:736311/143‑1 (MQ=255)
                                                                                   aTACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTGTCACACTATAAATACCTGCGGGCTTAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTc                                                       >  1:132450/1‑143 (MQ=255)
                                                                                     aCTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAACGGCAGGCTCCCGGTAAATGTCGATCTGGGTCACGCAATTACTTTATCGGTTCAg                                                     >  1:157280/1‑143 (MQ=255)
                                                                                             aTTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGACAGGCAGGCTCCCTGTACATATCGATCCGGGTCACACAATTACGTTCTCGTTTCAGCACCAAtt                                             >  2:833640/1‑143 (MQ=255)
                                                                                             aTTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAACGGCAGGCTCCCTGCCCCTCTCGCTCTGGCTCCCACACTTACTTTTTCGTTCCATCCCCAAtt                                             >  1:405760/1‑143 (MQ=255)
                                                                                             aTTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAAtt                                             >  1:175461/1‑143 (MQ=255)
                                                                                             aTTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATATGGGTCACACAATTACTGCACGGTTTCAGCATCAAtt                                             >  2:610454/1‑143 (MQ=255)
                                                                                              ttcttcAACCGCTCCGCCAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTg                                            <  1:70332/143‑1 (MQ=255)
                                                                                                 ttcAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAg                                         >  1:795151/1‑143 (MQ=255)
                                                                                                  tcAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGc                                        >  1:455726/1‑143 (MQ=255)
                                                                                                  tcAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGc                                        >  1:220936/1‑143 (MQ=255)
                                                                                                       tGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGCACATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGc                                   >  2:157166/1‑143 (MQ=255)
                                                                                                          tCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATCGCAGCGATGCCtt                                >  1:458943/1‑143 (MQ=255)
                                                                                                          tCCGCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACAACGCAGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTAAAGTTTCAGCACAAATTGCAGCGCTGCCtt                                >  1:228705/1‑143 (MQ=255)
                                                                                                             tcTACAGACAAGATAATTCCTTCTCACACTCTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCttttt                             <  2:786214/142‑1 (MQ=255)
                                                                                                             gCTAAAGTCAAAATAATTCTTTCTCACACTATAAATACCTGGGGGCATAACAACCAGGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACTATTACTTTATCGGTTCAGCACCAATTGCAGCGATGCCttttt                             <  2:122285/143‑1 (MQ=255)
                                                                                                                      aaaaTAACTCTTTCTCACACTATAAATACCTGGGGGCATAACAACACCGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTgg                    <  1:853230/143‑1 (MQ=255)
                                                                                                                              tctttcgCACACTATAGATACCTGGGATCCTTACAACAACGCTGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTa            <  1:441161/143‑1 (MQ=255)
                                                                                                                              tctctctcACACAATCAATACCCGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATTGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTa            <  1:627350/143‑1 (MQ=255)
                                                                                                                               ctttctCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTGATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTaa           <  2:13409/143‑1 (MQ=255)
                                                                                                                               ctgtctcACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTaa           <  2:611004/140‑1 (MQ=255)
                                                                                                                                tttctCACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAg          >  1:790706/1‑143 (MQ=255)
                                                                                                                                tttccgacacTATAAATCCCTGGCGGCATAACAACAACGCCGGAAAGGCATGCTCCCTGTAAATATCGAACTGGGTCACACAATCACTTTATCGTTTAAGCACCAATTGCAGCGATGCCTTCTTGCACGCTGGGCAAACTaac          >  2:185859/7‑142 (MQ=255)
                                                                                                                                   ctcACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCAGCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTAt       >  2:336622/1‑143 (MQ=255)
                                                                                                                                    tcACACTATACATACCTCGGCGCATAACAACAACGCCGGAAAGCCAGGCTCCCTGTAAATACCGACCTGGGTCACACAATGACTCTGCCGCTTCAGCAGCACGTGCAGCGAAGCCCGTGTGCACGCCGGGCCAACTAAGTGTc      >  2:1062529/1‑143 (MQ=255)
                                                                                                                                    tcACACTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCACACTAAGTATc      >  1:167060/1‑143 (MQ=255)
                                                                                                                                     cacaCTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCt     <  1:845367/143‑1 (MQ=255)
                                                                                                                                     cacaCTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATAGTGCTGGTTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCt     <  1:776729/143‑1 (MQ=255)
                                                                                                                                     cacaCGATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCt     <  1:501803/143‑1 (MQ=255)
                                                                                                                                      acacTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTg    >  1:589572/1‑143 (MQ=255)
                                                                                                                                      acacTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATAGCGATCTGGGCCACACACTTGCGTCATCGTTGCAGCGCCAATTGCAGCGAGACCTTTTTGCACGCTGGGCAAACTAAGTATcag    >  2:251052/1‑141 (MQ=255)
                                                                                                                                        acTATAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGAc  >  1:705068/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GTCTTAGCCTGCGAGATTTCACCACTGGCAATATCACACAGAGCAAGACGTGCGTTGGTGCCGCCCACATCACCGACTAATGCATACTTTGTCATTCTTCAACTGCTCCGCTAAAGTCAAAATAATTCTTTCTCACACTGTAAATACCTGGGGGCATAACAACAACGCCGGAAAGGCAGGCTCCCTGTAAATATCGATCTGGGTCACACAATTACTTTATCGTTTCAGCACCAATTGCAGCGATGCCTTTTTGCAAGCTGGGCAAACTAAGTATCTGAC  >  NC_000913/2509332‑2509610

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