Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,541,456 A→G 22.1% I75T (ATT→ACT)  cysW ← sulfate/thiosulfate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,541,4560AG22.1% 147.5 / 41.5 81I75T (ATT→ACT) cysWsulfate/thiosulfate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base G (6/12);  total (35/46)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

CAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACG  >  NC_000913/2541362‑2541554
                                                                                              |                                                                                                  
cAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAgc                                                                                               >  2:774765/1‑100 (MQ=255)
   tatcGTTAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGcgcg                                                                                            <  1:606997/97‑1 (MQ=255)
   tAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGcgcg                                                                                            <  1:1716641/100‑1 (MQ=255)
     gCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGAAATCAGCGCGat                                                                                          >  2:206867/1‑100 (MQ=255)
       gTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGatca                                                                                        <  1:1819238/100‑1 (MQ=255)
        tCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatcat                                                                                       <  1:57775/100‑1 (MQ=255)
           gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCAc                                                                                    >  1:1068854/1‑100 (MQ=255)
           gcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCAc                                                                                    <  1:2180667/100‑1 (MQ=255)
            cagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCAcc                                                                                   >  1:887830/1‑100 (MQ=255)
             agTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGTCAATCAGCGCGATCATCAc                                                                                    <  2:168255/98‑1 (MQ=255)
              gTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGTAGAATGCCGAACACCAGGTTTACCGGTACGGCAATTAGAGCGATCATCAcc                                                                                   <  2:831096/98‑1 (MQ=255)
              gTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGt                                                                                 >  2:1064341/1‑100 (MQ=255)
               tAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTc                                                                                >  2:1046780/1‑100 (MQ=255)
                naCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCa                                                                               >  1:988624/2‑100 (MQ=255)
                  cTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGc                                                                             <  1:1173104/100‑1 (MQ=255)
                   tGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGcc                                                                            <  1:2227544/100‑1 (MQ=255)
                      cGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGa                                                                         >  1:1524709/1‑100 (MQ=255)
                      cGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGa                                                                         >  1:1992884/1‑100 (MQ=255)
                      cGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGa                                                                         >  1:1446788/1‑100 (MQ=255)
                           aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCg                                                                    <  2:1501759/100‑1 (MQ=255)
                             ggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACTCCAGGTTTACCGGTACGGCAATCAGCGCGATCATCTCCGTCAGCCAGATGGCg                                                                    <  2:37853/98‑1 (MQ=255)
                              gAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgc                                                                 <  1:2352418/100‑1 (MQ=255)
                               aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                                <  1:1240286/100‑1 (MQ=255)
                               aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                                                <  1:244211/100‑1 (MQ=255)
                                aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                                               >  1:357445/1‑100 (MQ=255)
                                aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                                               <  1:1751766/100‑1 (MQ=255)
                                 gTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagc                                                              >  1:1908461/1‑100 (MQ=255)
                                  ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagca                                                             >  2:2150023/1‑100 (MQ=255)
                                   tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCAt                                                            >  2:1492392/1‑100 (MQ=255)
                                     aaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGt                                                          <  1:2048998/100‑1 (MQ=255)
                                      aGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTc                                                         >  1:1972979/1‑100 (MQ=255)
                                       gcgcGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTcc                                                        >  2:299177/1‑100 (MQ=255)
                                        cgcgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCg                                                       <  2:2237080/100‑1 (MQ=255)
                                          cgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGa                                                     >  1:2129111/1‑100 (MQ=255)
                                           gTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACTGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGAt                                                    <  1:1646576/100‑1 (MQ=255)
                                           gTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGAt                                                    >  2:2323870/1‑100 (MQ=255)
                                               ccagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                                                >  1:1450763/1‑100 (MQ=255)
                                                 agccagGCCAGCAGAATGCCGAACAGCAGGTGTACCGGTACGGCAGTCAGCGCGATCAGCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                               <  1:812961/99‑1 (MQ=255)
                                                 agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGcc                                               <  1:797349/99‑1 (MQ=255)
                                                 agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                              <  2:1068844/100‑1 (MQ=255)
                                                  gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                                              >  1:295022/1‑99 (MQ=255)
                                                      ngccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                         >  1:2150879/2‑100 (MQ=255)
                                                      ggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTc                                         >  1:1732424/1‑100 (MQ=255)
                                                       gccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCt                                        <  2:2058896/100‑1 (MQ=255)
                                                        ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTg                                       >  1:1530998/1‑100 (MQ=255)
                                                         cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTATGt                                      >  2:714500/1‑100 (MQ=255)
                                                         cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                      <  1:1277439/100‑1 (MQ=255)
                                                         cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                      <  2:24825/100‑1 (MQ=255)
                                                          agcagAATGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                     <  2:1077475/100‑1 (MQ=255)
                                                          agcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                      <  1:962644/99‑1 (MQ=255)
                                                          agcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa                                     >  2:726148/1‑100 (MQ=255)
                                                           gcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGt                                      <  1:1591709/98‑1 (MQ=255)
                                                               aaTGCCGAACACCAGGTTTACCGGTACGGCAGTCAGCACGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAcc                                <  2:1614441/100‑1 (MQ=255)
                                                               aaTGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAAcc                                <  1:124741/100‑1 (MQ=255)
                                                                aTGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCg                               >  2:2100958/1‑100 (MQ=255)
                                                                 tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCgg                              <  2:1058174/100‑1 (MQ=255)
                                                                   ccGAACACCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                            <  2:1450747/100‑1 (MQ=255)
                                                                   ccGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                            <  2:585449/100‑1 (MQ=255)
                                                                   ccGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                            <  1:1064351/100‑1 (MQ=255)
                                                                    cGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa                            <  2:728748/99‑1 (MQ=255)
                                                                    cGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt                           <  2:293130/100‑1 (MQ=255)
                                                                       acacCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAg                        <  1:614158/100‑1 (MQ=255)
                                                                         acCAGGTTTACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGcc                      <  2:2399900/100‑1 (MQ=255)
                                                                         acCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCCGCATCAGcc                      >  1:533497/1‑100 (MQ=255)
                                                                          ccAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                     >  1:2075438/1‑100 (MQ=255)
                                                                          ccAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                     >  1:709744/1‑100 (MQ=255)
                                                                          ccAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGccc                     >  2:651401/1‑100 (MQ=255)
                                                                           cAGGTGTGCCGGTGCGGCAATCAGCGCGTTCATCACCGACAGCCAGATGGGGTGCAGCATGTCTGGATCGGCCAAATTCTGGAAAACCGGCATCAGcc                      <  2:1666690/98‑1 (MQ=255)
                                                                            aGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCt                   <  2:1726949/100‑1 (MQ=255)
                                                                               tttACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGc                >  1:1330669/1‑100 (MQ=255)
                                                                               tttACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGc                <  1:2296157/100‑1 (MQ=255)
                                                                                ttACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCt               <  1:2332048/100‑1 (MQ=255)
                                                                                 tACCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTg              <  2:750679/100‑1 (MQ=255)
                                                                                  aCCGGTACGGCAGTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGa             <  1:1434373/100‑1 (MQ=255)
                                                                                     ggTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATg          >  1:1381207/1‑100 (MQ=255)
                                                                                       tACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAATCCGGCATCAGCCCCTTGCTGAATGcc        <  2:195032/100‑1 (MQ=255)
                                                                                       tACGGCAATCAGCGCGAGCAGCACCGTCAGCCAGATGGCGTGCAGCATGTACGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGcc        <  2:2361579/100‑1 (MQ=255)
                                                                                        aCGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCt       >  1:2069471/1‑100 (MQ=255)
                                                                                        aCGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCt       >  2:1642295/1‑100 (MQ=255)
                                                                                           gCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCa    >  2:1978228/1‑100 (MQ=255)
                                                                                            cAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCAc   <  1:1154347/100‑1 (MQ=255)
                                                                                             aaTCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACg  <  2:1503457/100‑1 (MQ=255)
                                                                                              |                                                                                                  
CAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATCAGCCCCTTGCTGAATGCCTGCACG  >  NC_000913/2541362‑2541554

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: