Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,197,365 G→A 100% pseudogene (892/2847 nt) yehH → DUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,197,3650GA97.8% 258.7 / ‑3.8 85pseudogene (892/2847 nt)yehHDUF4132 family pseudogene; defective yehI paralog; molybdate metabolism regulator, first fragment; molybdate metabolism regulator, second fragment 2; yehI' molR mutations have phenotypes and the N‑terminal molR fragment of yehI' may be expressed and may have retained or evolved MolR function
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (47/36);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTAC  >  NC_000913/2197269‑2197461
                                                                                                |                                                                                                
tACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCg                                                                                               >  1:2082358/1‑100 (MQ=255)
 aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAGAAATGTGCACAAAAACTTATTACCGc                                                                                              >  2:2553702/1‑100 (MQ=255)
 aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                              >  1:563930/1‑100 (MQ=255)
 aCTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGc                                                                                              >  1:554566/1‑100 (MQ=255)
      gAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGAAATTc                                                                                         >  1:1069376/1‑100 (MQ=255)
       aGCAATATTCCAGAGTTGAATGACGCCTACCTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTAAAAAAAAACTTATTACCGCAAtt                                                                                          <  2:2420588/98‑1 (MQ=255)
         ccatatTCCAGATTTAAATAACGACTACGTAAACATTTATCACTGACGGAAAAGTATCTCTGGCAAATATGTGCACAAAAAATTATTACCGCAATTCca                                                                                       >  2:29281/3‑99 (MQ=255)
         cAATATTCCAGATTTGACGGACGCCTACGTAAACATTTATTACTGACGGAATAGTCTCTCTGGCAAAAATGTGCACAAAATCTTATTACCGCAATTCcan                                                                                      <  2:2229289/100‑2 (MQ=255)
          aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                     >  2:179730/1‑100 (MQ=255)
          aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                     <  1:880411/100‑1 (MQ=255)
          aaTATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCcaca                                                                                     >  1:498532/1‑100 (MQ=255)
             atTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACatat                                                                                  >  2:425323/1‑100 (MQ=255)
             atTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACatat                                                                                  >  1:1270160/1‑100 (MQ=255)
              ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAtt                                                                                 >  2:455982/1‑100 (MQ=255)
              ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATAtt                                                                                 >  1:2350638/1‑100 (MQ=255)
              ttCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATACcaca                                                                                     >  2:999594/1‑96 (MQ=255)
                ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTcc                                                                               <  2:1728822/100‑1 (MQ=255)
                ccAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTcc                                                                               >  2:2295305/1‑100 (MQ=255)
                  aGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGTAATTCCACATATTCCAg                                                                             <  2:1437967/100‑1 (MQ=255)
                     tttGAATGACGCTTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCATTTCCACATATTCCAGAAt                                                                          <  2:16652/100‑1 (MQ=255)
                     tttGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCCCATATTCCAGAAt                                                                          <  1:2040161/100‑1 (MQ=255)
                       tGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATgg                                                                        >  1:2402279/1‑100 (MQ=255)
                          aTGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATggcgg                                                                     >  1:1947304/1‑100 (MQ=255)
                            gACGCCTACGTAAACATTTATCACTTACGTAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCa                                                                   <  2:1749633/100‑1 (MQ=255)
                              cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAc                                                                 <  1:1607897/100‑1 (MQ=255)
                              cGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAAGGGCGGCa                                                                   <  2:1445755/98‑1 (MQ=255)
                                ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAAGCa                                                               >  2:1838999/1‑100 (MQ=255)
                                ccTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                               >  2:1711083/1‑100 (MQ=255)
                                  tACGTAAACATTTATCACTGACGGAAGCGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCAAATATTCCAGAATGGCGGCAACAAtt                                                             >  1:1167745/1‑100 (MQ=255)
                                  tACGTAAACATTTATCACTGACGGAAGAGTTTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCa                                                               <  1:986139/98‑1 (MQ=255)
                                  tACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCAtt                                                             >  2:2287192/1‑100 (MQ=255)
                                    cGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAATCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTaa                                                           <  2:1936388/100‑1 (MQ=255)
                                      tAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAtt                                                         >  2:474638/1‑100 (MQ=255)
                                      tAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAAtt                                                         >  1:1907075/1‑100 (MQ=255)
                                       aaaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                        <  2:934402/100‑1 (MQ=255)
                                        aaCATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTg                                                        >  2:821196/1‑99 (MQ=255)
                                           aTTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCttt                                                    >  2:1607666/1‑100 (MQ=255)
                                           aTTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCttt                                                    <  1:2070016/100‑1 (MQ=255)
                                           aTTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCttt                                                    <  1:786524/100‑1 (MQ=255)
                                            tttATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTg                                                   <  2:87076/100‑1 (MQ=255)
                                            tttATCACTGACGGAAGAGGCTCTCGTGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCTGCTACCATTAATTGCTTTg                                                   <  2:487119/100‑1 (MQ=255)
                                             ttATCACTGACGGAAGAGTCTCTCTTGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGt                                                  <  1:58827/100‑1 (MQ=255)
                                                 cACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatt                                              >  1:1598993/1‑100 (MQ=255)
                                                 cACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatt                                              >  1:1469365/1‑100 (MQ=255)
                                                  aCTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGttatta                                             >  1:1036912/1‑100 (MQ=255)
                                                    tGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTAc                                            <  2:2810042/99‑1 (MQ=255)
                                                        ggAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCAGCAACCATTAATTGCTTTGTTATTACCCGaa                                       >  1:2970618/1‑100 (MQ=255)
                                                         gAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaa                                      >  2:415419/1‑100 (MQ=255)
                                                          aatagactatctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaaa                                     <  2:2096459/91‑1 (MQ=255)
                                                          aaGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGaaaa                                     <  1:2963400/100‑1 (MQ=255)
                                                            gagTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAc                                   <  2:53027/100‑1 (MQ=255)
                                                             agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                  <  2:2691619/100‑1 (MQ=255)
                                                             agTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAAcc                                  <  1:220269/100‑1 (MQ=255)
                                                              gTCTCTCTGGCAAAAATGTGCACAAAAACTTATTACCGCAATTTCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACTCGAAAAACCa                                 >  2:664875/1‑100 (MQ=255)
                                                                ctctctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGa                               >  1:2271735/1‑100 (MQ=255)
                                                                    ctGGCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAAt                            >  1:428362/1‑99 (MQ=255)
                                                                     tGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTg                          >  2:231359/1‑100 (MQ=255)
                                                                      ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                         <  1:2378709/100‑1 (MQ=255)
                                                                      ggCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                         <  1:1741868/100‑1 (MQ=255)
                                                                       gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                        <  2:563927/100‑1 (MQ=255)
                                                                       gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                        <  2:2781127/100‑1 (MQ=255)
                                                                       gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                        <  2:1569659/100‑1 (MQ=255)
                                                                       gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGca                        >  1:2083620/1‑100 (MQ=255)
                                                                       gCAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGc                         <  1:774653/99‑1 (MQ=255)
                                                                        ctaaaaTGTGCACAAAAACTTATTATCGTATTTCCGTATAATCCAGAATGGCGGCAATCATTACTTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                       <  2:1364447/98‑1 (MQ=255)
                                                                        cAAAAATGTGCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGcac                       <  1:425326/100‑1 (MQ=255)
                                                                               gtgCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAAt                >  1:1885557/1‑100 (MQ=255)
                                                                                 gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg              >  2:1334477/1‑100 (MQ=255)
                                                                                 gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTg              <  1:2482154/100‑1 (MQ=255)
                                                                                 gCACAAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAACTGCACATGAAATTg              >  2:691355/1‑100 (MQ=255)
                                                                                  cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATTGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc             >  2:151435/1‑100 (MQ=255)
                                                                                  cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc             <  1:1304745/100‑1 (MQ=255)
                                                                                  cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGc             >  2:176574/1‑100 (MQ=255)
                                                                                  cacaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTACTTAATTAATACCCGAAAAACCAGAAATTGCACATGAAATTGc             >  1:64904/1‑100 (MQ=255)
                                                                                   acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGcc            >  2:1405780/1‑100 (MQ=255)
                                                                                   acaAAAACTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTAcc            >  1:2197693/1‑100 (MQ=255)
                                                                                      aaaaCTTATTACCGCAATTCCACATATTCCAGCATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg         >  1:468507/1‑100 (MQ=255)
                                                                                      aaaaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg         >  2:2799181/1‑100 (MQ=255)
                                                                                       aaaCTTATTACCGCAATTCCACATATGCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAg         <  2:1337878/99‑1 (MQ=255)
                                                                                        aaCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCg       <  2:1397434/100‑1 (MQ=255)
                                                                                         aCTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGt      >  1:2245451/1‑100 (MQ=255)
                                                                                          cTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTc     >  1:2289217/1‑100 (MQ=255)
                                                                                          cTTATTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTc     >  1:1377512/1‑100 (MQ=255)
                                                                                            tATTACCGCAATTCCACAAATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTTCCCAGCGTCTa   >  2:2817243/1‑100 (MQ=255)
                                                                                             aTTACCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTAc  <  1:1055670/100‑1 (MQ=255)
                                                                                                |                                                                                                
TACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGCGGCAACCATTAATTGCTTTGTTATTACCCGAAAAACCAGAAATTGCACATGAAATTGCCCAGCGTCTAC  >  NC_000913/2197269‑2197461

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