Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,236,177 Δ1 bp intergenic (+239/‑45) ygjR → / → ygjT hypothetical protein/putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,236,1760G.98.4% 525.8 / ‑0.7 121intergenic (+238/‑46)ygjR/ygjThypothetical protein/putative transport protein
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base . (67/52);  total (67/54)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

ACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGC  >  NC_000913/3236101‑3236246
                                                                           |                                                                      
aCGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTg                                                                        >  1:5299433/1‑76 (MQ=255)
 cGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTga                                                                       <  2:3102856/76‑2 (MQ=255)
  gCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTgac                                                                      >  2:3817194/1‑74 (MQ=255)
  gCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTgac                                                                      <  1:4242914/76‑3 (MQ=255)
  gCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTACCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTgac                                                                      <  1:2570193/76‑3 (MQ=255)
   cAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTgacg                                                                     >  2:3292491/1‑73 (MQ=255)
     gACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTc                                                                  <  1:3902625/76‑1 (MQ=255)
      aCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCt                                                                 >  1:3324067/1‑76 (MQ=255)
      aCGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCt                                                                 >  1:3037280/1‑76 (MQ=255)
       cGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCtt                                                                >  1:1288386/1‑76 (MQ=255)
         tATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcc                                                              <  1:3258760/76‑1 (MQ=255)
         tATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcc                                                              >  2:4298147/1‑76 (MQ=255)
         tATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcc                                                              >  1:3060065/1‑76 (MQ=255)
         tATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcc                                                              <  2:206612/76‑1 (MQ=255)
          aTTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCg                                                             <  2:5480238/76‑1 (MQ=255)
          aTTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCg                                                             <  2:3897604/76‑1 (MQ=255)
           ntCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTcccc                                                            >  2:1351947/2‑74 (MQ=255)
            tCCTTTTTTGTTGTCAGTGCGACCTTGCCGGAAGGCTAGTTCTATGCATAAACCGCATCGGCTGGCGTCTTCCgcc                                                           <  1:2446087/76‑3 (MQ=255)
            tCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATTCATAAAAAGCAGCGGCTGACGTCTTCCGAc                                                           <  1:720977/76‑1 (MQ=255)
            tCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGAc                                                           <  2:1812433/76‑1 (MQ=255)
            tCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGAc                                                           <  2:4647222/76‑1 (MQ=255)
             ccTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCCACg                                                          >  2:2168181/1‑76 (MQ=255)
               ttttttGTTGTAGGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGtt                                                        <  2:846715/76‑1 (MQ=255)
               ttttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAGAAAGCAGCGGCTGACGTCTTCCGACGtt                                                        <  1:1134515/76‑1 (MQ=255)
               ttttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGtt                                                        <  2:1811091/76‑1 (MQ=255)
               ttttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGtt                                                        >  1:5609835/1‑76 (MQ=255)
               ttttttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTATATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGtt                                                        >  1:5666152/1‑76 (MQ=255)
                 tcttgttgTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTgg                                                      <  1:1414520/74‑1 (MQ=255)
                  tttGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGc                                                     >  1:4556517/1‑76 (MQ=255)
                   ttgttgTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGcc                                                    <  1:3095473/76‑1 (MQ=255)
                     gttgtAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGt                                                  >  2:4266518/1‑76 (MQ=255)
                      ttgtAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtt                                                 >  1:4909058/1‑76 (MQ=255)
                      ttgtAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtt                                                 <  1:1563377/76‑1 (MQ=255)
                      ttgtAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACCTTGGCCggt                                                 >  2:5084761/1‑74 (MQ=255)
                       tgtAAGTGAGACCTTGCCGTAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttt                                                <  1:5015473/76‑1 (MQ=255)
                       tgtAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttt                                                <  2:1456641/76‑1 (MQ=255)
                       tgtAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttt                                                >  1:3866576/1‑76 (MQ=255)
                          aaGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttttt                                             <  1:4880793/76‑1 (MQ=255)
                          aaGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGGTGGCCGgttttt                                             >  1:1702642/1‑76 (MQ=255)
                           aGTGAGGCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGttttttt                                            >  1:1224912/1‑76 (MQ=255)
                            gTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGtttttttt                                           >  1:3494129/1‑75 (MQ=255)
                             tGAGACGTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAt                                          <  2:3626301/76‑1 (MQ=255)
                             tGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAt                                          >  1:4217994/1‑76 (MQ=255)
                              gagaCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtt                                         >  2:2275703/1‑75 (MQ=255)
                              gagaCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtg                                         <  1:4571553/76‑1 (MQ=255)
                              gagaCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtg                                         >  2:2364110/1‑76 (MQ=255)
                               agcccTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgt                                        <  2:74329/73‑1 (MQ=255)
                                gcccTTGCGGGAAGTCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgtg                                       <  2:3350450/74‑1 (MQ=255)
                                gaCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgtg                                       >  1:1610246/1‑76 (MQ=255)
                                gaCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTAtgtg                                       >  2:4438188/1‑76 (MQ=255)
                                 aCCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGGTTTTTTATgtgt                                      >  2:332609/1‑76 (MQ=255)
                                  ccTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTa                                     <  2:2324606/76‑1 (MQ=255)
                                   cTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTaa                                    <  1:2289908/76‑1 (MQ=255)
                                   cTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTaa                                    >  2:3857811/1‑76 (MQ=255)
                                    ttGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTCTAAg                                   >  2:3533310/1‑76 (MQ=255)
                                     tGCCGGCAGGCGAGGTCTAGGCATAAAAAGCAGCGGCTGACGTCTTCCGCCGTTGGCCGTTTTTTTATGTGTAAgg                                  <  2:5580565/76‑1 (MQ=255)
                                     tGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAgg                                  >  1:1880008/1‑76 (MQ=255)
                                        cGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAAc                               >  2:3942815/1‑76 (MQ=255)
                                         ggAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACt                              <  1:2973978/76‑1 (MQ=255)
                                           aaGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTc                            >  1:3217298/1‑76 (MQ=255)
                                           aaGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTc                            >  1:2340446/1‑76 (MQ=255)
                                            cggCGAGGTCGATGCATAAAAAGCAGCGACTTACGTCTTCCGACGCTCGCCGTTTTTTCATTTTGCAGGAACctct                           >  2:4756375/2‑76 (MQ=255)
                                            aGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATTTGTAAGGAAATTCt                           >  2:4753697/1‑76 (MQ=255)
                                            aGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCt                           <  2:1731098/76‑1 (MQ=255)
                                            aGGCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCt                           >  2:3045600/1‑76 (MQ=255)
                                             ggCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTa                          <  2:5461505/76‑1 (MQ=255)
                                             ggCGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTa                          <  1:4266519/76‑1 (MQ=255)
                                             ggCGAGGTCTATGCATAAAAAGCAGCGGCTGAAGTCATAAGCCGTCGGCCGTATTTTTCTGTGTAAGGCACTTGTa                          >  2:2982501/1‑76 (MQ=255)
                                             ggCGAGGTCTATGCATAAAAAGCAGCAGCTGACGTCATCCGACGTCGCACGCTTTTTTATCTGTAAGGAACTCCtt                          >  2:5535597/1‑75 (MQ=255)
                                              gagAGGTCTATGCATCAAAAGCAGCCGCTGACGTCTTCCGACGTTGGCCGTTTTCTTCTGTGTCAGGAACTTCTAt                         >  1:3686394/3‑76 (MQ=255)
                                               cGAGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATg                        >  1:575264/1‑76 (MQ=255)
                                                gAGCTCTATGCACACCAAGCAGCGGCTGACGTCTTCCGGCGTTGGCCGTTTTTTTATGTGTACGGAACTTCTATGa                       <  2:1614795/76‑1 (MQ=255)
                                                 aGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGaa                      >  2:3472154/1‑76 (MQ=255)
                                                 aGGTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGaa                      >  2:3276256/1‑76 (MQ=255)
                                                   gTCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATa                    <  1:4438189/76‑1 (MQ=255)
                                                    tCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTTAGGAACTTCTATGAATAc                   >  1:4152088/1‑76 (MQ=255)
                                                    tCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATAc                   >  1:1425720/1‑76 (MQ=255)
                                                    tCTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATAc                   >  2:2583135/1‑76 (MQ=255)
                                                     cTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACt                  >  2:4591353/1‑76 (MQ=255)
                                                     cTATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGCAAAGAACCTCTATTAATCCt                  >  2:4698289/1‑76 (MQ=255)
                                                      tATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACt                  >  2:4942085/1‑75 (MQ=255)
                                                      tATGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTg                 <  2:2671094/76‑1 (MQ=255)
                                                       aTGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGt                >  1:1395658/1‑76 (MQ=255)
                                                       aTGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGCTTTTTTATGTGTAAGGAACTTATATGAATACTGt                >  1:4789528/1‑76 (MQ=255)
                                                        tGCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTc               >  1:5714900/1‑76 (MQ=255)
                                                         gCATAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATAATGTCg              >  1:3717285/1‑76 (MQ=255)
                                                           aTAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGc            <  1:2640859/76‑1 (MQ=255)
                                                            tAAAAAGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCa           <  1:1105982/76‑1 (MQ=255)
                                                             aaaaaGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGccc          >  2:722439/1‑74 (MQ=255)
                                                              aaaaGCAGCGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACg         <  1:2725103/76‑1 (MQ=255)
                                                                aaGCAGCGGCGGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGcc       <  2:5727407/76‑1 (MQ=255)
                                                                 agcagcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGACTACTGTCGGGccgccg      >  2:594008/1‑76 (MQ=255)
                                                                 agcagcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCg      >  2:4439986/1‑76 (MQ=255)
                                                                  gcagcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGt     <  2:5848698/76‑1 (MQ=255)
                                                                   cagcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGtt    <  2:4934717/76‑1 (MQ=255)
                                                                    agcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTg   <  1:5209934/76‑1 (MQ=255)
                                                                     gcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTATGGAACTTCTATGAATACTGTCGGCACGCCGTTGc  <  1:417113/76‑1 (MQ=255)
                                                                     gcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGc  <  2:4065903/76‑1 (MQ=255)
                                                                     gcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGc  >  1:4834259/1‑76 (MQ=255)
                                                                     gcGGCTGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGc  >  1:3330450/1‑76 (MQ=255)
                                                                           |                                                                      
ACGCAGACGTATTCCTTTTTTGTTGTAAGTGAGACCTTGCCGGAAGGCGAGGTCTATGCATAAAAAGCAGCGGCTGGACGTCTTCCGACGTTGGCCGTTTTTTTATGTGTAAGGAACTTCTATGAATACTGTCGGCACGCCGTTGC  >  NC_000913/3236101‑3236246

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: