Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,424,538 Δ1 bp 6.6% coding (89/915 nt) PP_4764 → histone deacetylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,424,5380A. 6.6% 182.6 / 6.5 61coding (89/915 nt)PP_4764histone deacetylase family protein
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base . (2/2);  total (29/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.27e-01

CCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCC  >  NC_002947/5424393‑5424684
                                                                                                                                                 |                                                                                                                                                   
ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAATATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGac                                                                                                                                                 <  2:330033/150‑1 (MQ=255)
ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCGGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGac                                                                                                                                                 >  1:529092/1‑150 (MQ=255)
ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGa                                                                                                                                                  >  2:240693/1‑149 (MQ=255)
ccATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGa                                                                                                                                                  <  1:240693/149‑1 (MQ=255)
   ttGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGTCTGCTGCACGaccat                                                                                                                                              <  2:189710/150‑2 (MQ=255)
             ccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCa                                                                                                                                   >  1:1556309/1‑150 (MQ=255)
             ccGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCa                                                                                                                                   >  1:1460487/1‑150 (MQ=255)
               gTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGaccac                                                                                                                                              <  1:736350/138‑1 (MQ=255)
               gTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGaccac                                                                                                                                              >  2:736350/1‑138 (MQ=255)
                  aTGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCggg                                                                                                                              <  2:1231742/150‑1 (MQ=255)
                    ggCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCt                                                                                                                            <  1:1091776/150‑1 (MQ=255)
                             aGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCCGACCAGGTGGTCgtgca                                                                                                                                   >  2:1268439/1‑130 (MQ=255)
                             aGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCCGACCAGGTGGTCgtgca                                                                                                                                   <  1:1268439/134‑5 (MQ=255)
                              gCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGa                                                                                                                   <  2:519703/149‑1 (MQ=255)
                              gCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGa                                                                                                                   >  1:519703/1‑149 (MQ=255)
                                 cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCggg                                                                                                                              <  1:1288394/135‑1 (MQ=255)
                                 cTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCggg                                                                                                                              >  2:1288394/1‑135 (MQ=255)
                                     cccTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGcc                                                                                                            <  2:242958/149‑1 (MQ=255)
                                     cccTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGcc                                                                                                            >  1:242958/1‑149 (MQ=255)
                                      ccTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCCGCTGCCGACCAGGTGGTCGTGCAGCAGGcggaac                                                                                                                        >  1:1060845/1‑131 (MQ=37)
                                      ccTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCCGCTGCCGACCAGGTGGTCGTGCAGCAGGCGGaac                                                                                                                        <  2:1060845/136‑3 (MQ=37)
                                         gATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgct                                                                                                       <  2:1299009/150‑1 (MQ=255)
                                           ttGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCtgctgc                                                                                                     >  2:806827/1‑150 (MQ=255)
                                                  cccATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGa                                                                                              <  2:1556309/150‑1 (MQ=255)
                                                    catccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACa                                                                                            >  1:1551648/1‑150 (MQ=255)
                                                    catccatGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCAATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGCCCAGGCCCTGCTGCGCCCGTAcc                                                                                            >  2:366628/1‑149 (MQ=255)
                                                     atccatGCCGTGGCCGCTGATCTACCACGAACATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACAt                                                                                           <  1:340338/150‑1 (MQ=255)
                                                         atGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGc                                                                                       <  1:975535/150‑1 (MQ=255)
                                                                   cGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACa                                                                             >  2:337967/1‑150 (MQ=255)
                                                                          ctcccccGAAGATTACAGCCCGGAGTGCCCCGCTGAACACCGCTTCCCGATGGACAACTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGc                                                                      <  1:77669/144‑1 (MQ=255)
                                                                           tACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACACCTGCCCCATCGACATCCTCGcc                                                                     <  2:120980/150‑1 (MQ=255)
                                                                                     aTTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATg                                                           >  1:1603632/1‑150 (MQ=255)
                                                                                        aCAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCCGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGAcc                                                        >  2:82111/1‑150 (MQ=255)
                                                                                            cccGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAAc                                                                                 <  1:1065354/121‑1 (MQ=255)
                                                                                            cccGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAAc                                                                                 >  2:1065354/1‑121 (MQ=255)
                                                                                               ggAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgc                                                                                                        <  1:1256752/95‑1 (MQ=255)
                                                                                               ggAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCctgc                                                                                                        >  2:1256752/1‑95 (MQ=255)
                                                                                                     ccccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGa                                           <  2:529092/150‑1 (MQ=255)
                                                                                                      cccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAg                                          >  1:492297/1‑150 (MQ=255)
                                                                                                      cccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAg                                          <  2:1551648/150‑1 (MQ=255)
                                                                                                       ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                         >  2:1637161/1‑150 (MQ=255)
                                                                                                       ccGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgc                                         <  1:337967/150‑1 (MQ=255)
                                                                                                         gCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgcgc                                       <  2:310076/150‑1 (MQ=255)
                                                                                                         gCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAgcgc                                       >  2:523436/1‑150 (MQ=255)
                                                                                                           tGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTa                                     <  2:1460487/150‑1 (MQ=255)
                                                                                                              acacCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCCGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACAt                                  <  1:82111/150‑1 (MQ=255)
                                                                                                                acCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACAt                                                                                           >  2:909105/1‑91 (MQ=255)
                                                                                                                acCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACAt                                                                                           <  1:909105/91‑1 (MQ=255)
                                                                                                                acCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGa                                >  1:1303422/1‑150 (MQ=255)
                                                                                                                       ccccATGGACAAGTTCCGCCTGCTGCACGACCCCCTCGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCACGACATCTGCCCCAACCACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGa                         <  1:658108/150‑1 (MQ=255)
                                                                                                                        ccGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAg                        <  2:1303422/150‑1 (MQ=255)
                                                                                                                               aCAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGa                                           <  2:972896/124‑1 (MQ=255)
                                                                                                                               aCAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGa                                           >  1:972896/1‑124 (MQ=255)
                                                                                                                                cAAGTTCCGCCTGCTGCACGACCACCTGGTCAGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTCTcgcg                >  2:898198/1‑150 (MQ=255)
                                                                                                                                 aaGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTcgcgc               <  2:1050124/150‑1 (MQ=255)
                                                                                                                                    ttCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAg            >  1:1369562/1‑150 (MQ=255)
                                                                                                                                       cGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAAc                                                                                 <  1:1319891/78‑1 (MQ=255)
                                                                                                                                       cGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAAc                                                                                 >  2:1319891/1‑78 (MQ=255)
                                                                                                                                       cGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGAc         <  2:1369562/150‑1 (MQ=255)
                                                                                                                                              tgcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGcgcc  >  1:1567381/1‑150 (MQ=255)
                                                                                                                                              tgcACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCAAGGACCAGcgcc  >  1:335369/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCATTGCAGTCACCCGTCATGGCAAGATCAGCACTGGCCCTGATTGCGACCCCATCCATGCCGTTGCCGCTGATCTACCACGAAGATTACAGCCCGGAGTTCCCCGCTGAACACCGCTTCCCGATGGACAAGTTCCGCCTGCTGCACGACCACCTGGTCGGCAGCGGGCTGACCACCGACCAGGCCCTGCTGCGCCCGGACATCTGCCCCAACGACATCCTCGCCCTGGCCCATGACCGCAGCTACATCGAGCGCTACATGAACGGCGAGCTGTCGCGCGAGGACCAGCGCC  >  NC_002947/5424393‑5424684

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: