Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 36,174 A→C 19.2% intergenic (+59/+24) PP_0032 → / ← PP_0033 hypothetical protein/undecaprenyl‑glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294736,1740AC19.2% 111.2 / 22.1 57intergenic (+59/+24)PP_0032/PP_0033hypothetical protein/undecaprenyl‑glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (23/22);  new base C (6/5);  total (29/28)
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.94e-01

AGCGATGACTCTTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGAT  >  NC_002947/36026‑36314
                                                                                                                                                    |                                                                                                                                            
aGCGATGACTCTTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAATAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTaaaa                                                                                                                                              <  1:116392/149‑1 (MQ=255)
        ctctTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtt                                                                                                                                      >  7:63621/1‑149 (MQ=255)
         tctTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCttt                                                                                                                                     >  5:204954/1‑149 (MQ=255)
         tctTAGGACAAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAACGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGTCGGTTAAATAGGATTTCAAATGCCGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCttt                                                                                                                                     >  2:255730/1‑149 (MQ=255)
             cGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                  <  7:241291/148‑1 (MQ=255)
             cGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCtttttt                                                                                                                                  <  8:63621/148‑1 (MQ=255)
                  cAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCtt                                                                                                                            <  6:204954/149‑1 (MQ=255)
                        aaacaaCAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAAtt                                                                                                                      >  1:64804/1‑149 (MQ=255)
                          acaacaAGTCGCTCAAGAGCAAAAGAAAGAAAGAAAGGATGGTAAGTAACAATCCGACAATTAATCGGCTGTAAAATAGGACTTCAAAACACGAAAGAAATCGAAAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGAGCTTGGAAta                                                                                                                      >  1:47477/1‑146 (MQ=255)
                          aaaacaaGTATCTAAAGAGCAAAAGAAATAAAGCAAGGATTGAAAGGAACAATCCGAAAATTAATCGGCGGTTAAATAGGATTTCAAATGACGAAATAAATCGNCAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTa                                                                                                                    <  1:255730/147‑1 (MQ=255)
                              caaGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCt                                                                                                                <  4:188465/149‑1 (MQ=255)
                                 gTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGTTTTACCCTGGCtttttt                                                                                                                                  >  7:316298/1‑128 (MQ=255)
                                 gTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGTTTTACCCTGGCtttttt                                                                                                                                  <  8:316298/128‑1 (MQ=255)
                                             cAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGTTTTACCCTGGCTTTTTtgtc                                                                                                                               >  7:154592/1‑118 (MQ=38)
                                             cAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGTTTTACCCTGGCTTTTTtgtc                                                                                                                               <  8:154592/119‑2 (MQ=38)
                                                aagaaagaaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGtt                                                                                              >  4:19968/1‑149 (MQ=255)
                                                      gaaagCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCa                                                                                        <  7:257346/149‑1 (MQ=255)
                                                         agCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCg                                                                                                             >  3:256919/1‑125 (MQ=255)
                                                         agCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCg                                                                                                             <  4:256919/125‑1 (MQ=255)
                                                           cAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGa                                                                                                                         >  3:351860/1‑111 (MQ=255)
                                                           cAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGa                                                                                                                         <  4:351860/111‑1 (MQ=255)
                                                                      aGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGcc                                                                         <  1:106009/148‑1 (MQ=255)
                                                                      aGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGc                                                                          >  3:29957/1‑147 (MQ=255)
                                                                      aGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCa                                                                        <  2:64804/149‑1 (MQ=255)
                                                                      aGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCa                                                                        <  8:228239/149‑1 (MQ=255)
                                                                        gAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTa                                                                                                                    >  7:322028/1‑103 (MQ=255)
                                                                        gAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTa                                                                                                                    <  8:322028/103‑1 (MQ=255)
                                                                            aaGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCt                                                                  >  5:109747/1‑149 (MQ=255)
                                                                              gCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGttt                                                                                             <  5:89895/120‑1 (MQ=255)
                                                                              gCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGttt                                                                                             >  6:89895/1‑120 (MQ=255)
                                                                                cGACAATTAATCGGCGGTTAAATAGGACTTCAAATAACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCgg                                                               >  6:56894/1‑148 (MQ=255)
                                                                                         aTCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGc                                                     >  8:104505/1‑149 (MQ=255)
                                                                                            ggcggTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCggg                                                  >  4:234536/1‑149 (MQ=255)
                                                                                            ggcggTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCgg                                                   <  7:104505/148‑1 (MQ=255)
                                                                                                gTTAAATATGACTTCAAATGACTATCGAAATCGCCAATAAAGCCAGCGCTAAACCTGGCTTTTTTGTGCTTGGATTTTCGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTCTGCGGGGCGt                                              <  6:329041/149‑1 (MQ=255)
                                                                                                   aaaTAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGttt                                            <  3:19968/148‑1 (MQ=255)
                                                                                                      tacgACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCaa                                        <  3:121637/146‑1 (MQ=255)
                                                                                                        ggACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATg                                      <  4:29957/149‑1 (MQ=255)
                                                                                                            ttCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGca                                   >  2:223359/1‑148 (MQ=255)
                                                                                                                  tGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAgg                            >  8:65543/1‑149 (MQ=255)
                                                                                                                      gaaagaaaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCAGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgct                        >  6:11129/1‑149 (MQ=255)
                                                                                                                            aaTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgctggctg                   <  3:234536/148‑1 (MQ=255)
                                                                                                                             aTCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCtggctggc                 <  6:305064/149‑1 (MQ=255)
                                                                                                                              tCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCg                <  5:56894/149‑1 (MQ=255)
                                                                                                                                caCAAAAAAGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGtt              >  3:49990/2‑149 (MQ=255)
                                                                                                                                 aCAAAAAAGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGTGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGtt              <  5:309440/148‑1 (MQ=255)
                                                                                                                                 aCAAAAAAGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAAATTAGCCTTCTTTCTTGCGAGGCGATTCCGCATTGCTTTATCCGGCCAGACAATCGGCAACTTTTGCGGGGCGTTTCCAATGGCACAGAAGGAGATGGCTGGCGtt              >  6:309440/1‑148 (MQ=255)
                                                                                                                                   aaaaaaGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCTGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgctggctg                   <  7:239148/141‑1 (MQ=255)
                                                                                                                                   aaaaaaGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCTGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAgctggctg                   >  8:239148/1‑141 (MQ=255)
                                                                                                                                   aaaaaaGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGc           >  2:101136/1‑149 (MQ=255)
                                                                                                                                     aaaaGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGTTGGCTGGCGTTCGCtt         >  8:296254/1‑149 (MQ=255)
                                                                                                                                       aaGCCATGGTAAATCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCTGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGGGTTCGCTTCg       <  2:47477/149‑1 (MQ=255)
                                                                                                                                       aaGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCATTCGCtt         >  2:319433/1‑147 (MQ=255)
                                                                                                                                       aaGCCAGGTTTTACCCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCATTCGCtt         <  1:319433/147‑1 (MQ=255)
                                                                                                                                         gCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGt     <  1:223359/149‑1 (MQ=255)
                                                                                                                                           cAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGa   >  6:314028/1‑149 (MQ=255)
                                                                                                                                            aGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGAt  >  8:334934/1‑149 (MQ=255)
                                                                                                                                            aGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGAATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCTGGGCGTTTCCAATGGCACAGCTGGAGCTGGCTGGCGTTCGCTTCGGTGAt  >  7:258644/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
AGCGATGACTCTTCGGACCAGAACAAACAACAAGTCGCTCAAGAGCAAAAGAAAGAAAGCAAGGATGGTAAGGAACAAGCCGACAATTAATCGGCGGTTAAATAGGACTTCAAATGACGAAAGAAATCGACAAAAAAGCCAGGGTAAAACCTGGCTTTTTTGTGCTTGGAATTTAGCCTTCGTTCTTGCGAGGCGTTTCCGCATTGCTTTCTCCGGCCAGACCCTCGGCACCTTGTGCGGGGCGTTTCCAATGGCACAGCAGGAGCTGGCTGGCGTTCGCTTCGGTGAT  >  NC_002947/36026‑36314

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