Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,164,658 A→T I395N (ATC→AAC)  yehU ← predicted sensory kinase in two‑component system with YehT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,164,6580AT91.6% 479.7 / NA 132I395N (ATC→AAC) yehUpredicted sensory kinase in two‑component system with YehT
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (22/98);  minor base G (11/0);  total (34/98)
Fisher's exact test for biased strand distribution p-value = 6.10e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

CCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGTTTGATT  >  REL606/2164585‑2164729
                                                                          |                                                                        
ccGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTc                                                                        >  1:1835302/1‑76 (MQ=255)
 cGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCa                                                                       <  1:4784064/76‑1 (MQ=255)
   cTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCAcc                                                                     <  2:4747350/76‑1 (MQ=255)
     ttAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGc                                                                   <  2:3596028/76‑1 (MQ=255)
     ttAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGc                                                                   <  2:5591273/76‑1 (MQ=255)
       aaGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCtt                                                                 <  2:5984251/76‑1 (MQ=255)
       aaGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTCACCGCtt                                                                 >  1:2803027/1‑76 (MQ=255)
        aGTTTTTGCGGAAAAAAGTGGAAAGATCCTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCttt                                                                <  1:5089440/76‑1 (MQ=255)
        aGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCttt                                                                <  1:1258586/76‑1 (MQ=255)
         gTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTa                                                               <  2:5410285/76‑1 (MQ=255)
         gTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTa                                                               <  2:2760160/76‑1 (MQ=255)
         gTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTCACCGCTTTa                                                               >  1:1545373/1‑76 (MQ=255)
          tttttGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTaa                                                              <  2:1897441/76‑1 (MQ=255)
          tttttGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTaa                                                              <  2:3924492/76‑1 (MQ=255)
           ttttGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTGACCGCTTTTAt                                                             >  2:5487431/1‑76 (MQ=255)
           ttttGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTTACCGCCTTTAt                                                             >  1:4179468/1‑76 (MQ=255)
           ttttGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTCACCGCTTTTAt                                                             >  1:3644305/1‑76 (MQ=255)
              tGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTCACCGCCTTAATggg                                                          >  1:3041954/1‑75 (MQ=255)
               gCGGAAAAAAGTGGAAAGATGCTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGgtg                                                         <  2:1410945/76‑1 (MQ=255)
                 ggAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGtt                                                       <  2:6508722/76‑1 (MQ=255)
                 ggAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCCCTTTAATGGTGtt                                                       >  1:3295163/1‑76 (MQ=255)
                  gAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTa                                                      <  1:3582009/76‑1 (MQ=255)
                    aaaaaGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTTATGGTGTTAAg                                                    >  1:6907781/1‑76 (MQ=255)
                    aaaaaGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAg                                                    <  2:1692577/76‑1 (MQ=255)
                     aaaaGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTGAACCCTTTAAATGGGTTaaag                                                   >  2:1481000/1‑74 (MQ=255)
                      aaaGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAgcg                                                  <  1:6796646/76‑1 (MQ=255)
                       aaGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTCACCCCTTTTATTGTGTTTAgggc                                                 >  2:2932034/1‑76 (MQ=255)
                       aaGGGGAAAGATACTGCACCATCTGGCGGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAgcgc                                                 <  1:2659965/76‑1 (MQ=255)
                        aGTGGAAAGATACTGCCCCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCa                                                <  2:5137946/76‑1 (MQ=255)
                         gTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCAt                                              <  1:3901760/76‑1 (MQ=255)
                          tGGAAAGATACTGCCCCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCAtt                                             <  1:1106994/76‑1 (MQ=255)
                          tGGAAAGATACTGCCCCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCAtt                                             <  1:3766639/76‑1 (MQ=255)
                           ggAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTa                                            <  1:3393048/76‑1 (MQ=255)
                           ggAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTCACCGCTTTTATTGGGTTAAAGGCATTa                                            >  1:5428856/1‑76 (MQ=255)
                            gAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTaa                                           <  1:4478919/76‑1 (MQ=255)
                            gAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTaa                                           <  1:6930650/76‑1 (MQ=255)
                            gAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTCACCGCCTTTATTGGGTTAAACGCACTaa                                           >  2:4322591/1‑76 (MQ=255)
                             aaaGATCCTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTaaa                                          <  2:3586608/76‑1 (MQ=255)
                             aaaGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTTACAGCTTTTATTGTGTTATGGGCATTaaa                                          >  2:5243131/1‑76 (MQ=255)
                              aaGATAGTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCAGTTCACCGCTTTAATGGTGTTAAGCGCATTAAAc                                         <  2:1001573/76‑1 (MQ=255)
                               aGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTTCCCCCTTTAATGGTGTTTAGCGCATTAAACa                                        >  2:3373290/1‑76 (MQ=255)
                                gATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGGCCCCGCCTTTATGGGGGTAAACGCAGTAAACaa                                       >  2:6743655/1‑76 (MQ=255)
                                 aTACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGGGGGTCTCCGATTTTATGGTGGTTAGCGGGTTTAACaaa                                      >  1:2061020/1‑76 (MQ=255)
                                  tACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTTCACCCCTTTATGGGGGTTAACGGCTTTAAACaaa                                      >  1:6426910/1‑76 (MQ=255)
                                   aCTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTCACCGCTTTAATGGGGTTTAGCGCATTAaaaaaaaa                                    >  2:865125/1‑76 (MQ=255)
                                      gCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATg                                 <  2:1938381/76‑1 (MQ=255)
                                      gCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTTACCCCTTTTATTGGGTTAAGCGCATTAAAAAAAAAATg                                 >  1:164228/1‑76 (MQ=255)
                                       cACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGGTTCCCCCCTTTATGGTGTTTAGCGCATTAAAAAAAAAAAGa                                >  2:5536135/1‑76 (MQ=255)
                                        aCCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGCTCACCCCTTTTATGGTGGTAAGCGCATTAAACAAAAAATGAg                               >  1:2580210/1‑76 (MQ=255)
                                         ccAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAgg                              <  1:557180/76‑1 (MQ=255)
                                           aGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTTATGGTGGTAAGCGCATTAAACAAAAAATGAgggt                            >  2:535363/1‑74 (MQ=255)
                                           aGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGAt                            >  1:1269535/1‑76 (MQ=255)
                                           aGCTGGCTGGCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGAt                            <  1:528195/76‑1 (MQ=255)
                                             ctggctggCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTc                          <  1:614408/76‑1 (MQ=255)
                                             ctggctggCCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTc                          <  2:2071752/76‑1 (MQ=255)
                                             ctggctggCCTGTTCGCTGTCGCGGCGGGTTACCGCTTTTATTGTGTTAAACGGATTTAAACAAAAAATAAGGAtt                           >  2:6231444/1‑76 (MQ=255)
                                             ctggctggCCTGTTCGCTGTCGCGGCGGGTTACCGCCTTTATGGGGTTTAACGCCTTAAAAAAAAAAAGAGGGTTc                          >  2:5363427/1‑76 (MQ=255)
                                              tggctggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCa                         <  2:4525424/76‑1 (MQ=255)
                                              tggctggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCa                         >  1:6743090/1‑76 (MQ=255)
                                              tggctggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCa                         <  2:3249329/76‑1 (MQ=255)
                                              tggctggcCTGTTCGCTGTCGCGGCGGGTCACCCCCTTAATTGTGTTTAACGCCTTTAACACAAAAAAAGGAttta                         >  2:3477259/1‑74 (MQ=255)
                                               ggctggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTTATGGTGTTTAGCGCATTTAACAAAAAATGAGGGTTCAc                        >  1:3625398/1‑76 (MQ=255)
                                               ggctggcCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTCAGCGCATTAAACAAAAAATGAGGATTCAc                        >  2:4390644/1‑76 (MQ=255)
                                                gctggcCTGTTCGCTGTCGCGGCGGGTTACCGCTTTTATGGTGTTAAGCGCCTTAAACAAAAAATGAGGATTCAcc                       >  2:1237476/1‑76 (MQ=255)
                                                gctggcCTGTTCGCTGTCGCGGCGGGTCACCGCTTTTATGGGGTTAAGCGCATTAAAAAAAAAATGAGGATTccac                       >  2:4803787/1‑73 (MQ=255)
                                                 ctggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCt                      <  2:5916757/76‑1 (MQ=255)
                                                 ctggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCt                      <  2:2758652/76‑1 (MQ=255)
                                                 ctggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCt                      <  2:2832590/76‑1 (MQ=255)
                                                  tggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTg                     <  2:626633/76‑1 (MQ=255)
                                                  tggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTg                     <  1:4420406/76‑1 (MQ=255)
                                                  tggcCTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTg                     <  1:1620026/76‑1 (MQ=255)
                                                  tggcCTGTTCGCTGTCGCGGCGGTTCACCCCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTg                     >  1:4257004/1‑76 (MQ=255)
                                                  tggcCTGTTCGCTGTCGCGGCGGGTCACCGCCTTAATGGTGTTAAGCGGATTAAACAAAAAATGAGGATTCACCTg                     >  2:297556/1‑76 (MQ=255)
                                                     ccTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGc                  >  1:2553403/1‑76 (MQ=255)
                                                      cTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCg                 >  1:2124339/1‑76 (MQ=255)
                                                      cTGTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCg                 <  2:4170442/76‑1 (MQ=255)
                                                        gTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTg               <  1:2439904/76‑1 (MQ=255)
                                                        gTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTg               <  2:2536052/76‑1 (MQ=255)
                                                        gTTCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTg               <  2:2279063/76‑1 (MQ=255)
                                                        gTTCGCTGTCGCGGCGGGTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTg               >  1:2849341/1‑76 (MQ=255)
                                                         ttCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGa              >  1:3365817/1‑76 (MQ=255)
                                                         ttCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGa              <  2:2576626/76‑1 (MQ=255)
                                                         ttCGCTGTCGCGGCGGGTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGa              >  2:6869234/1‑76 (MQ=255)
                                                          tCGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGaa             <  2:5745287/76‑1 (MQ=255)
                                                           cGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTTGGCGTGAAg            <  2:1674399/76‑1 (MQ=255)
                                                           cGCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAg            >  1:1791964/1‑76 (MQ=255)
                                                            gCTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGc           <  1:1725782/76‑1 (MQ=255)
                                                             cTGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCa          >  1:4786693/1‑76 (MQ=255)
                                                              tGTCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAg         >  2:1700233/1‑76 (MQ=255)
                                                              tGTCGCGGCGGGTTACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAg         >  2:306187/1‑76 (MQ=255)
                                                               gTCGCGGCGGGTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGt        >  2:6440188/1‑76 (MQ=255)
                                                                tCGCGGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGtt       <  1:6456363/76‑1 (MQ=255)
                                                                 cgcgGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGttt      <  1:4064547/76‑1 (MQ=255)
                                                                 cgcgGCGGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGttt      <  1:3711820/76‑1 (MQ=255)
                                                                 cgcgGCGGGTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGttt      >  2:4628346/1‑76 (MQ=255)
                                                                  gcggcgGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGTTTg     <  1:2205445/76‑1 (MQ=255)
                                                                  gcggcgGTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGTTTg     <  2:1839132/76‑1 (MQ=255)
                                                                  gcggcgGGTCACCCCCTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCAACTGGGCGTGAAGCAGGTTg     >  2:4683076/1‑76 (MQ=255)
                                                                   cggcggTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGTTTGa    <  2:4136335/76‑1 (MQ=255)
                                                                    ggcggTTCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGTTTGAt   >  2:328724/1‑76 (MQ=255)
                                                                          |                                                                        
CCGCTTTAAGTTTTTGCGGAAAAAAGTGGAAAGATACTGCACCAGCTGGCTGGCCTGTTCGCTGTCGCGGCGGATCACCGCTTTAATGGTGTTAAGCGCATTAAACAAAAAATGAGGATTCACCTGGGCGTGAAGCAGTTTGATT  >  REL606/2164585‑2164729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: