Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,803,104 A→G 20.6% pseudogene (11/426 nt) arpB → pseudogene, ankyrin repeats

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,803,1040AG20.6% 163.7 / 41.8 85pseudogene (11/426 nt)arpBpseudogene, ankyrin repeats
Reads supporting (aligned to +/- strand):  ref base A (27/39);  new base G (6/11);  total (34/51)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

ATAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCG  >  NC_000913/1803006‑1803201
                                                                                                  |                                                                                                 
ataACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAAc                                                                                                  <  1:1296527/100‑1 (MQ=255)
   aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGat                                                                                               <  1:1555215/100‑1 (MQ=255)
   aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGat                                                                                               >  2:88981/1‑100 (MQ=255)
   aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGat                                                                                               >  1:235545/1‑100 (MQ=255)
   aCAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGat                                                                                               <  1:2304960/100‑1 (MQ=255)
     aaTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGatat                                                                                             <  1:1213012/100‑1 (MQ=255)
      atgtatCGTGCATATTTCGTTGTGGCGTTAGTGATTAATCATCATCCTCAAAACCTCCATATATATATAAGGCTGTCGGATGATGAGTCAAAACGatata                                                                                            <  2:875426/95‑2 (MQ=255)
      aTTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATc                                                                                            <  2:306125/100‑1 (MQ=255)
      aTTTTTCCTGCACATTTCGTTATGGCGCTAATAATCAATCATCATCCTCAAAGCCTCCAGATATATATAAGGCTAACAGATTATGAGTCAAACCGATATc                                                                                            <  2:2106365/100‑1 (MQ=255)
       tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCa                                                                                           <  2:165432/100‑1 (MQ=255)
       tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCa                                                                                           <  2:165510/100‑1 (MQ=255)
       tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCa                                                                                           <  2:165290/100‑1 (MQ=255)
       tttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCa                                                                                           <  2:1542352/100‑1 (MQ=255)
        ttttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCAt                                                                                          <  1:88981/100‑1 (MQ=255)
         tttCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCAtt                                                                                         <  1:987017/100‑1 (MQ=255)
           tCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTAt                                                                                       <  2:900750/100‑1 (MQ=255)
             gTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCa                                                                                     <  2:423029/100‑1 (MQ=255)
              tGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAg                                                                                    <  2:2234025/100‑1 (MQ=255)
               gCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGa                                                                                   <  2:64905/100‑1 (MQ=255)
                 acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAAc                                                                                 <  1:2091795/100‑1 (MQ=255)
                 acaTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAAc                                                                                 <  1:115860/100‑1 (MQ=255)
                  caTTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACt                                                                                >  2:841486/1‑100 (MQ=255)
                    tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCa                                                                              >  2:1520941/1‑100 (MQ=255)
                    tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCa                                                                              >  2:2017997/1‑100 (MQ=255)
                    tttCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCa                                                                              >  2:871904/1‑100 (MQ=255)
                       cGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATta                                                                           <  1:572272/100‑1 (MQ=255)
                        gTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATtat                                                                          >  1:2274738/1‑100 (MQ=255)
                         ttCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATAGATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATtata                                                                         <  1:1626622/100‑1 (MQ=255)
                          tgtgtcGATAATAATTAATCATCAAACTAAAAACCTCTAGATATATATAAGGCTAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATtata                                                                         <  2:2227386/94‑1 (MQ=255)
                           cTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAg                                                                       <  2:2037772/100‑1 (MQ=255)
                            tGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTAATTATAAGt                                                                      >  2:1181279/1‑100 (MQ=255)
                             ggCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGtc                                                                     <  2:2153712/100‑1 (MQ=255)
                              gCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGtct                                                                    >  2:1082539/1‑100 (MQ=255)
                                    ataatTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCAt                                                              >  2:1798456/1‑100 (MQ=255)
                                     taattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATa                                                             >  1:1954833/1‑100 (MQ=255)
                                      aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAg                                                            >  2:1706440/1‑100 (MQ=255)
                                      aattaatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATCATAAGTCTCCTCATAg                                                            >  2:1376739/1‑100 (MQ=255)
                                         taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAtt                                                         <  2:1035263/100‑1 (MQ=255)
                                         taatCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGAt                                                          <  1:300080/99‑1 (MQ=255)
                                            tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                      >  1:1076520/1‑100 (MQ=255)
                                            tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                      <  2:788562/100‑1 (MQ=255)
                                            tcatcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGca                                                      <  1:1742261/100‑1 (MQ=255)
                                             catcatcCTCAAAACCTCCATATATATATATAGCGTACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                     <  2:1705142/100‑1 (MQ=255)
                                             catcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGcac                                                     >  2:273868/1‑100 (MQ=255)
                                               tcatcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCACCTCATAGATTGCACAt                                                   >  1:216724/1‑100 (MQ=255)
                                                catcCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATc                                                  <  2:2106461/100‑1 (MQ=255)
                                                   ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                               >  2:1936143/1‑100 (MQ=255)
                                                   ccTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGAt                                               >  2:152521/1‑100 (MQ=255)
                                                     tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                             >  2:709275/1‑100 (MQ=255)
                                                     tCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAg                                             >  1:576171/1‑100 (MQ=255)
                                                       aaaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCg                                           <  2:389354/100‑1 (MQ=255)
                                                        aaaCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                          <  1:2091299/100‑1 (MQ=255)
                                                        aaaCCACCCGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTACCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGa                                          <  1:106060/100‑1 (MQ=255)
                                                          aCCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACa                                        <  1:1099590/100‑1 (MQ=255)
                                                           ccTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACAt                                       >  1:903779/1‑100 (MQ=255)
                                                            ctacAGATATATATTAGTCGACCAGATTATGTGTCAAAACGATATCATTATTAAAATTCATTATAAGTCTCATCATATATTGCACATCGATAGCTTCATa                                      >  2:1005904/4‑100 (MQ=255)
                                                              ccAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                    <  1:1462929/100‑1 (MQ=255)
                                                              ccAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATAcc                                    <  2:2076364/100‑1 (MQ=255)
                                                              ccAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTAATAGATTGCACATCGATAGCGACATAcc                                    >  2:1448867/1‑100 (MQ=255)
                                                               cAGATATATATAAGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGTACATCGATAGCGACATACCa                                   <  2:1442188/100‑1 (MQ=255)
                                                                 gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                 <  2:1462087/100‑1 (MQ=255)
                                                                 gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                 <  1:1171083/100‑1 (MQ=255)
                                                                 gATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAac                                 <  1:711130/100‑1 (MQ=255)
                                                                  atatatatAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca                                <  1:1739516/100‑1 (MQ=255)
                                                                  atatatatAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAaca                                <  1:2327910/100‑1 (MQ=255)
                                                                   tatatataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTACTCATAGATTGCACATCGATAGCGACATACCAaca                                >  2:1090640/1‑99 (MQ=255)
                                                                      atataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCtt                            >  1:1599310/1‑100 (MQ=255)
                                                                       tataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTGTAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtc                           >  2:1343199/1‑100 (MQ=255)
                                                                        ataAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtca                          >  2:1139268/1‑100 (MQ=255)
                                                                         taAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcat                         >  1:267392/1‑100 (MQ=255)
                                                                          aaGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATGGATTGCACATCGATAGCGACATACCAACACCTtcatc                        <  2:1255638/100‑1 (MQ=255)
                                                                           aGGCGAACAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTtcatca                       <  2:332503/100‑1 (MQ=255)
                                                                            ggCGAACAGATTATGAGTCAAAACGATATCATTATCCGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCag                      >  1:1656635/1‑100 (MQ=255)
                                                                             gCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCaga                     <  1:1490542/100‑1 (MQ=255)
                                                                               gAACAGATTATGAGTCAAATCGATATTATTATCTGAACTCATTATATGTCGCTTCTTATATTGCACATCGATATCTATTTACTTTCTCCTTCATCAGAGc                   >  2:783317/1‑100 (MQ=255)
                                                                                aaCAGATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGcc                  <  1:1761252/100‑1 (MQ=255)
                                                                                  cAGATTACGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTa                >  1:1952529/1‑100 (MQ=255)
                                                                                   aGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAt               <  1:1192074/100‑1 (MQ=255)
                                                                                    gATTATGAGTCAAAGCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTAtt              <  1:2241611/100‑1 (MQ=255)
                                                                                       tATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAAt           <  2:576168/100‑1 (MQ=255)
                                                                                         tGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTACTCATAGATTGCACATCGATAGCGACATACCAACACATTCATCAGAGCCTATTAATCa         >  2:1145162/1‑100 (MQ=255)
                                                                                           aGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCa         <  2:1936201/98‑1 (MQ=255)
                                                                                            gTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAtt      >  2:2194977/1‑100 (MQ=255)
                                                                                             tCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAAttt     >  1:17719/1‑100 (MQ=255)
                                                                                              cAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTg    <  1:2145389/100‑1 (MQ=255)
                                                                                                aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACCTCGATAGCGACATACCAACACCTTCATCAGAGACTATTAATCAATTTgcg  >  1:889637/1‑100 (MQ=255)
                                                                                                aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg  <  2:2268460/100‑1 (MQ=255)
                                                                                                aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg  >  1:393927/1‑100 (MQ=255)
                                                                                                aaaCGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTgcg  >  1:858259/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
ATAACAATTTTTCGTGCACATTTCGTTCTGGCGATAATAATTAATCATCATCCTCAAAACCTCCAGATATATATAAGGCGAACAGATTATGAGTCAAAACGATATCATTATCAGAACTCATTATAAGTCTCCTCATAGATTGCACATCGATAGCGACATACCAACACCTTCATCAGAGCCTATTAATCAATTTGCG  >  NC_000913/1803006‑1803201

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