Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,149,6930TG37.0% 140.4 / 18.9 81G132G (GGT→GGGflgGflagellar component of cell‑distal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base T (11/40);  new base G (29/1);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 7.01e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 109 total aligned reads displayed.

CAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGCGTGGTC  >  REL606/1149620‑1149768
                                                                          |                                                                            
cAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCtggtgg                                                                            <  2:518647/76‑1 (MQ=255)
cAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCtggggg                                                                            >  2:4006028/1‑76 (MQ=255)
cAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCtggggg                                                                            >  1:4731047/1‑76 (MQ=255)
cAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCtggggg                                                                            >  2:122551/1‑76 (MQ=255)
cAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCtggggg                                                                            >  2:1434514/1‑76 (MQ=255)
 aGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCtgggggg                                                                           >  1:1037724/1‑75 (MQ=255)
 aGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCtgggggg                                                                           >  1:1679378/1‑75 (MQ=255)
 aGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTggtggt                                                                           <  2:5001080/76‑1 (MQ=255)
   aTGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGgtt                                                                         >  2:77498/1‑76 (MQ=255)
    tGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGgttt                                                                        >  1:4069669/1‑76 (MQ=255)
    tGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGgttt                                                                        >  1:4989852/1‑76 (MQ=255)
      gTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCa                                                                      <  2:4731047/76‑1 (MQ=255)
      gTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCa                                                                      >  1:3912258/1‑76 (MQ=255)
       ttCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGCGACGGCTGGGGGTTTTCAg                                                                     >  1:4948842/1‑76 (MQ=255)
       ttCATCAGCCTATACCCGTGACGGCTCTTTCCAGGGGGATCAGAACGGGCAGCTGGGGCCGGCTGGGTGGtttcaa                                                                      >  2:1026682/1‑73 (MQ=255)
        tcatcaGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTTAgg                                                                    >  2:4641435/1‑76 (MQ=255)
           tcaGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGTTTTCAGGTGc                                                                 >  1:642626/1‑76 (MQ=255)
           tcaGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGTTTTCAGGGGc                                                                 >  1:1697387/1‑76 (MQ=255)
           tcaGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTTAGGTGc                                                                 >  1:1372516/1‑76 (MQ=255)
            caGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGca                                                                <  1:4414494/76‑1 (MQ=255)
             aGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTggcg                                                               >  2:2767326/1‑73 (MQ=255)
              gCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagc                                                              <  1:294806/76‑1 (MQ=255)
              gCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGGTTTCAGGTGcagc                                                              >  2:4097757/1‑76 (MQ=255)
                cTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGcagccc                                                            >  2:1014392/1‑75 (MQ=255)
                cTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGGGACGGCTGGTGGGTTTTAAGTGgagcgc                                                            >  1:3599983/1‑74 (MQ=255)
                  ataCCCGGGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGcagccagc                                                          <  2:1977099/76‑1 (MQ=255)
                   taCCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCg                                                         <  2:3239851/76‑1 (MQ=255)
                   taCCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCg                                                         <  2:4905528/76‑1 (MQ=255)
                   taCCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCg                                                         <  2:3590633/76‑1 (MQ=255)
                   taCCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGGGGAGGCAGCg                                                         >  2:2247426/1‑76 (MQ=255)
                   taCCAGTGGTGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCg                                                         <  2:3385541/76‑1 (MQ=255)
                      ccGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGGTTTCAGGTGCAGCCAGCGGTc                                                     >  1:3422496/1‑76 (MQ=255)
                      ccGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCCGCGATc                                                     >  1:3698080/1‑76 (MQ=255)
                       cGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCa                                                    <  1:1929208/76‑1 (MQ=255)
                       cGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCTGGGGCAGCCAGCGATCa                                                    >  1:3808149/1‑76 (MQ=255)
                        gTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGCAGGCTGGGGGGTTTTGGGCGCAGGCCGCGATcca                                                   >  1:3243639/1‑74 (MQ=255)
                          gACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCa                                                 <  1:4285566/76‑1 (MQ=255)
                           aCGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCAGCGATCAccct                                                >  1:1181225/1‑74 (MQ=255)
                           aCGCCTCTTCCAGGGGGGATCAGACGGGGCAGCTGGGGACGGCTGGTGGTTTTCGGGTGCGGCAAGCGTTCCCAtt                                               <  2:2963594/76‑1 (MQ=255)
                            cgtctctTTCTAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAtt                                               <  2:2963595/73‑1 (MQ=255)
                            cGGTTCTTTCGAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAtt                                               <  2:4344663/76‑1 (MQ=255)
                            cGGCTGTTTCCAGGAGGACCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAtt                                               <  1:1722511/76‑1 (MQ=255)
                            cGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCAtt                                               <  2:1046629/76‑1 (MQ=255)
                            cGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTTAGGGGCAGCCAGGGATCACCact                                               >  2:4738750/1‑74 (MQ=255)
                             ggCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTGTTCAGGTGCAGCCCGCGCTCACCATTc                                              >  2:4109293/1‑76 (MQ=255)
                              gCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCAGCGAACACCATTcc                                             >  1:3054420/1‑76 (MQ=255)
                               ctctTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCg                                            >  1:2309807/1‑76 (MQ=255)
                               ctctTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGTTTTTAGGTGCAGCCAGCGATCACCCTTACg                                            >  1:1742539/1‑76 (MQ=255)
                               ctctTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCCGCGGTCACCATTCCg                                            >  1:811039/1‑76 (MQ=255)
                               ctctTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCAGCGATCACCATTTCg                                            >  2:3980748/1‑76 (MQ=255)
                                tctTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCgg                                           <  1:3816052/76‑1 (MQ=255)
                                tctTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGGGCAGCCAGCGATCACCATTCCgg                                           >  1:347977/1‑76 (MQ=255)
                                 ctTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCAGCGGTCACCATTCCGGc                                          >  2:800678/1‑76 (MQ=255)
                                 ctTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGGGCAGCCAGCGATCACCCTTCCGGc                                          >  2:2984102/1‑76 (MQ=255)
                                  tttCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCAGCGATCACCCTTCCGGCg                                         >  1:3058055/1‑76 (MQ=255)
                                  tttCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGGGGAGCCAGCGATCCCCATTCCGGCg                                         >  1:3940155/1‑76 (MQ=255)
                                   ttCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCCGGGGCAGACAGCGATCCACCTTTCCgggc                                         >  1:1116803/1‑74 (MQ=255)
                                     ccAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAAt                                      <  2:4839765/76‑1 (MQ=255)
                                     ccAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAAt                                      <  2:2108780/76‑1 (MQ=255)
                                     ccAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAAt                                      <  1:3601237/76‑1 (MQ=255)
                                       aGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCCGCGAATGc                                    >  1:1598029/1‑76 (MQ=255)
                                       aGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCAGCGATCACCCTTCCGGCGAATGc                                    >  1:2997303/1‑76 (MQ=255)
                                       aGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCAGCGATCACCATTCCCGCGAATGc                                    >  2:4765214/1‑76 (MQ=255)
                                        ggTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCg                                   >  1:2547472/1‑76 (MQ=255)
                                          tGGATCAGGACGGGCAGCTGGTGACGGCTGGTGTTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGtt                                 <  1:1224737/76‑1 (MQ=255)
                                          tGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGtt                                 <  2:4159084/76‑1 (MQ=255)
                                           ggATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCCACATTTCGGCGGAAGCGTTa                                >  2:3855244/1‑76 (MQ=255)
                                           ggATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTa                                <  2:1702152/76‑1 (MQ=255)
                                           ggATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTa                                <  2:2904046/76‑1 (MQ=255)
                                           ggATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTa                                <  2:558179/76‑1 (MQ=255)
                                              tCAGAACGGGCAGCTGGTGACGGCTGGGGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGt                             >  1:2497844/1‑76 (MQ=255)
                                               cAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTa                            <  1:517496/76‑1 (MQ=255)
                                               cAGAACGGGCAGCTGGTGACGGCTGGTGGGTTTCAGGTGCAGCCAGCGATCCCCATTCCCGCGAATGCGTTAAGTa                            >  1:3499014/1‑76 (MQ=255)
                                               cAGAACGGGCAGCTGGTGACGGCTGGGGGTTTTCAGGGGCAGCCAGCGATCACCATTCCGGCGGATGCGTTAAGTa                            >  1:2898013/1‑76 (MQ=255)
                                                aGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTAt                           <  2:299406/76‑1 (MQ=255)
                                                aGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTAt                           <  2:1216179/76‑1 (MQ=255)
                                                aGAACGGGCAGCTGGTGACGGCTGGGGGTTTTCCGGGGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTAt                           >  1:4333641/1‑76 (MQ=255)
                                                 gAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATc                          <  2:4236732/76‑1 (MQ=255)
                                                  aaCGGGCAGCTGGTGACGGCTGGGGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCa                         >  1:679379/1‑76 (MQ=255)
                                                   aCGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCAc                        <  1:2635658/76‑1 (MQ=255)
                                                     gggCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCa                      <  1:144573/76‑1 (MQ=255)
                                                     gggCAGCTGGTGACGGCTGGGGGGTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCa                      <  2:4295560/76‑1 (MQ=255)
                                                        cAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCg                   <  2:3184916/76‑1 (MQ=255)
                                                        cAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCg                   <  1:4187460/76‑1 (MQ=255)
                                                        cAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCg                   <  1:4976743/76‑1 (MQ=255)
                                                         aGCTGGTGACGGCTGGGGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCgg                  >  1:4850627/1‑76 (MQ=255)
                                                           cTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTc                <  2:2607221/76‑1 (MQ=255)
                                                             ggTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGt              <  2:201551/76‑1 (MQ=255)
                                                              gTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTg             >  1:5051111/1‑76 (MQ=255)
                                                                gACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGAt           <  1:4889188/76‑1 (MQ=255)
                                                                 aCGGCTGGGGGGTTTTAGGGGCAGCCCGCGATCACCATTTCGGCGAATGCGTTAAGTAACACCATCGGTCGGGATg          >  2:2064162/1‑76 (MQ=255)
                                                                  cGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATgg         <  2:1045937/76‑1 (MQ=255)
                                                                   ggCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGc        <  1:3655396/76‑1 (MQ=255)
                                                                    gCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGCg       <  1:833630/76‑1 (MQ=255)
                                                                     cTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGCGt      >  1:2901686/1‑76 (MQ=255)
                                                                     cTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGCGt      >  1:3174746/1‑76 (MQ=255)
                                                                     cTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGCGt      <  1:2417830/76‑1 (MQ=255)
                                                                      tggtggTTTTCAGGTGCAGCTAGCGATCGCCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGGGATGGCGTg     <  2:4466716/76‑1 (MQ=255)
                                                                      tggtggTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGCGTg     <  2:696194/76‑1 (MQ=255)
                                                                          tggtTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGCGTGGTc  >  2:1020046/1‑76 (MQ=255)
                                                                          |                                                                            
CAGATGGTTCATCAGCCTATACCCGTGACGGCTCTTTCCAGGTGGATCAGAACGGGCAGCTGGTGACGGCTGGTGGTTTTCAGGTGCAGCCAGCGATCACCATTCCGGCGAATGCGTTAAGTATCACCATCGGTCGTGATGGCGTGGTC  >  REL606/1149620‑1149768

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