Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 380,022 (G)10→11 intergenic (‑141/+47) frmR ← / ← yaiO regulator protein that represses frmRAB operon/outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,0121.G83.0% 172.2 / 25.4 59intergenic (‑131/+57)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (6/4);  new base G (29/20);  total (35/24)
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.00e+00

TCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTT  >  NC_000913/379883‑380143
                                                                                                                                  |                                                                                                                                        
taagagacagtcatcTGCATGCAATATACTATAGGGGGGTATTTTATATGTCAAGTAATAGCCCCCTATAGTATATAACATTGGAGAAAAAAAGTTCTCCGGTAAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑ggggggg‑‑gg                                                                                                                             <  1:1104202/129‑1 (MQ=255)
                     caATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCTCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGcctgtctctta                                                                                                        >  2:1104202/1‑129 (MQ=255)
                               taGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATa                                                                                              >  2:827708/1‑139 (MQ=255)
                                  gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                           >  2:765123/1‑139 (MQ=255)
                                  gggggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                           >  1:910274/1‑139 (MQ=255)
                                  aacgaacagttattTGTGAATGCATACCCCCCTATAGTATATAACTCTGGAGAATAAAATTTATCCGGTGAATGTGGGTGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                           <  2:1641488/129‑1 (MQ=255)
                                    gggTATTCTATATTTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                         <  1:827708/139‑1 (MQ=255)
                                    gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGG‑GGGGGGG‑‑GGGC‑AAATCGGCAGGGAAGGCCGCCCCGGATAgggg                                                                                          >  2:1490537/1‑139 (MQ=255)
                                    gggTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAgcggg                                                                                         >  2:896970/1‑139 (MQ=255)
                                       tATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc                                                                                     >  1:1675316/1‑139 (MQ=255)
                                                 gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑ggggggg‑‑ggggcta                                                                                                                        >  2:462426/1‑91 (MQ=255)
                                                 gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑ggggggg‑‑ggggcta                                                                                                                        <  1:462426/95‑5 (MQ=255)
                                                    aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCg                                                                                                                    >  2:869762/1‑96 (MQ=255)
                                                    aaTGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGGCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCg                                                                                                                    <  1:869762/96‑1 (MQ=255)
                                                      tGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACCAAGAGGAAA‑‑‑GGGGGGG‑‑GGGCTTAATGGGCAGGGAAGTCCGCCCCGGATAGCGGGCGGCCGAAgctgtctctt                                                                       >  1:1641488/1‑130 (MQ=255)
                                                            ccccccTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGCAAACGGGCGGCAGAAGGAATCAGAATTTCCa                                                                 >  1:1129504/1‑139 (MQ=255)
                                                                         tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCgg                                                                                                                   >  1:1194864/1‑76 (MQ=255)
                                                                         tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCgg                                                                                                                   >  1:1138782/1‑76 (MQ=255)
                                                                         tataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCgg                                                                                                                   >  2:1543772/1‑76 (MQ=255)
                                                                         tataACACTGGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCgg                                                                                                                   <  2:1138782/76‑1 (MQ=255)
                                                                         taaaaCACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCgg                                                                                                                   <  2:1194864/73‑1 (MQ=255)
                                                                         aataACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCgg                                                                                                                   <  1:1543772/75‑1 (MQ=255)
                                                                           taaacagggAGAATAAAATTTATCCGGTGAAGGGGGGGGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGc                                                  <  1:514333/132‑1 (MQ=255)
                                                                              cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                                  >  2:1370232/1‑104 (MQ=255)
                                                                              cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                                  <  1:1370232/104‑1 (MQ=255)
                                                                              cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGAGGGGCTGCAAg                                               >  2:1697680/1‑139 (MQ=255)
                                                                              cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                               >  2:340298/1‑139 (MQ=255)
                                                                              cacTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAg                                               >  1:1631244/1‑139 (MQ=255)
                                                                                 tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGGGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggg                                                      <  2:1766966/129‑1 (MQ=255)
                                                                                 tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGG‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt                                             >  2:475133/1‑139 (MQ=255)
                                                                                 tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACggg                                                      >  1:1766966/1‑129 (MQ=255)
                                                                                 tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                            >  2:514333/1‑139 (MQ=255)
                                                                                 tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                            >  1:1105092/1‑139 (MQ=255)
                                                                                 tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTg                                            >  1:509981/1‑139 (MQ=255)
                                                                                 tGGAGAATAAAATTTAACCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGG‑GGGGGGG‑‑GGGG‑TAATTGGTAGGGAAGGCCGCCCCGAAGAGCGGGCGGGAGAAGGAATAAGAATTTCCAGGTCAGACGGGCAGCAAGtt                                             >  2:231578/1‑139 (MQ=255)
                                                                                    agaATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGGGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAg                                         <  1:340298/139‑1 (MQ=255)
                                                                                        tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                        >  2:58723/1‑138 (MQ=255)
                                                                                        tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGa                                        <  1:58723/138‑1 (MQ=255)
                                                                                        tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCg                                     >  2:1044778/1‑139 (MQ=255)
                                                                                          aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGtt                                   <  2:910274/139‑1 (MQ=255)
                                                                                                  ccGGTGAATGTGGTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcat                           <  2:1129504/139‑1 (MQ=255)
                                                                                                   cGGTGAATGTGGTCGGAAAACAAAGAGGAAAGG‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcat                           <  1:475133/139‑1 (MQ=255)
                                                                                                      tGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGG‑‑GGGC‑AAATAGGCAGGGAAGGCCGCCCCGGAGAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGCCGGGCGGCAAGTTGCAGACCTTTAAAATCATCg                         >  2:1611283/1‑139 (MQ=255)
                                                                                                      tGAATGTGGTCGGAAAACAAAGAGGAAAGG‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAgg                                                                                                          >  2:253436/1‑57 (MQ=255)
                                                                                                      tGAATGTGGTCGGAAAACAAAGAGGAAAGG‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAgg                                                                                                          <  1:253436/57‑1 (MQ=255)
                                                                                                      tGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGtt                      >  1:138381/1‑139 (MQ=255)
                                                                                                      tGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGtt                      >  2:613933/1‑139 (MQ=255)
                                                                                                        aaTGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgg                    <  1:613933/139‑1 (MQ=255)
                                                                                                             gggcGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtg                <  2:1391679/136‑1 (MQ=255)
                                                                                                              gTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgg                    >  2:174865/1‑134 (MQ=255)
                                                                                                              gTCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgg                    <  1:174865/134‑1 (MQ=255)
                                                                                                              gTCGGAAAACAAAGAGGAAA‑‑‑‑GGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                  <  2:954345/118‑1 (MQ=255)
                                                                                                              gTCGGAAAACAAAGAGGAAA‑‑‑‑GGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                  >  1:954345/1‑118 (MQ=255)
                                                                                                               tCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCgg                                                                                                                   >  2:1782612/1‑38 (MQ=255)
                                                                                                               tCGGAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCgg                                                                                                                   <  1:1782612/38‑1 (MQ=255)
                                                                                                                  gAAAACAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATCTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa           >  1:1415373/1‑139 (MQ=255)
                                                                                                                  gAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                           <  2:960179/58‑1 (MQ=255)
                                                                                                                  gAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAgcg                                                                                           >  1:960179/1‑58 (MQ=255)
                                                                                                                  gAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc          <  2:1675316/139‑1 (MQ=255)
                                                                                                                   aaaaCAATGGGGAAA‑‑‑GGGGGGGGCGGGC‑TAAACGGCAGGGAAGGCCGCCCCGGATAGCGGCAGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGATGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa           <  1:231578/139‑1 (MQ=255)
                                                                                                                   aaaaCAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc          <  1:896970/139‑1 (MQ=255)
                                                                                                                       cAAAGAGGAAAG‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACcaca      <  2:509981/139‑1 (MQ=255)
                                                                                                                           gagGAAAG‑‑GGGGGGG‑‑GGGC‑TAACCGGCAGGGAAGGCCGCCCCGCAAACCGGGCGGCAGAAGGAATCAGATTTTCCGGGTCGGCCGGGCTGCAAGTTGCGGACCGTTAAAATCATCGGTTGGGGGGTCGTACCACACttt  >  2:1200863/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                        
TCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAA‑‑‑GGGGGGG‑‑GGGC‑TAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTT  >  NC_000913/379883‑380143

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