Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,087,769 C→G 7.1% R259P (CGC→CCC)  metR‑II ← homocysteine‑binding DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,087,7690CG 7.1% 176.7 / 3.5 56R259P (CGC→CCC) metR‑IIhomocysteine‑binding DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base G (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

CGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGG  >  NC_002947/3087624‑3087914
                                                                                                                                                  |                                                                                                                                                 
cGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCt                                                                                                                                                >  7:4656/1‑149 (MQ=255)
  ggCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTg                                                                                                                                              >  2:343205/1‑149 (MQ=255)
     cTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCg                                                                                                                                           <  1:211434/149‑1 (MQ=255)
        aCGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGtt                                                                                                                                        <  4:229589/149‑1 (MQ=255)
           ccACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACa                                                                                                                                     <  3:158745/149‑1 (MQ=255)
                   acacCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcgcg                                                                                                                             >  3:70929/1‑149 (MQ=255)
                    cacCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcgcgc                                                                                                                            <  7:6431/149‑1 (MQ=255)
                     acCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCcgcgc                                                                                                                            >  4:357188/1‑148 (MQ=255)
                       cTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCt                                                                                                                                                <  8:7603/126‑1 (MQ=255)
                       cTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCt                                                                                                                                                >  7:7603/1‑126 (MQ=255)
                         gCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGTCAGGCCGTGCTCGCCAAGCGCCTTGGCga                                                                                                                                          <  1:210948/131‑2 (MQ=255)
                         gCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCCTTGGCga                                                                                                                                          >  2:210948/1‑129 (MQ=255)
                          ccAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGtt                                                                                                                                        <  8:66532/131‑1 (MQ=255)
                          ccAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGtt                                                                                                                                        >  7:66532/1‑131 (MQ=255)
                          ccAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACt                                                                                                                       <  2:122639/148‑1 (MQ=255)
                            aGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTa                                                                                                                    <  3:121949/149‑1 (MQ=255)
                                  gtgAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGt                                                                                                               <  3:98207/148‑1 (MQ=255)
                                  gtgAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGt                                                                                                               <  8:4656/148‑1 (MQ=255)
                                     aaCGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATa                                                                                                                                   <  5:127064/125‑1 (MQ=255)
                                     aaCGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATa                                                                                                                                   >  6:127064/1‑125 (MQ=255)
                                     aaCGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCa                                                                                                            <  6:322390/148‑1 (MQ=255)
                                          gaTGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGc                                                                                                             >  1:64519/1‑142 (MQ=255)
                                          gaTGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGc                                                                                                             <  2:64519/142‑1 (MQ=255)
                                              tCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGtt                                                                                                  <  3:357188/149‑1 (MQ=255)
                                                cTTGGCCGTGAGCAGAAAATCACGCATGTAGGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTgg                                                                                                >  1:32470/1‑149 (MQ=255)
                                                               aaaaTCACGCATGTACGGCACCTCCAGCATGACACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCAAAGGCGCTTGGCGGTTACATAGCCGCGCGAAATGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGAcgc                                                                                 >  5:126739/1‑149 (MQ=255)
                                                                aaaTCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGAcgc                                                                                 <  6:146088/148‑1 (MQ=255)
                                                                 aaTCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGAcgccg                                                                               >  4:287885/1‑149 (MQ=255)
                                                                        cATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCCTTGGCgagca                                                                                                                                       >  1:251976/1‑82 (MQ=255)
                                                                        cATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCCTTGGCgagca                                                                                                                                       <  2:251976/86‑5 (MQ=255)
                                                                              cGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCaa                                                                  >  1:149604/1‑149 (MQ=255)
                                                                               ggCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAAc                                                                 <  4:70929/149‑1 (MQ=255)
                                                                                          cATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatca                                                      >  1:24706/1‑149 (MQ=255)
                                                                                                aCTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCAtcatcatca                                                >  1:192840/1‑149 (MQ=255)
                                                                                                 cTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCAc                                               >  1:62888/1‑149 (MQ=255)
                                                                                                  tGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACg                                              >  1:137833/1‑149 (MQ=255)
                                                                                                   gcgcAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatcatc                                                    <  7:306310/142‑1 (MQ=255)
                                                                                                   gcgcAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatcatc                                                    >  8:306310/1‑142 (MQ=255)
                                                                                                       aGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAg                                         >  3:232184/1‑149 (MQ=255)
                                                                                                         tGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGtt                                       <  2:137833/149‑1 (MQ=255)
                                                                                                         tGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGtt                                       <  4:232184/149‑1 (MQ=255)
                                                                                                           cagcaTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTc                                      <  2:32470/148‑1 (MQ=255)
                                                                                                              caTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGaa                                  <  1:343205/149‑1 (MQ=255)
                                                                                                               aTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGaa                                  <  2:192840/148‑1 (MQ=255)
                                                                                                                  gagGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgtac                              <  2:24706/149‑1 (MQ=255)
                                                                                                                     gTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgtacgta                           >  4:216644/1‑149 (MQ=255)
                                                                                                                            tGCAGGACGTGCTCGCAAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTACAGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGACAGTTCTGAAGTACGTACCTGTGc                    >  2:342205/1‑149 (MQ=255)
                                                                                                                                   ccGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTc              <  2:149604/149‑1 (MQ=255)
                                                                                                                                      tGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGAc           <  3:216644/149‑1 (MQ=255)
                                                                                                                                       gCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACg          >  4:17308/1‑149 (MQ=255)
                                                                                                                                         tCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTc        >  3:332489/1‑149 (MQ=255)
                                                                                                                                           gCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCt                        <  4:175256/131‑1 (MQ=255)
                                                                                                                                           gCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCt                        >  3:175256/1‑131 (MQ=255)
                                                                                                                                           gCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcg       <  2:62888/148‑1 (MQ=255)
                                                                                                                                             cAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACCGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcggcg    >  7:115686/1‑149 (MQ=255)
                                                                                                                                              aaGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcggcgg   >  5:89281/1‑149 (MQ=255)
                                                                                                                                                 gcgcTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCggg  >  4:59094/1‑147 (MQ=255)
                                                                                                                                                 gcgcTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCggg  <  3:59094/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGG  >  NC_002947/3087624‑3087914

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: