Predicted mutation
evidence seq id position mutation freq annotation gene description
RA USA300TCH1516_ALE 2,387,416 G→C 5.8% V121V (GTG→GTC)  ribZ_1 → Riboflavin transporter RibZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*USA300TCH1516_ALE2,387,4160GC 5.8% 301.6 / 7.3 86V121V (GTG→GTC) ribZ_1Riboflavin transporter RibZ
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base C (2/3);  total (43/43)
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.94e-01

AATTATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGTTTAGTATTA  >  USA300TCH1516_ALE/2387326‑2387506
                                                                                          |                                                                                          
aaTTATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTtgtg                                                                                            >  1:204525/1‑91 (MQ=255)
 aTTATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGa                                                                                           <  1:1092792/91‑1 (MQ=255)
 aTTATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGa                                                                                           <  1:2051594/91‑1 (MQ=255)
 aTTATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGa                                                                                           >  1:506553/1‑91 (MQ=255)
  ttATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGAc                                                                                          <  1:1132296/91‑1 (MQ=255)
  ttATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGAc                                                                                          <  1:18776/91‑1 (MQ=255)
  ttATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGAc                                                                                          <  1:18875/91‑1 (MQ=255)
   tATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTcaca                                                                                         <  1:1301285/91‑1 (MQ=255)
    aTCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTCACaa                                                                                        <  1:1937549/91‑1 (MQ=255)
     tCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAAc                                                                                       >  1:3101950/1‑91 (MQ=255)
     tCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTCACAAc                                                                                       <  1:2080801/91‑1 (MQ=255)
      cACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACt                                                                                      >  1:1341721/1‑91 (MQ=255)
        cGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGc                                                                                    <  1:547622/91‑1 (MQ=255)
        cGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGc                                                                                    >  1:415376/1‑91 (MQ=255)
          cAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCaa                                                                                  >  1:2494220/1‑91 (MQ=255)
          cAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTCACAACTGCag                                                                                  <  1:884601/91‑2 (MQ=255)
              tACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAgg                                                                              >  1:898304/1‑91 (MQ=255)
                cTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTa                                                                            >  1:1530040/1‑91 (MQ=255)
                cTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTa                                                                            >  1:2663118/1‑91 (MQ=255)
                 tttattttatTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAGACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTAc                                                                           <  1:2673021/91‑1 (MQ=255)
                 tttattttatTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTAc                                                                           >  1:2725367/1‑91 (MQ=255)
                 tttattttatTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTAc                                                                           <  1:1124247/91‑1 (MQ=255)
                 tttattttatTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTAc                                                                           <  1:1124246/91‑1 (MQ=255)
                    attttattGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTAt                                                                        >  1:1903846/1‑91 (MQ=255)
                            gATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCAt                                                                <  1:1090392/91‑1 (MQ=255)
                                  tCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTa                                                          >  1:3520767/1‑91 (MQ=255)
                                   cTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTAc                                                         <  1:2334943/91‑1 (MQ=255)
                                   cTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTAc                                                         <  1:131601/91‑1 (MQ=255)
                                    tttAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTAcc                                                        <  1:4024502/91‑1 (MQ=255)
                                     ttAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTAccc                                                       <  1:2386192/91‑1 (MQ=255)
                                     ttAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTAccc                                                       >  1:1708207/1‑91 (MQ=255)
                                      tAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTAcccc                                                      >  1:528846/1‑91 (MQ=255)
                                      tAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTAcccc                                                      >  1:1935812/1‑91 (MQ=255)
                                      tAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTAcccc                                                      <  1:2752254/91‑1 (MQ=255)
                                          ttttttAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGtt                                                  >  1:1206333/1‑91 (MQ=255)
                                            ttttAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTaa                                                >  1:658503/1‑91 (MQ=255)
                                            ttttAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTaa                                                >  1:658351/1‑91 (MQ=255)
                                             tttAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAAt                                               >  1:2021840/1‑91 (MQ=255)
                                              ttAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATa                                              >  1:1995397/1‑91 (MQ=255)
                                              ttAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATa                                              >  1:1616014/1‑91 (MQ=255)
                                                aaTCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGa                                            <  1:2902042/91‑1 (MQ=255)
                                                   cTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGaaaa                                         >  1:787945/1‑91 (MQ=255)
                                                     ttGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAACAAg                                       >  1:2965101/1‑91 (MQ=255)
                                                     ttGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAg                                       >  1:2965492/1‑91 (MQ=255)
                                                     ttGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAg                                       >  1:2218025/1‑91 (MQ=255)
                                                     ttGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAg                                       <  1:2360238/91‑1 (MQ=255)
                                                      tGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATCCGTTACCCCTGTTAATAGAAAAAGt                                      >  1:3847593/1‑91 (MQ=255)
                                                      tGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGt                                      >  1:925489/1‑91 (MQ=255)
                                                      tGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGt                                      >  1:3225797/1‑91 (MQ=255)
                                                      tGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGt                                      <  1:3324740/91‑1 (MQ=255)
                                                       gTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTg                                     <  1:1637019/91‑1 (MQ=255)
                                                        tACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGa                                    >  1:357389/1‑91 (MQ=255)
                                                             ctcATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCa                               >  1:263653/1‑91 (MQ=255)
                                                               cATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATc                             <  1:3817579/91‑1 (MQ=255)
                                                                aTAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCa                            >  1:133844/1‑91 (MQ=255)
                                                                 tAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAg                           <  1:808911/91‑1 (MQ=255)
                                                                 tAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAg                           <  1:2384069/91‑1 (MQ=255)
                                                                 tAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAg                           <  1:2384089/91‑1 (MQ=255)
                                                                 tAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAg                           <  1:810258/91‑1 (MQ=255)
                                                                 tAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAg                           <  1:144119/91‑1 (MQ=255)
                                                                 tAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAg                           <  1:808887/91‑1 (MQ=255)
                                                                    aCGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGttt                        >  1:2752405/1‑91 (MQ=255)
                                                                     cGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTa                       >  1:1705655/1‑91 (MQ=255)
                                                                     cGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTa                       >  1:1705686/1‑91 (MQ=255)
                                                                      ggTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGCGAAGGCATCAGTTTAt                      >  1:2803394/1‑91 (MQ=255)
                                                                         aTCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTAtttt                   >  1:2836228/1‑91 (MQ=255)
                                                                         aTCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTAtttt                   >  1:1857339/1‑91 (MQ=255)
                                                                           cGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTcc                 <  1:1179939/91‑1 (MQ=255)
                                                                           cGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTcc                 <  1:609146/91‑1 (MQ=255)
                                                                           cGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTcc                 <  1:2265068/91‑1 (MQ=255)
                                                                           cGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTcc                 >  1:2450143/1‑91 (MQ=255)
                                                                           cGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTcc                 <  1:2711012/91‑1 (MQ=255)
                                                                            gCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCt                >  1:2214236/1‑91 (MQ=255)
                                                                             cAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTTCCCCTGTTAATGGAAAAAGTGAAGGCATCAGTTTATTTTCCtt               <  1:1674979/91‑1 (MQ=255)
                                                                             cAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCtt               <  1:3720897/91‑1 (MQ=255)
                                                                             cAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCtt               <  1:3336905/91‑1 (MQ=255)
                                                                             cAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCtt               <  1:2986356/91‑1 (MQ=255)
                                                                             cAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCtt               <  1:2986302/91‑1 (MQ=255)
                                                                               aCTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTaa             >  1:1093058/1‑91 (MQ=255)
                                                                                 tGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGt           <  1:963392/91‑1 (MQ=255)
                                                                                 tGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGt           <  1:963544/91‑1 (MQ=255)
                                                                                 tGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGt           <  1:2223851/91‑1 (MQ=255)
                                                                                  gCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGtt          <  1:544230/91‑1 (MQ=255)
                                                                                   cAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGttt         >  1:519538/1‑91 (MQ=255)
                                                                                     gTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGTTTAg       <  1:787143/91‑1 (MQ=255)
                                                                                     gTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGTTTAg       >  1:2271758/1‑91 (MQ=255)
                                                                                      ttGTCACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAACGCATCAGTTTATTTTCCTTATGTTTaat      >  1:3323307/1‑89 (MQ=255)
                                                                                      ttGTCACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCCGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGTTTAGt      >  1:3899708/1‑91 (MQ=255)
                                                                                         tgACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGTTTAGTAtt   <  1:3882614/91‑1 (MQ=255)
                                                                                          gACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGTTTAGTATTa  >  1:2668911/1‑91 (MQ=255)
                                                                                          |                                                                                          
AATTATCACGCAATTACTTTATTTTATTGATGGATCTTTAGCTTTTTTAATCTTTGTACGTCTCATAAACGGTATCGCAACTGCAGTTGTGACAACTGCAACAGGTACTATTGCGGCATACGTTACCCCTGTTAATAGAAAAAGTGAAGGCATCAGTTTATTTTCCTTAAGTTTAGTATTA  >  USA300TCH1516_ALE/2387326‑2387506

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