Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,927,062 +G coding (2/1284 nt)
coding (169/183 nt)
yieM ←
b4564 ←
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,927,0621.G99.1% 304.2 / NA 115M1A (GTG→GCG) †
C57R (TGC→CGC) 
yieM
b4564
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (59/55);  minor base T (0/1);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

AGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCACGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGA  >  NC_000913/3926987‑3927130
                                                                            |                                                                    
agatGGCTTTGTTGATAGCACATCACCTCTTCGGTGAGTTCGGGCGGGACACGGGCATCTTTCAGCCGACTGCGCa                                                                       <  2:2734512/73‑1 (MQ=255)
agCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCa                                                                       <  1:1210314/76‑1 (MQ=255)
agCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCa                                                                       <  2:4418158/76‑1 (MQ=255)
  cTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagc                                                                     <  2:3706612/76‑3 (MQ=255)
   tGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagcg                                                                    >  1:2939087/1‑73 (MQ=255)
   tGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagcg                                                                    >  1:2151408/1‑73 (MQ=255)
   tGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCagcg                                                                    <  2:1858738/76‑4 (MQ=255)
       tttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTc                                                                <  1:784653/76‑1 (MQ=255)
       tttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTc                                                                >  1:4613339/1‑76 (MQ=255)
       tttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTc                                                                >  1:4066584/1‑76 (MQ=255)
        ttGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCa                                                               <  2:2745732/76‑1 (MQ=255)
         tGCTGATAGCACATCGCCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCAc                                                              <  2:4157915/76‑1 (MQ=255)
         tGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTccc                                                              >  2:4203467/1‑74 (MQ=255)
         tGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCAc                                                              >  2:1104055/1‑76 (MQ=255)
         tGCTGATAGCACATCACCTCTTCGGTGAGCTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCAc                                                              <  1:4334984/76‑1 (MQ=255)
           cTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGGCTGCGCAGCGCTTCACGc                                                            <  2:3061243/76‑1 (MQ=255)
           cTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCCCGc                                                            >  2:156836/1‑76 (MQ=255)
           cTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGc                                                            <  1:4698891/76‑1 (MQ=255)
            tGATAGCACTTCACCTCTTGGGTGGGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGcc                                                           <  2:761277/76‑1 (MQ=255)
            tGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGcc                                                           <  1:1709363/76‑1 (MQ=255)
              aTAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAg                                                         <  2:526420/76‑1 (MQ=255)
              aTAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAg                                                         <  1:2459159/76‑1 (MQ=255)
               tAGCGCATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                        <  2:4181636/76‑1 (MQ=255)
               tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                        >  1:2690752/1‑76 (MQ=255)
               tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                        <  2:4139955/76‑1 (MQ=255)
               tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                        <  2:2866702/76‑1 (MQ=255)
               tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                        <  1:4801220/76‑1 (MQ=255)
               tAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGc                                                        <  1:3561872/76‑1 (MQ=255)
                aGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCg                                                       >  1:4089567/1‑76 (MQ=255)
                 gcgcaTCACCTCTTCGGTGAGTTCCGGCGGGACTAGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgg                                                      <  2:2350375/73‑1 (MQ=255)
                 gCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgg                                                      <  1:2797115/76‑1 (MQ=255)
                   acaTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgggg                                                    <  1:668902/76‑1 (MQ=255)
                   acaTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgggg                                                    >  1:1951967/1‑76 (MQ=255)
                   acaTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgggg                                                    >  1:2583475/1‑76 (MQ=255)
                   acaTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACCACAGCgcgg                                                    >  2:3522550/1‑76 (MQ=255)
                   acaTCACCTCGTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCgggg                                                    <  2:1533656/76‑1 (MQ=255)
                     aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGaa                                                  >  1:3406825/1‑76 (MQ=255)
                     aTCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGaa                                                  >  2:514430/1‑76 (MQ=255)
                      tCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAc                                                 <  1:2477597/76‑1 (MQ=255)
                      tCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAc                                                 <  1:2512154/76‑1 (MQ=255)
                       cACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGACca                                                >  2:65687/1‑76 (MQ=255)
                       cACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAca                                                >  2:1209976/1‑76 (MQ=255)
                       cACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAca                                                <  1:4562180/76‑1 (MQ=255)
                         ccTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcatc                                              >  1:4862966/1‑76 (MQ=255)
                          ctctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcatca                                             >  1:1642209/1‑76 (MQ=255)
                          ctctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcatca                                             >  1:1456951/1‑76 (MQ=255)
                           tctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAAcatcat                                            >  1:4020329/1‑76 (MQ=255)
                            ctTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACatcatc                                           >  2:449989/1‑76 (MQ=255)
                             ttCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACAtcatca                                          >  1:374507/1‑76 (MQ=255)
                             tgcGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACAtcatca                                          <  2:2584653/74‑1 (MQ=255)
                                ggTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTg                                       >  2:3162056/1‑76 (MQ=255)
                                ggTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTg                                       >  2:3198417/1‑76 (MQ=255)
                                ggTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGCAACATCATCACTg                                       >  2:3015946/1‑76 (MQ=255)
                                  tGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGAt                                     >  2:2012207/1‑76 (MQ=255)
                                   gAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGAtt                                    >  2:2486262/1‑76 (MQ=255)
                                    aGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTg                                   <  1:4270425/76‑1 (MQ=255)
                                    aGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTg                                   >  1:1426542/1‑76 (MQ=255)
                                     gTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgc                                  <  2:1029428/76‑1 (MQ=255)
                                     gTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgc                                  >  2:886604/1‑76 (MQ=255)
                                      ttCCGGCGGGCCTCGGGCATCTTTCAGCCGACTGCGCCGCGCTCCACGCCAGCGGGGACCATCACCAGTGATTgcc                                 >  1:4464056/1‑76 (MQ=255)
                                      ttCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgcc                                 <  1:834271/76‑1 (MQ=255)
                                       tCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccg                                >  2:2084408/1‑76 (MQ=255)
                                       tCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccc                                >  2:981142/1‑75 (MQ=255)
                                        ccGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgc                               >  2:1872803/1‑76 (MQ=255)
                                        ccGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgc                               <  1:791883/76‑1 (MQ=255)
                                        ccGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgc                               >  2:4293178/1‑76 (MQ=255)
                                        ccGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgc                               <  1:4639504/76‑1 (MQ=255)
                                         cggcggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgcc                              >  1:2898145/1‑76 (MQ=255)
                                         cggcggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTgccgcc                              >  1:2124612/1‑76 (MQ=255)
                                          ggcggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCt                             >  1:2150809/1‑76 (MQ=255)
                                          ggcggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCt                             >  1:3695798/1‑76 (MQ=255)
                                          ggcggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGCGATTGCCGCCt                             >  2:757614/1‑76 (MQ=255)
                                           gcggGCCGCGGGCATCTTTCAGCCGACTGCGCAGCGCTGCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCtt                            <  1:3803932/76‑1 (MQ=255)
                                            tgggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTc                           >  2:3729047/2‑76 (MQ=255)
                                            cggGACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTc                           <  1:2320735/76‑1 (MQ=255)
                                            agggAATCGGGCATCTTTCACCAGAAAGCCCAGCGCTTCACGACAGCAGCGAACATCATCAGTAATTGACGCCTTc                           >  1:581554/2‑76 (MQ=255)
                                             gggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCa                          <  1:2441597/76‑1 (MQ=255)
                                             gggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCa                          <  2:3331346/76‑1 (MQ=255)
                                              ggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAg                         <  1:2796115/76‑1 (MQ=255)
                                              ggACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAg                         <  1:1454903/76‑1 (MQ=255)
                                               gccTCGGGCATCTTTCAGCCGCCTGCGCAGCGCTTCACGCCCGCGGGGAACATCATCAGTGCTTGCCGCCTTCAgg                        >  1:2608757/3‑76 (MQ=255)
                                               gACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAgg                        >  1:3274327/1‑76 (MQ=255)
                                               gACTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAgg                        >  1:2576033/1‑76 (MQ=255)
                                                aCTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGc                       <  2:2838859/76‑1 (MQ=255)
                                                aCTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGc                       <  1:2182206/76‑1 (MQ=255)
                                                 cTCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCg                      >  2:2615698/1‑76 (MQ=255)
                                                  tCGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGt                     >  2:2510261/1‑76 (MQ=255)
                                                   cGGGCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTg                    <  2:489153/76‑1 (MQ=255)
                                                    gggCATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTgg                   >  1:1611138/1‑76 (MQ=255)
                                                       cATCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGaa                >  2:412105/1‑76 (MQ=255)
                                                        aTCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCGGGCGTGGGAAt               <  2:4200910/76‑1 (MQ=255)
                                                         tCTTTCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAAtt              >  1:322470/1‑76 (MQ=255)
                                                             tCAGCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTc          <  1:96613/76‑1 (MQ=255)
                                                                gCCGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCaaa       >  1:354965/1‑76 (MQ=255)
                                                                 ccGCCTGCGCAGCGCTTCGCGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaag      <  1:4574279/76‑1 (MQ=255)
                                                                 ccGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaag      <  1:4573482/76‑1 (MQ=255)
                                                                 ccGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaag      <  2:2163023/76‑1 (MQ=255)
                                                                  cGACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaaga     >  2:373555/1‑76 (MQ=255)
                                                                   gACTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaagaa    <  1:1292660/76‑1 (MQ=255)
                                                                    aCTGCGCAGCGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAaagaag   <  2:2251172/76‑1 (MQ=255)
                                                                            |                                                                    
AGCTGGCTTTGCTGATAGCACATCACCTCTTCGGTGAGTTCCGGCGGGACTCGGGCATCTTTCAGCCGACTGCGCACGCTTCACGCCAGCGGGGAACATCATCAGTGATTGCCGCCTTCAGGCGTGGGAATTTTTCAAAGAAGA  >  NC_000913/3926987‑3927130

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