Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,787,644 A→T 6.9% L50Q (CTG→CAG)  PP_5070 ← ParA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,787,6440AT 6.9% 175.2 / 2.9 58L50Q (CTG→CAG) PP_5070ParA family protein
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base T (2/2);  total (32/26)
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.97e-01

GTCGGCCAGCTCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCACGCG  >  NC_002947/5787496‑5787786
                                                                                                                                                     |                                                                                                                                               
gTCGGCCAGCTCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCa                                                                                                                                                 <  6:189930/149‑1 (MQ=255)
 tCGGCCAGCTCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAg                                                                                                                                                <  4:275164/149‑1 (MQ=255)
      cAGCTCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCgg                                                                                                                                           >  8:252647/1‑149 (MQ=255)
          tCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAg                                                                                                                                      >  4:86530/1‑149 (MQ=255)
                 tGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTgg                                                                                                                               <  3:179931/149‑1 (MQ=255)
                   gCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTggg                                                                                                                              <  3:86530/148‑1 (MQ=255)
                          ccACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCACCGGCCTGACCGGCGAggac                                                                                                                                   <  8:118719/136‑3 (MQ=255)
                          ccACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCACCGGCCTGACCGGCGAggac                                                                                                                                   >  7:118719/1‑134 (MQ=255)
                               tGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTg                                                                                                                 <  7:67877/149‑1 (MQ=255)
                                    gTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCAt                                                                                                            >  5:31207/1‑149 (MQ=255)
                                       ctCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATcc                                                                                                          <  7:252647/147‑1 (MQ=255)
                                        tCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGa                                                                                                                                       <  8:159850/118‑1 (MQ=255)
                                        tCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGa                                                                                                                                       >  7:159850/1‑118 (MQ=255)
                                          gAACGGGGTTTCGTAGATATGTACCTTGTTCTTCTTGCTGAACGGGCCGCTGGAAAGGCTTTGGTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGGCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGt                                                                                                      <  5:312392/149‑1 (MQ=255)
                                           aaCGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTc                                                                                                     >  7:321918/1‑149 (MQ=255)
                                                      gTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAgg                                                                                          >  1:336259/1‑149 (MQ=255)
                                                        aGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTc                                                                                        >  4:364294/1‑149 (MQ=255)
                                                        aGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTc                                                                                        >  3:109009/1‑149 (MQ=255)
                                                            atCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCgg                                                                                     >  3:222817/1‑148 (MQ=255)
                                                               gACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGc                                                                                 >  5:4327/1‑149 (MQ=255)
                                                                acattGTTCTTATTGCTGAACGGCCCTCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATTCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCTCATGGGCATCCAGGTCGATAAACAAGTTCCGGTAGcc                                                                                >  2:164986/4‑149 (MQ=255)
                                                                     gTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCTATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCAt                                                                           >  1:279033/1‑149 (MQ=255)
                                                                         ttcttGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCACCGGCCTGACCGGCGAg                                                                                                                                      <  2:367826/86‑1 (MQ=255)
                                                                         ttcttGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCACCGGCCTGACCGGCGAg                                                                                                                                      >  1:367826/1‑86 (MQ=255)
                                                                                tGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAcc                                                                >  8:37533/1‑149 (MQ=255)
                                                                                tGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAcc                                                                >  3:192481/1‑149 (MQ=255)
                                                                                   aCGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAccgcc                                                             <  6:4327/149‑1 (MQ=255)
                                                                                    cGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAccgc                                                              <  2:205031/147‑1 (MQ=255)
                                                                                    cGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGAccgc                                                              >  1:205031/1‑147 (MQ=255)
                                                                                       cccGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGTGGATGTCCTCGCCGGTCAGGACGGTGAGGTATTGGGTCGAGTTCGCCTGGTCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACcgccgc                                                           >  2:86080/1‑147 (MQ=255)
                                                                                       cccGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACcgccgc                                                           >  1:50606/1‑147 (MQ=255)
                                                                                         cGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATACAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAgg                                                       >  5:221783/1‑149 (MQ=255)
                                                                                              gAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCa                                                  <  4:109009/149‑1 (MQ=255)
                                                                                                aaGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGGCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTc                                                                                                     >  1:337824/1‑96 (MQ=255)
                                                                                                aaGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGCCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTc                                                                                                     <  2:337824/96‑1 (MQ=255)
                                                                                                 aGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGAACTGACCGCCGCCAGGTTGCAGGc                                               <  2:279033/149‑1 (MQ=255)
                                                                                                  ggCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCa                                              <  1:370114/149‑1 (MQ=255)
                                                                                                   gCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCaa                                             >  7:349456/1‑149 (MQ=255)
                                                                                                   gCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCaa                                             <  4:222817/149‑1 (MQ=255)
                                                                                                    cTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCAGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAAt                                            >  8:285162/1‑149 (MQ=255)
                                                                                                     tttGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATg                                           >  6:148971/1‑149 (MQ=255)
                                                                                                        gCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGTGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCt                                         <  7:285162/148‑1 (MQ=255)
                                                                                                        gCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCt                                         <  3:364294/148‑1 (MQ=255)
                                                                                                          ttGAAGAAATCGGCGATGCCCATGGGGATGGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTc                                                                                        >  1:216044/1‑100 (MQ=255)
                                                                                                          ttGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTc                                                                                        <  2:216044/99‑1 (MQ=255)
                                                                                                              agaaATCGGCGATGCCCATGGGGATGTCCTCGCTGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGt                                                                                                      <  5:346210/81‑1 (MQ=255)
                                                                                                              agaaATCGGCGATGCCCATGGGGATGTCCTCGCTGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGt                                                                                                      >  6:346210/1‑81 (MQ=255)
                                                                                                               gaaATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGc                                 <  7:37533/149‑1 (MQ=255)
                                                                                                                        gATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCAc                        <  8:349456/149‑1 (MQ=255)
                                                                                                                            cccATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGGCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCtt                     <  1:164986/148‑1 (MQ=255)
                                                                                                                               aTGGGGATGTCCTCGCCGGTCAGGCCGGTTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTc                   >  1:27367/1‑148 (MQ=255)
                                                                                                                               aTGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTc                   <  2:27367/147‑1 (MQ=255)
                                                                                                                                  gggATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGAtt              <  6:31207/149‑1 (MQ=255)
                                                                                                                                    gATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCAGCTTTCTGATTg             >  4:328240/1‑148 (MQ=255)
                                                                                                                                         tCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAacac        >  8:196785/1‑149 (MQ=255)
                                                                                                                                          ccTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTc                   >  2:121622/1‑137 (MQ=255)
                                                                                                                                          ccTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTc                   <  1:121622/137‑1 (MQ=255)
                                                                                                                                           cTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAAcacca      >  3:98672/1‑149 (MQ=255)
                                                                                                                                               ccGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCAcgcg  >  3:23154/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GTCGGCCAGCTCGGCAGTGGCGGTGACCACATGCAGGTTGTCGAACGGGGTTTCGTAGATATCGACCTTGTTCTTCTTGCTGAACGGCCCGCTGGAAAGGCTTTGCTTGAAGAAATCGGCGATGCCCATGGGGATGTCCTCGCCGGTCAGGCCGGTGAGGTATTGGGTCGAGTTCGCCTGGGCATCCAGGTCGATCAACAAGGTCCGGTAGCCTTCATTGGCACTGACCGCCGCCAGGTTGCAGGCAATGCTCGACTTGCCCACGCCACCTTTCTGATTGAACACCACGCG  >  NC_002947/5787496‑5787786

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