Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,021,030 A→T 9.1% Q523L (CAG→CTG)  PP_2638 → cellulose synthase operon protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,021,0300AT 9.1% 151.4 / 4.9 55Q523L (CAG→CTG) PP_2638cellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base T (3/2);  total (28/27)
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.41e-01

CGAGGACTGGTCGACGGCCAGCCACGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGCACGCGCTGGAACAGCGTGTT  >  NC_002947/3020884‑3021175
                                                                                                                                                  |                                                                                                                                                 
cgAGGACTGGTCGACGGCCAGCCACGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTgcagc                                                                                                                                                 <  6:87/149‑1 (MQ=255)
 gAGGACTGGTCGACGGCCAGCCACGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTGTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGAATTTGCAAGGCTGGTgcggcc                                                                                                                                                >  6:171288/1‑148 (MQ=255)
        ggTCGACGGCCAGCCACGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCTGCACCAGcct                                                                                                                                         <  1:276910/149‑2 (MQ=255)
         gTCGACGGCCAGCCACGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCAc                                                                                                                                        <  7:224908/149‑1 (MQ=255)
         gTCGACGGCCAGACACGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCGTTTGCCAGGCTGGTGCAGCACCAGCCAc                                                                                                                                        >  2:170039/1‑149 (MQ=255)
               ggCCAGCCACGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATc                                                                                                                                  <  3:227721/149‑1 (MQ=255)
                 ccagccaCGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCCATGAGCCTTGGCTGGTGTATCAGCTGGCCACCACTCTGCGTCACCAAGGGCGCACTGGCGAGGCCCACGAGGCATTTGCCAGGCTGGTGCAGCCCCAGCCACCTGATCCg                                                                                                                                >  1:373948/1‑149 (MQ=255)
                     ccaCGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACTTGATCCGGCaa                                                                                                                            >  6:272449/1‑149 (MQ=255)
                       aCGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAAcg                                                                                                                          <  7:362779/149‑1 (MQ=255)
                               ggCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCGTTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGc                                                                                                                  <  1:170039/149‑1 (MQ=255)
                                   aGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAAc                                                                                                              <  2:374312/149‑1 (MQ=255)
                                    ggCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACg                                                                                                             <  3:208274/149‑1 (MQ=255)
                                     gCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACgg                                                                                                            >  4:291644/1‑149 (MQ=255)
                                     gCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACgg                                                                                                            >  8:10500/1‑149 (MQ=255)
                                      cTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGc                                                                                                           <  7:356639/149‑1 (MQ=255)
                                      cTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGc                                                                                                           <  8:327292/149‑1 (MQ=255)
                                            cccTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTc                                                                                                     >  8:217762/1‑149 (MQ=255)
                                            cccTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCAGGCAACGCGCTATGCAAACGGCCTGtg                                                                                                      >  8:205690/1‑147 (MQ=255)
                                                   ccccGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGt                                                                                              >  6:158388/1‑149 (MQ=255)
                                                       gATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAg                                                                                          <  5:230238/149‑1 (MQ=255)
                                                        aTGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGc                                                                                         <  3:133091/149‑1 (MQ=255)
                                                                  gCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATg                                                                               >  6:137221/1‑149 (MQ=255)
                                                                    tGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGcc                                                                             >  4:353437/1‑149 (MQ=255)
                                                                        gtATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTgg                                                                         >  3:113775/1‑149 (MQ=255)
                                                                         tATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTgg                                                                         <  7:217762/148‑1 (MQ=255)
                                                                          aTCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGct                                                                       >  5:14811/1‑149 (MQ=255)
                                                                             aGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGctct                                                                     <  6:394047/148‑1 (MQ=255)
                                                                               cTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGAt                                                                                <  8:312825/135‑1 (MQ=255)
                                                                               cTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGAt                                                                                >  7:312825/1‑135 (MQ=255)
                                                                                  gCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGaaa                                                                <  7:205690/148‑1 (MQ=255)
                                                                                              cGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCg                                                     <  2:247607/147‑1 (MQ=255)
                                                                                               gTCACCATGGGCGCACGGGCGAGGCCGACGCGGCATTTGCCTGGGTTGTGCAGCACCAGGCAGGTGATCCGGCATCGCTCTATGCCAACTGCCTGTTCCTGGAGTCCAGCGAGCGTGATGCCTTGGGGCTGGAAAGTCTGAGGGCGGTa                                                  <  2:373948/149‑1 (MQ=255)
                                                                                                tCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTAc                                                 >  6:69617/1‑149 (MQ=255)
                                                                                                  gcgaaGGGGGCGGTGGCGAGGCCGACGCGGCATTTGCCGGGCTGGTGCCGCGCCAGCGCGGGGCTCCGGGGGCGCGCTATGCCGGGGGCCTGTTGGTGGAGTCCCGCGGGCGGGATGCGGTGGGTGGGGGGAGGCGGAGGGCGGtgcgt                                               >  4:167634/4‑145 (MQ=255)
                                                                                                  aCCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCt                                               <  5:272449/149‑1 (MQ=255)
                                                                                                    cAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCt                                               >  5:19192/1‑147 (MQ=255)
                                                                                                     aaGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAg                                            >  3:272010/1‑149 (MQ=255)
                                                                                                          cgcACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGAc                                       >  3:163814/1‑149 (MQ=255)
                                                                                                             aCTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTg                                    <  3:167634/149‑1 (MQ=255)
                                                                                                             aCTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTg                                    >  5:287769/1‑149 (MQ=255)
                                                                                                               tGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGc                                                                                                           <  7:37775/76‑1 (MQ=255)
                                                                                                               tGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGc                                                                                                           >  8:37775/1‑76 (MQ=255)
                                                                                                               tGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTg                                    >  8:318321/1‑147 (MQ=255)
                                                                                                                 gCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAgc                                >  3:387361/1‑149 (MQ=255)
                                                                                                                   gAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAgcgc                              >  7:63069/1‑149 (MQ=255)
                                                                                                                      gCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGa                           <  4:272010/149‑1 (MQ=255)
                                                                                                                      gCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGa                           <  7:318321/149‑1 (MQ=255)
                                                                                                                      gCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGa                           >  8:360005/1‑149 (MQ=255)
                                                                                                                         gACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACAt                        >  7:18558/1‑149 (MQ=255)
                                                                                                                           cgcgGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGc                      <  8:63069/149‑1 (MQ=255)
                                                                                                                               gCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGCAcgc                  <  7:10500/149‑1 (MQ=255)
                                                                                                                                  tttGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGCACGCGCt               <  5:69617/149‑1 (MQ=255)
                                                                                                                                       gtggCTGGTGCTGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGCACGCGCTGGaa           >  1:390857/3‑148 (MQ=255)
                                                                                                                                       gtggCTGGTGCTGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGCACGCGCTGGaa           >  7:294136/3‑148 (MQ=255)
                                                                                                                                         ggCTGGTGCTGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGCGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGCCCTGGAGCGCGGACATGCACGCGCTGGAACAg        >  5:116695/1‑149 (MQ=255)
                                                                                                                                          gCTGGTGCTGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGCACGCGCTGGAACAGc       <  8:294136/149‑1 (MQ=255)
                                                                                                                                               tGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGCACGCGCTGGAACAGCGTGtt  <  4:387361/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGAGGACTGGTCGACGGCCAGCCACGTGCTGGGCCAGGCTTCGGCCCTGGCCCCCGATGAGCCTTGGCTGGTGTATCAGCTGGCCAACAGTCTGCGTCACCAAGGGCGCACTGGCGAGGCCGACGCGGCATTTGCCAGGCTGGTGCAGCACCAGCCACGTGATCCGGCAACGCGCTATGCCAACGGCCTGTTCCTGGAGTCCAGCGACCGTGATGCCTTGGCTCTGGAAAGTCTGAGGGCGGTACCTCAGGCGACCTGGAGCGCGGACATGCACGCGCTGGAACAGCGTGTT  >  NC_002947/3020884‑3021175

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