Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,615,859 T→G T403P (ACC→CCC)  yjjI ← DUF3029 family protein, putative glycine radical enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,615,8590TG97.6% 290.6 / ‑1.1 82T403P (ACC→CCC) yjjIDUF3029 family protein, putative glycine radical enzyme
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (42/38);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTT  >  NC_000913/4615762‑4615955
                                                                                                 |                                                                                                
gCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGa                                                                                                <  2:7774/100‑1 (MQ=255)
  cTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATc                                                                                              >  2:1382431/1‑100 (MQ=255)
  cTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATc                                                                                              <  1:426306/100‑1 (MQ=255)
   tGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCg                                                                                             >  2:2165774/1‑100 (MQ=255)
    gCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCgg                                                                                            <  2:385395/100‑1 (MQ=255)
       ggTTACGTTTGATGATTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAAAAATAAGCATTATGCGGGGCGACAGTTTGCAGATGGGGGATCg                                                                                             <  2:1983128/96‑1 (MQ=255)
         ttACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTg                                                                                       <  2:1395752/100‑1 (MQ=255)
          tACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTgg                                                                                      <  1:620486/100‑1 (MQ=255)
          tACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTgg                                                                                      <  2:1723438/100‑1 (MQ=255)
           aCGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGc                                                                                     <  2:174791/100‑1 (MQ=255)
             gTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTc                                                                                   <  2:713289/100‑1 (MQ=255)
              tttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCGGATGGGGGATCGGATCTGGCTCa                                                                                  >  2:1647938/1‑100 (MQ=255)
              tttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCa                                                                                  <  2:651949/100‑1 (MQ=255)
              tttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCa                                                                                  <  1:1917618/100‑1 (MQ=255)
               ttGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCAt                                                                                 >  2:493844/1‑100 (MQ=255)
                tGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATc                                                                                >  2:512643/1‑100 (MQ=255)
                  aTGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGc                                                                              <  2:647551/100‑1 (MQ=255)
                  aTGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGc                                                                              <  2:1802669/100‑1 (MQ=255)
                     gTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCAt                                                                           <  1:1289577/100‑1 (MQ=255)
                       ttCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATAc                                                                         <  2:1389121/100‑1 (MQ=255)
                            cgAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCaa                                                                    >  2:934418/1‑100 (MQ=255)
                                  tCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgc                                                              >  2:1826411/1‑100 (MQ=255)
                                   cAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAAcgcgcg                                                             >  2:1921026/1‑100 (MQ=255)
                                    aGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACgcgcgc                                                            >  1:823/1‑100 (MQ=255)
                                       aTGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCg                                                         >  1:429501/1‑100 (MQ=255)
                                         gTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGc                                                       >  2:65859/1‑100 (MQ=255)
                                            gCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                    >  2:1871083/1‑100 (MQ=255)
                                            gCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtg                                                    >  1:1630865/1‑100 (MQ=255)
                                             cTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtgg                                                   >  1:714962/1‑100 (MQ=255)
                                             cTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtgg                                                   >  2:1094711/1‑100 (MQ=255)
                                              tGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggt                                                  >  1:1216523/1‑100 (MQ=255)
                                              tGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggt                                                  <  2:429499/100‑1 (MQ=255)
                                               gATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCgtggtg                                                 >  1:1438205/1‑100 (MQ=255)
                                                  gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCg                                              >  1:1024868/1‑100 (MQ=255)
                                                  gCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTCGCAGATGGGGGATCGGATCTGGCTCATCGCCCTACGGCAAACGCGCGCCCGGCGTGGTGCCg                                              >  2:1965450/1‑100 (MQ=255)
                                                   ccGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGa                                             <  1:934420/100‑1 (MQ=255)
                                                    cGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGat                                            <  1:1526997/100‑1 (MQ=255)
                                                    cGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGat                                            >  2:2336213/1‑100 (MQ=255)
                                                      gAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGatat                                          >  2:1608830/1‑100 (MQ=255)
                                                        ataataataAGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCg                                        >  1:819946/1‑100 (MQ=255)
                                                         taataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                       <  2:368058/100‑1 (MQ=255)
                                                         taataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCgg                                       <  2:431297/100‑1 (MQ=255)
                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                     >  2:117314/1‑100 (MQ=255)
                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                     >  1:298748/1‑100 (MQ=255)
                                                           ataataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGaa                                     >  1:2256466/1‑100 (MQ=255)
                                                            taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                    <  1:1382434/100‑1 (MQ=255)
                                                            taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                    >  2:1652950/1‑100 (MQ=255)
                                                            taataaGCATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAAc                                    <  1:512645/100‑1 (MQ=255)
                                                                   cATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATccc                             >  1:436465/1‑100 (MQ=255)
                                                                   cATGATGCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATccc                             >  1:1627991/1‑100 (MQ=255)
                                                                    atgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccc                            <  2:1622672/100‑1 (MQ=255)
                                                                    atgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATcccc                            <  1:1608835/100‑1 (MQ=255)
                                                                     tgatgCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                           >  1:650996/1‑100 (MQ=255)
                                                                     tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTgcccgcct                           >  2:337877/1‑93 (MQ=255)
                                                                     tgatgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCg                           >  2:1631574/1‑100 (MQ=255)
                                                                       atgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGAc                         >  2:309622/1‑100 (MQ=255)
                                                                        tgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACt                        <  2:1958819/100‑1 (MQ=255)
                                                                        tgCGGGGCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACt                        >  2:1477333/1‑100 (MQ=255)
                                                                           ggggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgtg                     >  1:364758/1‑100 (MQ=255)
                                                                           ggggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACtgtg                     <  2:1736628/100‑1 (MQ=255)
                                                                             ggCGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCg                   <  1:395668/100‑1 (MQ=255)
                                                                               cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                 <  1:1384883/100‑1 (MQ=255)
                                                                               cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                 >  1:827290/1‑100 (MQ=255)
                                                                               cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                 <  2:1403683/100‑1 (MQ=255)
                                                                               cGACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                 <  2:2256458/100‑1 (MQ=255)
                                                                               cGACAGTTTGCAGATGGGGGATCGGATCTAGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtg                 <  2:1405904/100‑1 (MQ=255)
                                                                                gACAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCgtgt                <  2:1438201/100‑1 (MQ=255)
                                                                                 aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa               <  2:2166616/100‑1 (MQ=255)
                                                                                 aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa               <  1:309624/100‑1 (MQ=255)
                                                                                 aCAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTa               <  1:652275/100‑1 (MQ=255)
                                                                                  cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa              >  2:1499007/1‑100 (MQ=255)
                                                                                  cAGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTaa              >  1:936828/1‑100 (MQ=255)
                                                                                   aGTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAAc             <  1:729439/100‑1 (MQ=255)
                                                                                    gTTTGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACa            >  2:2311963/1‑100 (MQ=255)
                                                                                      ttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACAt           >  2:587482/1‑99 (MQ=255)
                                                                                      ttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATg          <  2:1767400/100‑1 (MQ=255)
                                                                                      ttGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATg          <  2:818569/100‑1 (MQ=255)
                                                                                       tGCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATgg         <  2:823/100‑1 (MQ=255)
                                                                                        gCAGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgc        >  1:1971574/1‑100 (MQ=255)
                                                                                          aGATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGgcgc      >  2:568896/1‑100 (MQ=255)
                                                                                           gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGcgcg     >  1:2160226/1‑100 (MQ=255)
                                                                                           gATGGGGGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATAACGGAACTGATCACCGACTGTGCGTAAAACATGgcga      >  1:429872/1‑98 (MQ=255)
                                                                                              ggggGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGttt  <  2:429261/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GCCTGCGGGTTACGTTTGATGGTTTCGTCGAGCGTCAGAATGTCGCTGATGCCGGAATAATAATAAGCATGATGCGGGGCGACAGTTTGCAGATGGGTGATCGGATCTGGCTCATCGCCATACGGCAAACGCGCGCCCGGCGTGGTGCCGATATCGGAACTGATCCCCGACTGTGCGTGTAACATGGCGCGTTT  >  NC_000913/4615762‑4615955

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