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.G88.1% 136.6 / 11.4 43intergenic (‑131/+57)frmR/yaiOregulator protein that represses frmRAB operon/outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (4/1);  new base G (21/16);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCA  >  NC_000913/379865‑380154
                                                                                                                                                       |                                                                                                                                                    
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGaa‑‑a                                                                                                                                                       <  2:1198640/148‑1 (MQ=255)
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑ggg                                                                                                                                                 >  2:221028/1‑151 (MQ=255)
tcCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gg                                                                                                                                                  <  1:1059876/150‑1 (MQ=255)
 ccGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑ggg                                                                                                                                                 <  2:913777/150‑1 (MQ=255)
 ccGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑ggg                                                                                                                                                 <  2:374757/150‑1 (MQ=255)
  cGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑g                                                                                                                                                   <  2:347828/147‑1 (MQ=255)
           gggCGTTTCGCACCCCTTCACGTGGATGTAAGATTCTAGAGGGGGGTATTCTATATGTCAATGCATCCCCCCCTATAGTACAATACACGGGAGAATAAAATTTATCCGGTGAATGCGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gggggg‑‑ggggg                                                                                                                                       <  1:1050290/148‑2 (MQ=255)
            ggCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gggggg‑‑ggggggct                                                                                                                                    >  2:723633/1‑146 (MQ=255)
                tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCgg                                                                                                                              >  2:469283/1‑151 (MQ=255)
                tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gggggg‑‑ggggggctaatcg                                                                                                                               >  2:972538/1‑142 (MQ=255)
                tttCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑gggggg‑‑ggggggctaatcg                                                                                                                               >  1:577816/1‑142 (MQ=255)
                             cctcTTCATTCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGgcc                                                                                                                   >  2:877786/3‑151 (MQ=255)
                              atcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGgccgc                                                                                                                 >  1:856937/1‑151 (MQ=255)
                              atcTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAAAGTGGTCGGAAAACAAAGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAAAGGGAGGGGAAGgccccc                                                                                                                >  2:821124/1‑151 (MQ=255)
                                     atgaattaacctaaagggggTATTCTTATAGTTAAATCCTACCCCCCCTAAGTAAAAAACCCGGGGAAAAAAAATTTTTCCGGGGGAAGGGGGCGGAAAAAAAAAGGGAA‑‑AAGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                          <  1:336024/136‑1 (MQ=255)
                                        aaaatccTAAAGGGGGGTTTTTTAATTTTAAAACCATCCCCCCCTATAGTTAAAAACCGTGGAAAAAAAAATTTATCCGGGGAATGTGGTCGGAAAAAAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGAt                                                                                                          <  2:20392/142‑1 (MQ=255)
                                         atatACTATAGGGGGGTATTATAAATGTCAATGCATCCCCCCCTATAGTATAAAACACAAGAGAATAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGAAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgg                                                                                                     <  2:321147/151‑1 (MQ=255)
                                           ataCTATAGGGGGGTATTCTATATGTCAAGGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATGGGCAGGGAAGGACGCCCCGGATAgcgg                                                                                                     >  1:823374/1‑149 (MQ=255)
                                                ataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGggcggc                                                                                                >  2:870411/1‑150 (MQ=255)
                                                ataGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGAAGGGAAGGCCGCCCCGGATAGCGggcggc                                                                                                >  2:688760/1‑150 (MQ=255)
                                                   ggggtgTAGTCTTTATTTAAAGGCTTCCCCCCTTTTAATACAAAAACCAGGAGAAAAAAATTTATCCGGCAAAGGGGCTCGGAAAACAAAGGGAAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc                                                                                                   <  1:1033070/144‑1 (MQ=255)
                                                   ggggggTTTTCTATATTTCAATGCATACCCCCCTATAGTATAAAACACTGGAAAAAAAAATTTATCCGGTGAATGTGGGCGGAAAAAAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAg                                                                                              <  2:50050/149‑1 (MQ=255)
                                                   ggggggTTTTCAAATTGTCAATGAATACCCCCCTATTATATAAAAACGTGGAAAAAAAAATTTTTCCGGGGAATGGGGGCGGAAAAAAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGa                                                                                             <  2:646925/150‑1 (MQ=255)
                                                    gggggTATTCTATAGGGCAATACATCCCCCCCTATAATAAATAACACTGGAGAAAAAAATTTATCCGGGGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAgcgggc                                                                                                   <  1:233016/143‑1 (MQ=255)
                                                          aTTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑GGGGGG‑‑GGGGGCCTATTCGCGAGGGAAGGCCGCCCCGGATGGCGGGCGGCAGAGGGAAt                                                                                       >  2:1104775/1‑149 (MQ=255)
                                                            tCTATATGTCAATGCAGACCCCCCTATAGTAAATAAAACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAAGGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAg                                                                                    <  1:326440/150‑1 (MQ=255)
                                                            tCTATATGACAATGCAACCCCCCCTAAATTATAAAACAGGGAAAAAAAAAATTTAGCGGGGGAATGTGGGGGGAAAAAAAAGGGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATc                                                                                      <  2:725654/148‑1 (MQ=255)
                                                                tatGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAAtt                                                                                >  2:536718/1‑150 (MQ=255)
                                                                 atGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑GGGGGG‑‑GGGGGCCTAACCGGCGGGGAAGGCCGCCCCGGGTACGGGGCCGAAGAAGGAGTCAGAttt                                                                                >  1:1084852/1‑149 (MQ=255)
                                                                    gTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTcc                                                                             >  2:1158366/1‑150 (MQ=255)
                                                                           cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGAc                                                                    >  2:900501/1‑151 (MQ=255)
                                                                           cATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAg                                                                      >  2:359745/1‑149 (MQ=255)
                                                                            aaaCCCCCTTATATTATTAAACACTGGAGAAAAAAATTTATCCGGGGAATGGGGTGGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCa                                                                       <  1:869114/146‑1 (MQ=255)
                                                                                         gTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt                                                        >  2:365614/1‑150 (MQ=255)
                                                                                            tataACACGGGAGAATAAAATTTATCCGGTTAATGTGGTCGGAAAACAAGGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGtt                                                        <  1:972000/147‑1 (MQ=255)
                                                                                                           tAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGGGGGCGGCAGAAGGAATCAGAATTTCCAGGACAGACGGGCTGCAAGTTGCAGACGGTTAAAatcat                                      >  2:1087632/1‑150 (MQ=255)
                                                                                                            aaaaTTTACCCGGTAAATGGGGTCGAAAAACAAAGAAAAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAatcatc                                     <  2:1163006/150‑1 (MQ=255)
                                                                                                             aaaTTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCAGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCgg                                   >  1:968885/1‑151 (MQ=255)
                                                                                                              aaTTTATCCGGTGAATGTGGTCGGAAAACAACGTGGAA‑‑A‑‑‑GGGGGG‑‑GGAC‑‑ATAATCGGCCGGGAAGGCGTGCATGGATAGCGGCTGGCCGAGGGGCAGAAAATTTCGAGTTCTATCGGGGTGAAAGTTGCGGCGCGTTAAAATCATCGGTTg                                >  1:654595/1‑151 (MQ=255)
                                                                                                                  tATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTg                                >  2:780337/1‑149 (MQ=255)
                                                                                                                    tCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATAGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTgggg                             >  1:889049/1‑150 (MQ=255)
                                                                                                                     ccGGTGAATGGGGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgt                            <  1:59199/150‑1 (MQ=255)
                                                                                                                        gTGAGTGTGGTGGAAAAAAAAAGAGGAAGGG‑‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGgtgt                          <  2:141403/149‑1 (MQ=255)
                                                                                                                         tGAATGGGGTCGAAAAAAAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCg                        <  2:751247/150‑1 (MQ=255)
                                                                                                                              gtgGTCGGAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAcc                    >  2:443560/1‑149 (MQ=255)
                                                                                                                                    ggAAAACAAAGAGGCA‑‑AGGGGGGGG‑‑GGGG‑‑CAAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGCATTTCCAGGTAAGACGGGCTGCAAGTGTCAGACCGTTAAAATCTTCGGTTGGCGTGTCGTACCACACttt             >  1:650826/1‑150 (MQ=255)
                                                                                                                                     gAAAACAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTAc           >  1:97617/1‑151 (MQ=255)
                                                                                                                                      aaaaCAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTa            <  1:916267/149‑1 (MQ=255)
                                                                                                                                        aaCAAAGAGGAA‑‑AGGGGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGGGGGCGGCAGAAGGAATAAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt         >  2:796670/1‑150 (MQ=255)
                                                                                                                                           aaaGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCg     >  1:276459/1‑150 (MQ=255)
                                                                                                                                            aaGAGGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc   <  2:763111/151‑1 (MQ=255)
                                                                                                                                            aaGAGGAA‑‑A‑‑‑GGGGGGGCGGGG‑‑CTAAACGGCAGGGAAGGCCCCCCCGGATAGCGGGCGGGAGAAGGAACCAGAATTTCCAGGTCAGCCGGGCCGCAAGTTGCAGACCGTTAAAATCATCGGGTGGGGTGCCGCACCCCCCTTTACCTGcc      >  1:32678/1‑149 (MQ=255)
                                                                                                                                              gagGAA‑‑AG‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCa  <  1:469269/150‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
TCCGGAGTACTGGGCATTTCGCACCTCATCATCTGCATGCAATATACTATAGGGGGGTATTCTATATGTCAATGCATACCCCCCTATAGTATATAACACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAA‑‑A‑‑‑GGGGGG‑‑GGGG‑‑CTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCA  >  NC_000913/379865‑380154

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