Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,976,789 C→G 7.5% R151R (CGC→CGG PP_1770 → 3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,976,7890CG 7.5% 168.5 / 3.7 53R151R (CGC→CGGPP_17703‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (26/23);  new base G (2/2);  total (28/25)
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.95e-01

CGTCGGCAGCGCCAAGGGCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGCCACCAGCCATCTGC  >  NC_002947/1976651‑1976937
                                                                                                                                           |                                                                                                                                                    
cGTCGGCAGCGCCAAGGGCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAgg                                                                                                                                            <  5:262947/149‑1 (MQ=255)
cGTCGGCAGCGCCAAGGGCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAgg                                                                                                                                            <  4:288847/149‑1 (MQ=255)
cGTCGGCAGCGCCAAGGGCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAgg                                                                                                                                            <  4:126252/149‑1 (MQ=255)
 gtcgtcAGCGCCAAGGGCAACGTCCTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCTGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTGCCGGCGGAACAAGGt                                                                                                                                           <  2:177150/149‑1 (MQ=255)
    ggCAGCGCCAAGGGCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCAt                                                                                                                                        >  7:33305/1‑149 (MQ=255)
            cAAGGGCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATc                                                                                                                                       <  6:134054/142‑1 (MQ=255)
            cAAGGGCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATc                                                                                                                                       >  5:134054/1‑142 (MQ=255)
               gggCAACGTCGTGCGTGAGGCGCGTGCCTTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGCGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCATCTTCCGCCGGCACAACGTCATCCTCAccccg                                                                                                                              >  6:42189/1‑148 (MQ=255)
                ggCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGCGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCAcgccgc                                                                                                                             <  1:286186/149‑1 (MQ=255)
                       tcgtGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTc                                                                                                                                          >  6:187696/1‑128 (MQ=255)
                                  gcgcGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCccgc                                                                                                          >  8:34375/1‑149 (MQ=255)
                                          cGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTgg                                                                                                  >  8:219812/1‑149 (MQ=255)
                                            tCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACACCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAg                                                                                                                                             >  4:270911/1‑104 (MQ=255)
                                            tCTTCGGTGCCCGCCTGCAGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAg                                                                                                                                             <  3:270911/104‑1 (MQ=255)
                                                cGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTcc                                                                                               <  1:265610/146‑1 (MQ=255)
                                                 ggTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTgg                                                                                            <  3:345025/148‑1 (MQ=255)
                                                 ggTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTgg                                                                                            <  3:64045/148‑1 (MQ=255)
                                                     cccGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGATGAc                                                                                       <  3:130596/149‑1 (MQ=255)
                                                          cTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCt                                                                                  >  2:342156/1‑149 (MQ=255)
                                                             ccGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCtgtg                                                                               >  6:24631/1‑149 (MQ=255)
                                                             ccGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCtgtg                                                                               <  7:219812/149‑1 (MQ=255)
                                                               gcgcTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTgg                                                                                            <  6:210728/134‑1 (MQ=255)
                                                               gcgcTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGc                                                                             <  7:34375/149‑1 (MQ=255)
                                                               gcgcTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCCCTGGTCCTgg                                                                                            >  5:210728/1‑134 (MQ=255)
                                                                   ttCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTgc                                                                         >  5:75870/1‑149 (MQ=255)
                                                                      gTCCCGGGCCACCCCATCGCCGGTTCCGTGCAGAGGGGTGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCTTCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTgcgc                                                                       <  5:42189/148‑1 (MQ=255)
                                                                      gTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCt                                                                      <  5:165143/149‑1 (MQ=255)
                                                                       tCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTg                                                                     >  5:151902/1‑149 (MQ=255)
                                                                       tCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTg                                                                     >  3:148117/1‑149 (MQ=255)
                                                                                  cccATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGAtgt                                                          <  6:98017/149‑1 (MQ=255)
                                                                                  cccATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGAtgt                                                          >  5:98017/1‑149 (MQ=255)
                                                                                  cccATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGAtgt                                                          >  6:263748/1‑149 (MQ=255)
                                                                                   ccATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGAtgtg                                                         <  1:342156/149‑1 (MQ=255)
                                                                                      tCGCCGGGTCCGCGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAg                                                      <  5:370618/149‑1 (MQ=255)
                                                                                        gCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGc                                                     >  2:13644/1‑148 (MQ=255)
                                                                                                         cGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGcc                                   <  6:75870/149‑1 (MQ=255)
                                                                                                             ggTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAg                                                      <  7:313969/127‑1 (MQ=255)
                                                                                                             ggTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAg                                                      >  8:313969/1‑127 (MQ=255)
                                                                                                              gTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTgcg                                                                        <  4:85584/108‑1 (MQ=255)
                                                                                                              gTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGcac                                                   >  3:173373/1‑129 (MQ=255)
                                                                                                              gTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGcac                                                   <  4:173373/129‑1 (MQ=255)
                                                                                                                 gagGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACg                             >  4:68843/1‑148 (MQ=255)
                                                                                                                  agGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACg                             >  6:146021/1‑147 (MQ=255)
                                                                                                                          aaTGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGgc                   >  5:374076/1‑149 (MQ=255)
                                                                                                                          aaTGCCACCTTGTTCCGCCGCCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGGCCTGGTTGACCGCCACTGGCGTCCCCTGGACGCCGAAGTCGACCACAAGTACGAAGAGCGCCATGCCGAAGTACTggg                   >  8:234612/1‑148 (MQ=255)
                                                                                                                               gaCCTTGTGCCGGCGGAACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGc                                                                             <  2:309565/85‑1 (MQ=39)
                                                                                                                               gaCCTTGTGCCGGCGGAACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGc                                                                             >  1:309565/2‑86 (MQ=39)
                                                                                                                                aCCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGccac             >  6:351501/1‑149 (MQ=255)
                                                                                                                                 ccTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGccac             >  5:373841/1‑148 (MQ=255)
                                                                                                                                  cTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGccacca           <  5:24631/149‑1 (MQ=255)
                                                                                                                                  cTTGTGCCGGCGGAACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGccac             >  6:207467/1‑147 (MQ=255)
                                                                                                                                     gTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGCCACCAGc         >  7:131295/1‑148 (MQ=255)
                                                                                                                                      ttCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGCCACCAGCCa       >  8:26578/1‑149 (MQ=255)
                                                                                                                                        ccggcgGAACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGt                                                                                           >  8:138320/1‑63 (MQ=255)
                                                                                                                                        ccggcgGAACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGt                                                                                           <  7:138320/63‑1 (MQ=255)
                                                                                                                                          gccgGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATgtcggtc                                          <  3:150868/110‑1 (MQ=255)
                                                                                                                                          gccgGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATgtcggtc                                          >  4:150868/1‑110 (MQ=255)
                                                                                                                                          gccgGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGCCACCAGCCATCTg   <  8:131295/149‑1 (MQ=255)
                                                                                                                                           ccgGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGCCACCAGCCATCTGc  <  7:26578/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CGTCGGCAGCGCCAAGGGCAACGTCGTGCGTGAGGCGCGTGCCGTCTTCGGTGCCCGCCTGCCGCGCTTCGTCCCGGGCCACCCCATCGCCGGTTCCGAGCAGAGCGGGGTAGAGGCGTCCAATGCCACCTTGTTCCGCCGGCACAAGGTCATCCTCACGCCGCTGGCCGAAACCGACCCCGCCGCGCTGGTCCTGGTTGACCGCCTGTGGCGTGCGCTGGACGCCGATGTGGAGCACATGTCGGTCGAGCGCCATGACGAAGTACTGGCCGCCACCAGCCATCTGC  >  NC_002947/1976651‑1976937

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: