Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,129,394 G→C 18.6% intergenic (+20/+20) murJ → / ← flgN putative lipid II flippase/export chaperone for FlgK and FlgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,129,3940GC18.6% 172.6 / 36.8 70intergenic (+20/+20)murJ/flgNputative lipid II flippase/export chaperone for FlgK and FlgL
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base C (6/7);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.89e-01

GGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGC  >  NC_000913/1129296‑1129492
                                                                                                  |                                                                                                  
ggCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagt                                                                                                   <  2:744205/100‑1 (MQ=255)
ggCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagt                                                                                                   >  2:2396545/1‑100 (MQ=255)
  cGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgc                                                                                                 <  1:1181349/100‑1 (MQ=255)
  cGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgc                                                                                                 <  1:975444/100‑1 (MQ=255)
  cGGGAATTGCCGCGTACTTCGCTGCTCTGGCGGTACTGGGCTTCACAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGTATTCCGGCCtgcagtgc                                                                                                 <  2:1270819/100‑1 (MQ=255)
    gtgATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCtgcagtgcag                                                                                               <  1:1322106/98‑1 (MQ=255)
      gATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGc                                                                                             <  2:366257/100‑1 (MQ=255)
        ttGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCg                                                                                           <  1:319880/100‑1 (MQ=255)
         ttccgatctcTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCgg                                                                                          <  1:1765111/91‑1 (MQ=255)
         tGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCgg                                                                                          >  2:1383241/1‑100 (MQ=255)
         tGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCgg                                                                                          <  1:1731231/100‑1 (MQ=255)
         tGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCgg                                                                                          <  2:1377879/100‑1 (MQ=255)
          gCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                         <  1:983295/100‑1 (MQ=255)
          gCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGga                                                                                         <  1:1285106/100‑1 (MQ=255)
          gCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGga                                                                                         <  1:514205/100‑1 (MQ=255)
           ccGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGgaa                                                                                        >  1:1636988/1‑99 (MQ=255)
             gcgTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGAt                                                                                      >  2:242155/1‑100 (MQ=255)
             gcgTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGgaatg                                                                                      <  1:742787/100‑4 (MQ=255)
                tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGga                                                                                         <  1:1408606/94‑1 (MQ=255)
                tACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAAt                                                                                   >  2:1207920/1‑100 (MQ=255)
                  cTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCACTGCAGGCCGGAGAtcggaa                                                                                 >  2:1765090/1‑95 (MQ=255)
                  cTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGACTGCAGTGCAGGTCGGAGATAAt                                                                                   >  2:2234708/1‑98 (MQ=255)
                  cTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCAGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCt                                                                                 <  1:1606102/100‑1 (MQ=255)
                    tCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTc                                                                               <  1:2268623/100‑1 (MQ=255)
                     cGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCa                                                                              <  2:1693576/100‑1 (MQ=255)
                       cTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGa                                                                            >  1:2163524/1‑100 (MQ=255)
                         gCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTATCAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATc                                                                          <  2:1330843/100‑1 (MQ=255)
                         gCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATc                                                                          <  2:1788111/100‑1 (MQ=255)
                          cACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCg                                                                         >  1:736813/1‑100 (MQ=255)
                              ggCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAAt                                                                     <  1:1192088/100‑1 (MQ=255)
                               gCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATc                                                                    <  2:1077437/100‑1 (MQ=255)
                                   tACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                <  1:871509/100‑1 (MQ=255)
                                   tACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCtttt                                                                >  2:1803675/1‑100 (MQ=255)
                                      tGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTAc                                                             >  1:600561/1‑100 (MQ=255)
                                       gggCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAAGCTTTTTAcc                                                            >  2:1630150/1‑100 (MQ=255)
                                          cTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTAccgcc                                                         >  1:1319873/1‑100 (MQ=255)
                                          cTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTAccgcc                                                         >  1:1629092/1‑100 (MQ=255)
                                            tCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACcgccgc                                                       >  1:1063404/1‑100 (MQ=255)
                                                aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATg                                                   <  2:253904/100‑1 (MQ=255)
                                                aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATg                                                   >  1:1735364/1‑100 (MQ=255)
                                                aGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATg                                                   >  1:128306/1‑100 (MQ=255)
                                                 gTTAAAGAATTTGCCCGCCGGACGATGGAACAATGCATTCCGGCCTGCAGGGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATg                                                   <  2:169/99‑1 (MQ=255)
                                                   tAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgt                                                <  1:505778/100‑1 (MQ=255)
                                                    aaaGAATTTGCCCGCCGGACGGAGAAACAATGCATTCCGGCCTGCAGTGCAGGCCGCAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATggg                                                 >  2:1534107/1‑98 (MQ=255)
                                                     aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgttg                                              <  2:1797348/100‑1 (MQ=255)
                                                     aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgttg                                              >  1:2238196/1‑100 (MQ=255)
                                                       gAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGgttgtt                                            >  1:902416/1‑100 (MQ=255)
                                                           ttGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAgg                                        <  2:1248452/100‑1 (MQ=255)
                                                             tcctgccaGACGGTGTAACAATTCATTCCGGCCTACAGTGCAGGCCGGAGAGAATCTTCAGATCGAAATCTTTTTACCTCCGCGATGGATTGTTGAg                                         >  2:946239/5‑97 (MQ=255)
                                                             gCCCGCCGGACGGTGTAACAATGCATTCCTGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTc                                      <  2:2034349/100‑1 (MQ=255)
                                                                cgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGa                                   >  1:1887994/1‑100 (MQ=255)
                                                                cgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGa                                   >  1:2104662/1‑100 (MQ=255)
                                                                 gccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAc                                  <  1:90995/100‑1 (MQ=255)
                                                                 gccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAc                                  <  2:736805/100‑1 (MQ=255)
                                                                  ccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAcc                                 <  1:593184/100‑1 (MQ=255)
                                                                   cgGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCg                                <  1:1399155/100‑1 (MQ=255)
                                                                    ggATGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCGTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTGGAGGTCTGACCg                                <  2:526085/99‑1 (MQ=255)
                                                                    ggACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGAcc                                 <  1:722169/98‑1 (MQ=255)
                                                                       cGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCg                            <  1:1639484/100‑1 (MQ=255)
                                                                         gtgtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGcc                          >  2:1249926/1‑100 (MQ=255)
                                                                         gtgtAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGcc                          >  2:2298343/1‑100 (MQ=255)
                                                                            tAACAGTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCa                       <  1:867060/100‑1 (MQ=255)
                                                                            tAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCa                       <  2:1242351/100‑1 (MQ=255)
                                                                            tAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCa                       <  1:2350424/100‑1 (MQ=255)
                                                                            tAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTAAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACTGTTCGCCCCa                       >  2:1932153/1‑100 (MQ=255)
                                                                                 aTGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAg                  >  1:2383702/1‑100 (MQ=255)
                                                                                    cATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGt               >  1:1280237/1‑100 (MQ=255)
                                                                                     tctCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTc              >  1:1693852/3‑100 (MQ=255)
                                                                                     tctCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAAATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGt               >  1:1693866/3‑99 (MQ=255)
                                                                                      ttCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCg             >  1:963161/1‑100 (MQ=255)
                                                                                       tCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCgg            >  2:809855/1‑100 (MQ=255)
                                                                                       tCCGGCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCg             >  1:2102532/1‑99 (MQ=255)
                                                                                           gCCTGCACTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTcc        >  1:41236/1‑100 (MQ=255)
                                                                                            ccTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt       <  1:1719648/100‑1 (MQ=255)
                                                                                             cTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCt       <  1:1508735/99‑1 (MQ=255)
                                                                                             cTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTg      >  1:539699/1‑100 (MQ=255)
                                                                                               gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGAt    >  2:917408/1‑100 (MQ=255)
                                                                                               gcagtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGAt    >  1:1804740/1‑100 (MQ=255)
                                                                                                 agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGc  <  1:620066/100‑1 (MQ=255)
                                                                                                 agtgcagGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGc  <  1:1960984/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
GGCGGGGATTGCCGCGTACTTCGCTGCACTGGCGGTACTGGGCTTCAAAGTTAAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGGTTGTTGAGGTCTGACCGTTCGCCCCATATAGCGTCGGTTCCTGATGC  >  NC_000913/1129296‑1129492

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