Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,512,324 A→T 6.7% Y331F (TAC→TTC)  PP_3108 → rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,512,3240AT 6.7% 173.6 / 2.5 60Y331F (TAC→TTC) PP_3108rhs family protein
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base T (2/2);  total (28/32)
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.18e-01

GAGCCTGCGCCGCCAGGGGGCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCGACGGCACACTGACCGACATCA  >  NC_002947/3512184‑3512470
                                                                                                                                            |                                                                                                                                                  
gagCCTGCGCCGCCAGGGGGCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGc                                                                                                                                            <  5:12218/149‑1 (MQ=255)
   ccTGCGCCGCCAGGGGGCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGa                                                                                                                                         <  8:256642/149‑1 (MQ=255)
   ccTGCGCCGCCAGGGGGCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGa                                                                                                                                         >  5:119403/1‑149 (MQ=255)
     tGCGCCGCCAGGGGGCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGAcc                                                                                                                                       >  4:45814/1‑149 (MQ=255)
              aGGGGGCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAAc                                                                                                                              <  1:112121/149‑1 (MQ=255)
               gggggCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAAcc                                                                                                                             <  4:376483/149‑1 (MQ=255)
                 gggCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTg                                                                                                                           <  6:119403/149‑1 (MQ=255)
                     aTTCATCTAACTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAAACGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGttct                                                                                                                       >  5:151067/1‑149 (MQ=255)
                      ttCATCTACCTCCCCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGttctt                                                                                                                      <  4:264206/149‑1 (MQ=255)
                        cATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCCACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTg                                                                                                                                             <  8:289618/124‑1 (MQ=255)
                        cATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCCACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTg                                                                                                                                             >  7:289618/1‑124 (MQ=255)
                             aCCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACttt                                                                                                               >  2:154130/1‑149 (MQ=255)
                              ccTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTg                                                                                                              <  6:21640/149‑1 (MQ=255)
                                tCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTTCCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGCCAACCTGTTCTTCTACTTTGGc                                                                                                            >  3:199678/1‑149 (MQ=255)
                                 cACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTgg                                                                                                             <  4:166782/147‑1 (MQ=255)
                                 cACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTgg                                                                                                             >  3:166782/1‑147 (MQ=255)
                                    cGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTg                                                                                                              <  4:50492/143‑1 (MQ=255)
                                    cGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTg                                                                                                              >  3:50492/1‑143 (MQ=255)
                                       cAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTTGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGc                                                                                                     <  3:45814/149‑1 (MQ=255)
                                         aCCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCa                                                                                                   <  5:146584/149‑1 (MQ=255)
                                                  gcgAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCaa                                                                                           >  8:297924/1‑148 (MQ=255)
                                                        tCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGaa                                                                                    <  7:297924/149‑1 (MQ=255)
                                                          ccGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGt                                                                                  <  1:341129/149‑1 (MQ=255)
                                                                gTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGc                                                                             <  1:154130/148‑1 (MQ=255)
                                                                  gACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTg                                                                           >  3:161587/1‑148 (MQ=255)
                                                                   aCCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTg                                                                           >  2:196545/1‑147 (MQ=255)
                                                                      cTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTTCATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCg                                                                      >  4:121968/1‑149 (MQ=255)
                                                                      cTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTTCATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCg                                                                      >  4:116683/1‑149 (MQ=255)
                                                                       tGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCg                                                                      >  1:330229/1‑148 (MQ=255)
                                                                       accAACCGCCGCAGCGGAAATACAACCCGAACAAGAAGGTGTACCTGCTGTGCAACGAGGGCAGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTACTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCg                                                                      >  4:373016/3‑148 (MQ=255)
                                                                         cAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTAt                                                                   >  2:208535/1‑149 (MQ=255)
                                                                          aaCCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGaa                                                                                    <  6:125675/131‑1 (MQ=255)
                                                                          aaCCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGaa                                                                                    >  5:125675/1‑131 (MQ=255)
                                                                          aaCCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAg                                                                        <  3:369392/143‑1 (MQ=255)
                                                                          aaCCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAg                                                                        >  4:369392/1‑143 (MQ=255)
                                                                          aaCCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATc                                                                  <  1:208535/149‑1 (MQ=255)
                                                                             cGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCg                                                                 >  4:5860/1‑147 (MQ=255)
                                                                             cGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCg                                                                 <  3:5860/147‑1 (MQ=255)
                                                                                       cATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCg                                                     <  1:196545/149‑1 (MQ=255)
                                                                                             cAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACt                                               >  2:325680/1‑149 (MQ=255)
                                                                                                  cGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTg                                          <  4:161587/149‑1 (MQ=255)
                                                                                                    cACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTcc                                            <  7:282560/145‑1 (MQ=255)
                                                                                                    cACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTcc                                            >  8:282560/1‑145 (MQ=255)
                                                                                                    cACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCa                                        <  1:325680/149‑1 (MQ=255)
                                                                                                         gCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCa                                   >  1:318363/1‑149 (MQ=255)
                                                                                                              tgtACCTGGTGTGCACCGAGGGCGGTCACTTCATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGc                              <  3:116683/149‑1 (MQ=255)
                                                                                                              tgtACCTGGTGTGCACCGAGGGCGGTCACTTCATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGc                              <  3:121968/149‑1 (MQ=255)
                                                                                                              tgtACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGc                              <  2:318363/149‑1 (MQ=255)
                                                                                                                 aCCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACg                                                           <  3:197954/117‑1 (MQ=255)
                                                                                                                 aCCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACg                                                           >  4:197954/1‑117 (MQ=255)
                                                                                                                                aGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCGACGGCACACTg            >  8:82742/1‑149 (MQ=255)
                                                                                                                                  ggCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCGACGTCACACTGAc          <  6:151067/149‑1 (MQ=255)
                                                                                                                                     ggTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCGACGGCACACTGACCGa       >  6:20111/1‑149 (MQ=255)
                                                                                                                                       tCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCGACGGca                 >  7:38039/1‑137 (MQ=255)
                                                                                                                                       tCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCGACGGca                 <  8:38039/137‑1 (MQ=255)
                                                                                                                                         aCTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTTGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCTACGGCACACTGACCGACAt    <  2:330229/148‑1 (MQ=255)
                                                                                                                                          cTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCGACGGCACACTGACCGACATCa  <  4:199678/149‑1 (MQ=255)
                                                                                                                                          cTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGGATCGAGAACGTCCTCGGCCACTTCCTGCACTTCCCCCGCACGCCCGACGGCACACTGCACGACAt    >  8:13416/1‑147 (MQ=255)
                                                                                                                                           tACATCCTGCAGACCCTCTACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCGGCGTATCGAGAACGTCCTCGGCCACTTCCTg                                          <  6:144625/108‑1 (MQ=255)
                                                                                                                                           tACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCCGCGTATCGAGAACGTCCTCGGCCACTTCCTg                                          >  5:144625/1‑108 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
GAGCCTGCGCCGCCAGGGGGCATTCATCTACCTCACCGACAACCAGGGCCGCGAAGTCCCGTTCGTGACCCTGCAACCGGGGCAGCGCATCTACAACCCGCACGAGCAGGTGTACCTGGTGTGCACCGAGGGCGGTCACTACATCCTGCAGACCCTCGACAACCTGTTCTTCTACTTTGGCGAAGTGCCAGACACCAACACCGAAGTGCCGCTGCAGCGTATCGAGAACGTCCTCGGCCACTTCCTGCACTTCACCCGCACGCCCGACGGCACACTGACCGACATCA  >  NC_002947/3512184‑3512470

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: ‑