Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 1,416,458 A→T 100% I216N (ATT→AAT)  oppD_4 ← Oligopeptide transport ATP‑binding protein OppD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE1,416,4580AT99.0% 353.8 / ‑3.6 96I216N (ATT→AAT) oppD_4Oligopeptide transport ATP‑binding protein OppD
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (39/56);  total (39/57)
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.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAAT  >  USA300TCH1516_ALE/1416374‑1416544
                                                                                    |                                                                                      
cGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACt                                                                                  <  1:1023164/91‑1 (MQ=255)
cGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACt                                                                                  <  1:868118/91‑1 (MQ=255)
 gTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTg                                                                                 >  1:975292/1‑91 (MQ=255)
 gTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTg                                                                                 <  1:2641389/91‑1 (MQ=255)
  ttCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGa                                                                                <  1:3362175/91‑1 (MQ=255)
   tCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGAc                                                                               >  1:3769471/1‑91 (MQ=255)
    cAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACa                                                                              >  1:3058685/1‑91 (MQ=255)
    cAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACa                                                                              >  1:1895502/1‑91 (MQ=255)
    cAGTATACGGACGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACa                                                                              <  1:3753760/91‑1 (MQ=255)
      gTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAat                                                                            <  1:120023/91‑1 (MQ=255)
      gTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAat                                                                            <  1:3626498/91‑1 (MQ=255)
         taCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATAtt                                                                         <  1:1066276/91‑1 (MQ=255)
         taCGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATAtt                                                                         <  1:3164773/91‑1 (MQ=255)
             gATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCaa                                                                     <  1:752619/91‑1 (MQ=255)
               tGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATa                                                                   >  1:2157232/1‑91 (MQ=255)
               tGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATa                                                                   >  1:1792307/1‑91 (MQ=255)
                gCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATAc                                                                  <  1:3358627/91‑1 (MQ=255)
                gCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATAc                                                                  <  1:2293735/91‑1 (MQ=255)
                gCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATAc                                                                  <  1:3360693/91‑1 (MQ=255)
                gCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATAc                                                                  <  1:4023897/91‑1 (MQ=255)
                   aaTGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATg                                                               <  1:955777/91‑1 (MQ=255)
                     tGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGAc                                                             >  1:348693/1‑91 (MQ=255)
                     tGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGAc                                                             >  1:2182564/1‑91 (MQ=255)
                      gCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACt                                                            <  1:3588295/91‑1 (MQ=255)
                      gCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACt                                                            <  1:550891/91‑1 (MQ=255)
                      gCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACt                                                            <  1:553630/91‑1 (MQ=255)
                       cTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACtt                                                           <  1:265883/91‑1 (MQ=255)
                         cATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTaa                                                         >  1:2696653/1‑91 (MQ=255)
                          aTACTTGAAATGATTAATATTGCCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAg                                                        <  1:590814/91‑1 (MQ=255)
                           tACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGt                                                       >  1:3078987/1‑91 (MQ=255)
                           tACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGt                                                       <  1:2946732/91‑1 (MQ=255)
                           tACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGt                                                       >  1:1940378/1‑91 (MQ=255)
                           tACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGt                                                       >  1:1940379/1‑91 (MQ=255)
                            aCTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGtt                                                      >  1:4043817/1‑91 (MQ=255)
                               tGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATAACTTAAGTCAtg                                                   <  1:555205/91‑1 (MQ=255)
                                 aaaTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAtgtg                                                 >  1:2047309/1‑91 (MQ=255)
                                   aTGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTAtgtgtg                                               <  1:2052524/91‑1 (MQ=255)
                                     gATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGat                                             >  1:3290182/1‑91 (MQ=255)
                                     gATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGat                                             >  1:920587/1‑91 (MQ=255)
                                     gATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGat                                             >  1:2807646/1‑91 (MQ=255)
                                      aTTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGata                                            >  1:1492897/1‑91 (MQ=255)
                                      aTTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGata                                            <  1:888857/91‑1 (MQ=255)
                                       ttAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGataa                                           >  1:2267447/1‑91 (MQ=255)
                                       ttAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGataa                                           <  1:2273079/91‑1 (MQ=255)
                                        tAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGataat                                          >  1:839739/1‑91 (MQ=255)
                                          atatTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGAtaataa                                        >  1:1957330/1‑91 (MQ=255)
                                          atatTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGAtaataa                                        >  1:105660/1‑91 (MQ=255)
                                            atTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                      >  1:2621534/1‑91 (MQ=255)
                                            atTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                      <  1:2497269/91‑1 (MQ=255)
                                            atTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                      <  1:2497209/91‑1 (MQ=255)
                                            atTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                      >  1:178466/1‑91 (MQ=255)
                                            atTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaat                                      >  1:179647/1‑91 (MQ=255)
                                             ttACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaata                                     <  1:3872929/91‑1 (MQ=255)
                                             ttACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATaataaata                                     <  1:3312964/91‑1 (MQ=255)
                                               aCCTCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAaataat                                   >  1:2277830/1‑91 (MQ=255)
                                                  tCGTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGtt                                <  1:1236480/91‑1 (MQ=255)
                                                    gTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTg                              <  1:1036173/91‑1 (MQ=255)
                                                    gTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTg                              <  1:1852195/91‑1 (MQ=255)
                                                    gTTCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTg                              <  1:310089/91‑1 (MQ=255)
                                                     ttCAATGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGa                             <  1:1410270/91‑1 (MQ=255)
                                                     ttCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGa                             <  1:30758/91‑1 (MQ=255)
                                                        aaTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCa                          <  1:4046230/91‑1 (MQ=255)
                                                        aaTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCa                          <  1:1135816/91‑1 (MQ=255)
                                                        aaTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCa                          <  1:3475855/91‑1 (MQ=255)
                                                        aaTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCa                          <  1:1135774/91‑1 (MQ=255)
                                                        aaTGATTTGACCTTCTTTTAAAACATAATTGGACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCa                          >  1:3432018/1‑91 (MQ=255)
                                                         aTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCAt                         <  1:4105112/91‑1 (MQ=255)
                                                         aTGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCAt                         <  1:1933090/91‑1 (MQ=255)
                                                          tGATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATg                        >  1:3080917/1‑91 (MQ=255)
                                                           gATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATgg                       <  1:2895532/91‑2 (MQ=255)
                                                           gATTTGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACATGtt                      >  1:217000/1‑91 (MQ=255)
                                                             tttGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTc                     <  1:3128998/91‑1 (MQ=255)
                                                             tttGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTc                     >  1:969742/1‑91 (MQ=255)
                                                               tGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTa                   <  1:3775229/91‑1 (MQ=255)
                                                               tGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTa                   <  1:2831221/91‑1 (MQ=255)
                                                               tGACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTa                   >  1:3981629/1‑91 (MQ=255)
                                                                gACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTaa                  <  1:4096083/91‑1 (MQ=255)
                                                                gACCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTaa                  <  1:4096407/91‑1 (MQ=255)
                                                                 aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAAt                 <  1:1245458/91‑1 (MQ=255)
                                                                 aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAAt                 <  1:1245459/91‑1 (MQ=255)
                                                                 aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAAt                 >  1:248291/1‑91 (MQ=255)
                                                                 aCCTTCTTTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAAt                 >  1:3058258/1‑91 (MQ=255)
                                                                  ccTTCTTTTAAAACATAATTGTACTGACAACATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATa                >  1:2761049/1‑91 (MQ=255)
                                                                     tcttTTAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACaa             <  1:4027987/91‑1 (MQ=255)
                                                                       ttttAAAACATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAata           <  1:927366/91‑1 (MQ=255)
                                                                           aaaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCt       <  1:29213/91‑1 (MQ=255)
                                                                           aaaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCt       <  1:29059/91‑1 (MQ=255)
                                                                           aaaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCt       <  1:3997122/91‑1 (MQ=255)
                                                                           aaaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCt       <  1:1701811/91‑1 (MQ=255)
                                                                            aaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTg      >  1:1688062/1‑91 (MQ=255)
                                                                            aaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTg      >  1:914056/1‑91 (MQ=255)
                                                                            aaaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTg      >  1:3335512/1‑91 (MQ=255)
                                                                             aaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGt     >  1:2532859/1‑91 (MQ=255)
                                                                             aaCATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGt     <  1:1930800/91‑1 (MQ=255)
                                                                               cATAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTTTAATACAATATGCTGtaa   <  1:1245980/91‑1 (MQ=255)
                                                                                aTAATTGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGtaat  >  1:742615/1‑91 (MQ=255)
                                                                                    |                                                                                      
CGTTCAGTATACGGATGCAAATGCTCATACTTGAAATGATTAATATTACCTCGTTCAATGATTTGACCTTCTTTTAAAACATAAATGTACTGACAATATTTCAATACATGACTTAAGTTATGTGTGATAATAAATAATGTTTGACCATGTTCTAATACAATATGCTGTAAT  >  USA300TCH1516_ALE/1416374‑1416544

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