Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 3,167,544 C→T S227S (AGC→AGT PP_2781 → beta‑ketoacyl synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,167,5440CT94.6% 165.9 / ‑0.4 56S227S (AGC→AGTPP_2781beta‑ketoacyl synthase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (26/27);  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 = 9.68e-01

GTGTCGCCGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGCGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGCCGGTCCTCAGGCGTGCCTCGCCGA  >  NC_002947/3167402‑3167675
                                                                                                                                              |                                                                                                                                   
gtgtCGCCGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAAc                                                                                                                               <  5:280254/149‑1 (MQ=255)
  gtCGCCGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGCGCCAACgg                                                                                                                             >  8:172292/1‑149 (MQ=255)
   tCGCCGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGc                                                                                                                            >  4:269138/1‑149 (MQ=255)
     gCCGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTc                                                                                                                          >  7:52269/1‑149 (MQ=255)
     gCCGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTc                                                                                                                          >  7:140213/1‑149 (MQ=255)
      ccGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCg                                                                                                                         <  3:61040/149‑1 (MQ=255)
       cGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGa                                                                                                                        <  4:72264/149‑1 (MQ=255)
         gCAACCTGCCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGTCTTGCCGAATTGACGCGGCAATCGTTCTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGTTGTAATCAGTGCCAACGGCTCGAac                                                                                                                      <  6:29799/149‑1 (MQ=255)
              cTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACgct                                                                                                                 >  2:100737/1‑149 (MQ=255)
               tGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACgctg                                                                                                                >  2:308490/1‑149 (MQ=255)
                ggCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACgctgc                                                                                                               >  3:221726/1‑149 (MQ=255)
                gccAGCCTTGAGTGAGGCTACGTCCGATATGCAGCCTGGCCGAATTGACGCCGCAATCGTAGTCGGCGCGGGCAAAGAGCTGGATGCCCTGGCGCTGTCAGGCGAGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGCACACgctgc                                                                                                               >  8:259837/3‑149 (MQ=255)
                         gCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGtt                                                                                                        <  5:30613/147‑1 (MQ=255)
                          cGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCg                                                                                                      <  2:312710/148‑1 (MQ=255)
                              aGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACCCGCTGCTGCCGTTCGACCtt                                                                                                 >  8:287850/1‑149 (MQ=255)
                               ggCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTc                                                                                                <  8:140213/149‑1 (MQ=255)
                                gCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCg                                                                                               <  8:252655/149‑1 (MQ=255)
                                gCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCg                                                                                               >  4:106340/1‑149 (MQ=255)
                                gCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCg                                                                                               >  8:117177/1‑149 (MQ=255)
                                 cTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGCGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCgg                                                                                              <  5:21028/149‑1 (MQ=255)
                                     ggCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCa                                                                                          >  1:151465/1‑149 (MQ=255)
                                        cGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACcg                                                                                       >  6:317201/1‑149 (MQ=255)
                                               atggCTGGACGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGGGTGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAAGACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCa                                                                                 <  2:151465/146‑1 (MQ=255)
                                                  gCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAAAGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGc                                                                             >  7:101212/1‑149 (MQ=255)
                                                      gCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCgg                                                                          <  5:211066/148‑1 (MQ=255)
                                                      gCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCgg                                                                          >  6:211066/1‑148 (MQ=255)
                                                      gCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGCGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCgg                                                                          <  7:172292/148‑1 (MQ=255)
                                                        cGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCggg                                                                         >  2:164133/1‑147 (MQ=255)
                                                        cGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCACGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGa                                                                       <  1:100737/149‑1 (MQ=255)
                                                            ttGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGgc                                                                   >  2:293776/1‑149 (MQ=255)
                                                                  aggCAATCGTTGTCGGCGCGGGCAGCGAGGGGGATGCCCTGGCGCTGTCCTGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCa                                                             <  6:84805/148‑1 (MQ=255)
                                                                   ggCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAg                                                            >  6:78673/1‑149 (MQ=255)
                                                                       aTCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCGCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGtt                                                        <  3:194475/149‑1 (MQ=255)
                                                                           ttGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCggg                                                                         >  5:277593/1‑128 (MQ=255)
                                                                           ttGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCggg                                                                         <  6:277593/128‑1 (MQ=255)
                                                                           ttGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGtg                                                    >  7:219931/1‑149 (MQ=255)
                                                                             gTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGtgct                                                  >  4:103688/1‑149 (MQ=255)
                                                                               cGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGtgctgc                                                <  5:317201/149‑1 (MQ=255)
                                                                               cGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGtgctgc                                                <  4:71082/149‑1 (MQ=255)
                                                                                ggCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGgcgc                                                                 >  7:192930/1‑131 (MQ=255)
                                                                                ggCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGgcgc                                                                 <  8:192930/131‑1 (MQ=255)
                                                                                     gggCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGgc                                          <  4:282485/149‑1 (MQ=255)
                                                                                     aggCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGgc                                          <  2:77388/148‑1 (MQ=255)
                                                                                          gCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCGCCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGaa                                     <  1:308490/149‑1 (MQ=255)
                                                                                             aGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGAc                                  >  4:5373/1‑149 (MQ=255)
                                                                                             aGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCAGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGAc                                  <  2:151071/149‑1 (MQ=255)
                                                                                              gCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGAcc                                 <  1:293776/149‑1 (MQ=255)
                                                                                                 ggATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGt                              >  6:135059/1‑149 (MQ=255)
                                                                                                    tGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCgg                           <  7:287850/149‑1 (MQ=255)
                                                                                                     gCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGc                          >  8:137058/1‑149 (MQ=255)
                                                                                                       ccTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGCCg                        <  7:117177/149‑1 (MQ=255)
                                                                                                                 tCCGGGTTGGGCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGCCGGTCCTCAg                <  7:173054/147‑1 (MQ=255)
                                                                                                                   cGGGTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGCCGGTCCTCAGGCGt            >  5:110143/1‑149 (MQ=255)
                                                                                                                     ggTTGGTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGCCGGTCCTCAGGCGTGc          <  8:52269/149‑1 (MQ=255)
                                                                                                                          gTCCAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGACGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGCCGGTCCTCAGGCGTGCCTCGc     >  8:16221/1‑149 (MQ=255)
                                                                                                                             cAGGCGGGTGTCATCAGTGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGCCGGTCCTCAGGCGTGCCTCGCCGa  <  5:78673/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                   
GTGTCGCCGGCAACCTGGCAGCCTTGCGTGAGGCTACGGCCGCTTTGCAGGCTGGCCGAATTGACGCGGCAATCGTTGTCGGCGCGGGCAGCGAGCTGGATGCCCTGGCGCTGTCCGGGTTGGTCCAGGCGGGTGTCATCAGCGCCAACGGCTCGAACACGCTGCTGCCGTTCGACCTTCGGGGCACCGGCGGCATTGCCGGGGAGGGCGCGGCAGCGTTGGTGCTGCGCCGGCGCGAAGACCTGTCGGCCGGTCCTCAGGCGTGCCTCGCCGA  >  NC_002947/3167402‑3167675

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: