Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,128,197 T→C 68.8% T67A (ACC→GCC)  metJ ← transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,128,1970TC68.8% 101.7 / 58.0 77T67A (ACC→GCC) metJtranscriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (15/9);  new base C (18/35);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 2.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTA  >  NC_000913/4128098‑4128295
                                                                                                   |                                                                                                  
ggTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGt                                                                                                    >  2:2171782/1‑100 (MQ=255)
   tAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCaaa                                                                                                 <  1:1221299/100‑1 (MQ=255)
     aTCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAgg                                                                                               >  1:1264271/1‑100 (MQ=255)
      tCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGc                                                                                              <  2:1031115/100‑1 (MQ=255)
      tCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGc                                                                                              >  2:1340715/1‑100 (MQ=255)
       ccccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCa                                                                                             >  2:2048005/1‑100 (MQ=255)
       ccccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCGGCGCTCTTGACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAgg                                                                                               <  2:1868517/98‑1 (MQ=255)
         ccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATg                                                                                           <  2:1463352/100‑1 (MQ=255)
          cATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGc                                                                                          <  1:21356/100‑1 (MQ=255)
             ctcACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGa                                                                                       <  1:1434482/100‑1 (MQ=255)
               cACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGaaa                                                                                     >  2:1087458/1‑100 (MQ=255)
                   cATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCt                                                                                 >  1:1544849/1‑100 (MQ=255)
                    aTGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCtt                                                                                <  1:2071965/100‑1 (MQ=255)
                     tGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTc                                                                               <  1:1737593/100‑1 (MQ=255)
                      gATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCg                                                                              >  1:1216204/1‑100 (MQ=255)
                      gATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCg                                                                              <  1:1340721/100‑1 (MQ=255)
                       aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGc                                                                             <  2:664773/100‑1 (MQ=255)
                       aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGc                                                                             <  1:1286477/100‑1 (MQ=255)
                          tctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGcaca                                                                          <  2:32945/100‑1 (MQ=255)
                           ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCAcag                                                                         >  1:960187/1‑100 (MQ=255)
                            ttttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCAcagc                                                                        <  1:1399675/100‑1 (MQ=255)
                             tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCAcagca                                                                       >  2:2232207/1‑100 (MQ=255)
                               tGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACagcagc                                                                     >  1:1524729/1‑100 (MQ=255)
                                 ccGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTc                                                                   <  1:688020/100‑1 (MQ=255)
                                 ccGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTc                                                                   <  1:1369669/100‑1 (MQ=255)
                                  cGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCg                                                                  <  2:1010092/100‑1 (MQ=255)
                                  cGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCg                                                                  >  1:1256560/1‑100 (MQ=255)
                                   gCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGc                                                                 <  2:423335/100‑1 (MQ=255)
                                     ttCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTg                                                               <  1:1863183/100‑1 (MQ=255)
                                       ccGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGtt                                                             <  2:485631/100‑1 (MQ=255)
                                          ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTgg                                                           <  1:1054006/99‑1 (MQ=255)
                                          ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGt                                                          <  1:629403/100‑1 (MQ=255)
                                          ggATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGt                                                          <  1:988785/100‑1 (MQ=255)
                                             tttCGTCGCTGCGCTCTTTGCGCAGATCGTTATCTTCCGGCAAAGGTTGCCGGGCAAAGGCATGCAGAAACGCTTCGTATAGCAGCTCGCTGtgt                                                            <  2:816065/95‑3 (MQ=255)
                                             tttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGc                                                       >  1:1453585/1‑100 (MQ=255)
                                               tcgtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGt                                                     <  1:1691566/100‑1 (MQ=255)
                                                     cTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCa                                               <  2:228375/100‑1 (MQ=255)
                                                      tGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAg                                              >  1:542886/1‑100 (MQ=255)
                                                        cgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgt                                            <  2:217129/100‑1 (MQ=255)
                                                           tctGTACGCAGATCCGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCAGGCAGCAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgt                                            <  2:2153921/97‑1 (MQ=255)
                                                            ctTTACGAAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgttgtt                                        <  1:814856/100‑1 (MQ=255)
                                                             tttACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTc                                       <  1:1182344/100‑1 (MQ=255)
                                                               tACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAc                                     >  1:76991/1‑100 (MQ=255)
                                                                aCGCAGATCGGCATCTTCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAcc                                    <  1:2215267/100‑1 (MQ=255)
                                                                aCGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAcc                                    <  1:285671/100‑1 (MQ=255)
                                                                aCGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAcc                                    <  2:1720165/100‑1 (MQ=255)
                                                                aCGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAcc                                    <  1:1786846/100‑1 (MQ=255)
                                                                 cGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAcc                                    >  1:1730377/1‑99 (MQ=255)
                                                                 cGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCt                                   >  1:1881237/1‑100 (MQ=255)
                                                                   cAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCga                                 >  2:26634/1‑100 (MQ=255)
                                                                    aGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAcaaaa                                 >  1:855661/1‑95 (MQ=255)
                                                                    aGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                >  2:2221601/1‑100 (MQ=255)
                                                                    aGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                <  1:1770575/100‑1 (MQ=255)
                                                                    aGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                <  1:744094/100‑1 (MQ=255)
                                                                     gATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacg                               <  2:1522826/100‑1 (MQ=255)
                                                                     gATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacg                               >  1:1463131/1‑100 (MQ=255)
                                                                       tCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacgac                             >  2:1299957/1‑100 (MQ=255)
                                                                         ggCATCATCCGGCAAAGGTTGACCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGAGTGACGCAGATTGTTCACCCGACGAcgc                           >  2:1894957/1‑100 (MQ=255)
                                                                          gCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGAcgcg                          <  1:1041121/100‑1 (MQ=255)
                                                                            atcatcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtg                        <  2:397056/100‑1 (MQ=255)
                                                                              catcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgcg                      >  2:339712/1‑100 (MQ=255)
                                                                               atcCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACgcgtgcgt                     >  1:1436211/1‑100 (MQ=255)
                                                                               atcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgcgt                     >  2:1241303/1‑100 (MQ=255)
                                                                                  cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCa                  >  1:592122/1‑100 (MQ=255)
                                                                                  cGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCa                  >  2:34872/1‑100 (MQ=255)
                                                                                    gCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATc                <  2:858219/100‑1 (MQ=255)
                                                                                    gCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATc                <  1:719444/100‑1 (MQ=255)
                                                                                      aaaGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCCCGCTGTTGGTAGCGTTACGCAGGTTGTTAACCTGACGACACGTGCGTTCATCgg              >  1:33105/1‑100 (MQ=255)
                                                                                       aagtttttccGGCAATGGCTTGCAGAAACGCTTCGCACAGCAGCTTGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCg               <  1:45540/90‑1 (MQ=255)
                                                                                        aGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTg            <  1:1299963/100‑1 (MQ=255)
                                                                                           ttGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgg         >  2:337301/1‑100 (MQ=255)
                                                                                            tGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCCGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgg         >  2:1267614/1‑99 (MQ=255)
                                                                                             gCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTCCAACCGACGACGCGTGCGTTCATCGGTGAGGAt       <  2:1475098/100‑1 (MQ=255)
                                                                                              cccGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCCCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtt      >  2:12292/1‑100 (MQ=255)
                                                                                               ccGGCACAGGCATGCAGAAACGCTTCGCACAACAGCTAGCTGTTGGTAGCGTGACGAAGGTTCTTACCCTGAAGACCCGTGCGTTCATCGGTGa           >  2:714305/1‑94 (MQ=255)
                                                                                                cGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGAtttt    >  1:1297098/1‑100 (MQ=255)
                                                                                                 ggTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGGTCATCGGTGAGGAttttt   >  1:1571187/1‑100 (MQ=255)
                                                                                                  gTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTa  <  1:1378730/100‑1 (MQ=255)
                                                                                                   tAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCCGACGACGCGTGCGTTCATCGGTGAGGATTTTTa  <  1:1953686/99‑1 (MQ=255)
                                                                                                   |                                                                                                  
GGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTA  >  NC_000913/4128098‑4128295

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