Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,353,121 C→T *310* (TAG→TAA yhcC ← putative Fe‑S oxidoreductase, Radical SAM superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,353,1210CT100.0% 233.1 / NA 72*310* (TAG→TAAyhcCputative Fe‑S oxidoreductase, Radical SAM superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/34);  total (38/34)

GGTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACG  >  NC_000913/3352986‑3353255
                                                                                                                                       |                                                                                                                                      
ggTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTAc                                                                                                                                     >  1:992638/1‑139 (MQ=255)
 gTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    <  2:547277/139‑1 (MQ=255)
           taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                          >  1:1012776/1‑139 (MQ=255)
           taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                          <  2:992638/139‑1 (MQ=255)
           taCTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAgg                                                                                                                          <  1:1142119/139‑1 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCg                                                                                                                                 <  2:127727/127‑1 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCg                                                                                                                                 >  1:127727/1‑127 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  1:250989/1‑139 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  2:1107924/1‑139 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGCTGGAGGTAGcc                                                                                                                     >  2:43801/1‑139 (MQ=255)
                cgcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGCTATAAATTATGTGAGAAATGCAGAATTCCCCTCAACAACTTACTCTGTTGGAGGTAGcc                                                                                                                     >  2:24472/1‑139 (MQ=255)
                 gcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           >  2:723210/1‑132 (MQ=255)
                 gcCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           <  1:723210/132‑1 (MQ=255)
                   cagcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGtt                                                                                                                               >  1:1170569/1‑126 (MQ=255)
                   cagcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGtt                                                                                                                               <  2:1170569/126‑1 (MQ=255)
                     gcagACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGa                                                                                                                <  2:250989/139‑1 (MQ=255)
                                    tCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 <  1:1307153/139‑1 (MQ=255)
                                               ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGa                                                                                                                            >  1:627885/1‑101 (MQ=255)
                                               ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGa                                                                                                                            <  2:627885/101‑1 (MQ=255)
                                               ttCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGcc                                                                                      >  2:1173257/1‑139 (MQ=255)
                                                tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    <  2:1230228/92‑1 (MQ=255)
                                                tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACt                                                                                                                                    >  1:1230228/1‑92 (MQ=255)
                                                tCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCa                                                                                     <  1:1107924/139‑1 (MQ=255)
                                                          cAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGttcattca                                                                           <  1:43801/139‑1 (MQ=255)
                                                             ccTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGAt                                                                        >  2:1234488/1‑139 (MQ=255)
                                                                         ataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    <  2:90864/99‑1 (MQ=255)
                                                                         ataCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    >  1:90864/1‑99 (MQ=255)
                                                                           aCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGt                                                                                                                                <  1:212202/69‑1 (MQ=255)
                                                                           aCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGt                                                                                                                                >  2:212202/1‑69 (MQ=255)
                                                                           aCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 <  1:1251703/116‑1 (MQ=255)
                                                                           aCCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGtt                                                                                 >  2:1251703/1‑116 (MQ=255)
                                                                               aaTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCa                                                      >  1:562825/1‑139 (MQ=255)
                                                                                aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    <  2:1180557/92‑1 (MQ=255)
                                                                                aTTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGcc                                                                                                    >  1:1180557/1‑92 (MQ=255)
                                                                                 ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGtt                                                                                                                               <  1:740996/64‑1 (MQ=255)
                                                                                 ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGtt                                                                                                                               >  2:740996/1‑64 (MQ=255)
                                                                                 ttCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCAtt                                                    >  1:1195314/1‑139 (MQ=255)
                                                                                  tCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTc                                                   <  2:1012776/139‑1 (MQ=255)
                                                                                  tCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTc                                                   <  1:1115304/139‑1 (MQ=255)
                                                                                   cTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTcc                                                  <  1:1173257/139‑1 (MQ=255)
                                                                                      aCAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGt                                               >  2:445110/1‑139 (MQ=255)
                                                                                           tACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATc                                          >  1:649694/1‑139 (MQ=255)
                                                                                             cGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGa                                        >  1:960835/1‑139 (MQ=255)
                                                                                                  taAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCAt                                                                                    <  1:1267144/90‑1 (MQ=255)
                                                                                                  taAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCAt                                                                                    >  2:1267144/1‑90 (MQ=255)
                                                                                                    aaTTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     >  1:888473/1‑55 (MQ=255)
                                                                                                    aaTTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGcc                                                                                                                     <  2:888473/55‑1 (MQ=255)
                                                                                                     aTTGTGTGAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGca                                <  1:1234488/139‑1 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTcc                                                                                                           <  1:907360/61‑1 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTcc                                                                                                           >  2:907360/1‑61 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGAc                                                        >  1:1055743/1‑112 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGAc                                                        <  2:1055743/112‑1 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcacca                             <  1:445110/139‑1 (MQ=255)
                                                                                                        gtgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcacca                             <  2:562825/139‑1 (MQ=255)
                                                                                                         tgtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGa                                                                   >  1:964122/1‑100 (MQ=255)
                                                                                                         tgtgAGAAATGCGGAATTCCCCTAAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGa                                                                   <  2:964122/100‑1 (MQ=255)
                                                                                                          gtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           >  1:520062/1‑43 (MQ=255)
                                                                                                          gtgAGAAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAg                                                                                                                           <  2:520062/43‑1 (MQ=255)
                                                                                                             agaAATGCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGTCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACGGCg                        <  1:1258291/139‑1 (MQ=255)
                                                                                                                   gCGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACGGCGGagcaa                  >  1:782164/1‑139 (MQ=255)
                                                                                                                    cGGAATTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACGGCGGagcaag                 >  1:973498/1‑139 (MQ=255)
                                                                                                                       aaTTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcaccacggc                         >  2:609359/1‑125 (MQ=255)
                                                                                                                       aaTTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGcaccacggc                         <  1:609359/128‑4 (MQ=255)
                                                                                                                       aaTTCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACGGCGGAGCAagcag              >  2:396596/1‑139 (MQ=255)
                                                                                                                         ttCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACGGCGGAGCAAGCAGCg            >  2:651617/1‑139 (MQ=255)
                                                                                                                          tCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 >  2:211586/1‑53 (MQ=255)
                                                                                                                          tCCCCTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGAt                                                                                                 <  1:211586/53‑1 (MQ=255)
                                                                                                                             ccTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACGGCGGAGCAAGCAGCGTCgg        >  1:374751/1‑139 (MQ=255)
                                                                                                                              cTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTcc                                                  >  1:1135090/1‑96 (MQ=255)
                                                                                                                              cTCAACGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTcc                                                  <  2:1135090/96‑1 (MQ=255)
                                                                                                                                  aCGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCCTTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACGGCGGAGCAAGCAGCGTCGGGCGAc   <  2:960835/139‑1 (MQ=255)
                                                                                                                                   cGACTTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACGGCGGAGCAAGCAGCGTCGGGCGACg  >  2:1090968/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GGTCAACATAATACTGCGCCAGCAGACGAATTTGCTTCATTAGGGAATTCCTTCACGACAACCTGCCCCAATAATACCCAATTCTGACAGTTACGTTATAAATTGTGTGAGAAATGCGGAATTCCCCTCAACGACCTACTCCGTTGGAGGTAGCCACGGACGTCCCAGCGCCGATCCCTGTACGCCATGTTCATTCAGATAGCGATCCAGCTCGACCATTCCCGTCCATCGATTTTCGCACCACAGCGGAGCAAGCAGCGTCGGGCGACG  >  NC_000913/3352986‑3353255

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