Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,806,762 T→C intergenic (+617/‑876) proX → / → b2682 high‑affinity transport system for glycine betaine and proline/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,806,7620TC99.3% 349.5 / NA 137intergenic (+617/‑876)proX/b2682high‑affinity transport system for glycine betaine and proline/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (59/77);  minor base A (1/0);  total (60/77)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 145 total aligned reads displayed.

TCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGG  >  NC_000913/2806689‑2806834
                                                                         |                                                                        
tCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTctgctgctg                                                                        >  2:3567247/1‑76 (MQ=255)
 ccGCGGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTgg                                                                       <  2:797877/76‑1 (MQ=255)
  cGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTggg                                                                      <  1:2240925/76‑1 (MQ=255)
   gCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTgggg                                                                     >  1:4902700/1‑76 (MQ=255)
    cgggCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGa                                                                    <  2:1332023/74‑1 (MQ=255)
     tGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGAt                                                                   >  1:100816/1‑76 (MQ=255)
     tGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGGTGGCACTATTATGAGCGGTCTGCTGCTGGGGAt                                                                   >  2:4854771/1‑76 (MQ=255)
     tGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGGGGGTGGCACTATTATGAGCGGTCTGCTGCTGGGGAt                                                                   >  1:513949/1‑76 (MQ=255)
       gCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCt                                                                 <  1:3969140/76‑1 (MQ=255)
         aGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTg                                                               >  1:84600/1‑76 (MQ=255)
         aGCGAAGCTGGCTCCAACGCATAAACGCGGCAAAGCGGCTGGAACTATTACGAGCCGCCCGGTGCTGGCGATCTTg                                                               >  2:798570/1‑76 (MQ=255)
          gCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGc                                                              >  1:3504553/1‑76 (MQ=255)
          gCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGc                                                              <  2:648766/76‑1 (MQ=255)
           cGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCt                                                             <  1:5012349/76‑1 (MQ=255)
            gACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTg                                                            <  2:4008452/76‑1 (MQ=255)
            gACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTg                                                            >  2:2486601/1‑76 (MQ=255)
            gACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTg                                                            >  1:2907934/1‑76 (MQ=255)
            gACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTg                                                            <  2:753901/76‑1 (MQ=255)
            gACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATActgctg                                                            >  2:2455244/1‑76 (MQ=255)
            gACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGGTGGCACTATTATTAGCGGTCTGCTGCTGGGCAGCTTGCTg                                                            >  1:3327456/1‑76 (MQ=255)
              cGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGc                                                          >  1:5504656/1‑76 (MQ=255)
               gCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCa                                                         >  1:3323622/1‑76 (MQ=255)
               gCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCa                                                         <  2:2773143/76‑1 (MQ=255)
               gCTGGCGTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCa                                                         <  2:4789055/76‑1 (MQ=255)
                  ggCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACgg                                                      <  1:3192243/76‑1 (MQ=255)
                  ggCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACgg                                                      >  2:78567/1‑76 (MQ=255)
                   gCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGa                                                     <  1:2117921/76‑1 (MQ=255)
                    gttCACCGGAGAAACGCGGCAAAGTGGTTGGTACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGAc                                                    <  2:5647715/75‑1 (MQ=255)
                    cTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGAc                                                    <  1:3489266/76‑1 (MQ=255)
                    cTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGAc                                                    >  1:5100015/1‑76 (MQ=255)
                    cTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGAc                                                    <  2:4276318/76‑1 (MQ=255)
                     ttCAGCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACa                                                   <  2:916776/76‑1 (MQ=255)
                     ttCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACa                                                   <  2:4449868/76‑1 (MQ=255)
                      tCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAg                                                  <  1:3761187/76‑1 (MQ=255)
                         ccGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTg                                               >  2:4785652/1‑76 (MQ=255)
                         ccGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTg                                               >  1:5101585/1‑76 (MQ=255)
                         ccGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTg                                               >  1:5441185/1‑76 (MQ=255)
                          cGGATAAACGCGGCAAAGCGGTTGGCACTATTATGAGCGGTCTGCTGATGGGGATCACGCTGGCACGGACCGTCGc                                              >  2:3209712/1‑76 (MQ=255)
                           ggATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGcc                                             <  2:2985983/76‑1 (MQ=255)
                           ggATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGcc                                             >  2:3291728/1‑76 (MQ=255)
                           ggATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGcc                                             >  2:3462209/1‑76 (MQ=255)
                            gATAAACGCGTCGAAGTGGTTGGTGCGATGGTGAGCGGTCTGCTGCGGGGGATCTTTTTGGTACGGCCTGTTGCCg                                            <  1:1711800/76‑1 (MQ=255)
                            gATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCg                                            <  2:4084875/76‑1 (MQ=255)
                            gATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCg                                            <  1:1714778/76‑1 (MQ=255)
                               aaaCGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGAt                                         <  2:611058/76‑1 (MQ=255)
                               aaaCGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGAt                                         <  1:203670/76‑1 (MQ=255)
                                aaCGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGAtt                                        >  1:2450694/1‑76 (MQ=255)
                                 aCGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCCGATTg                                       >  2:2693334/1‑76 (MQ=255)
                                   gcgGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCt                                     >  1:2392470/1‑76 (MQ=255)
                                   gcgGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCt                                     <  1:1920364/76‑1 (MQ=255)
                                    cgGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTg                                    <  1:2427442/76‑1 (MQ=255)
                                     ggCAAAGTGGTTGGCCCTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTgg                                   <  1:5097503/76‑1 (MQ=255)
                                     ggCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTgg                                   >  2:1622734/1‑76 (MQ=255)
                                      gCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGc                                  <  1:2461869/76‑1 (MQ=255)
                                      gCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGc                                  <  1:3317041/76‑1 (MQ=255)
                                      gCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGc                                  <  1:4716409/76‑1 (MQ=255)
                                       cAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCg                                 <  2:1522812/76‑1 (MQ=255)
                                        aaaGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGa                                >  1:1376471/1‑76 (MQ=255)
                                         aaGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGaa                               >  1:1097736/1‑76 (MQ=255)
                                          aGTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAt                              <  2:3467283/76‑1 (MQ=255)
                                           gTGGTTGGGACTATGATGGGGGGTCTGCGGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAtc                             <  2:116698/76‑1 (MQ=255)
                                           gTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAtc                             >  2:5073898/1‑76 (MQ=255)
                                           gTGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAtc                             <  1:3198188/76‑1 (MQ=255)
                                            tGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACa                                                   <  1:5137002/53‑1 (MQ=255)
                                            tGGTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAtct                            <  1:853570/76‑1 (MQ=255)
                                             ggTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAAtctc                           <  1:2860004/76‑1 (MQ=255)
                                              gTTGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGCATCTcg                          >  1:5808161/1‑76 (MQ=255)
                                              gTTGGCACTATCATGAGCGGTGTTCTGCTGGGGATCTTGCTGGGGCGGACAGTTGGCGGATTGCTGGGGAATCTcg                          <  2:5203284/76‑1 (MQ=255)
                                               ttGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTcgg                         <  1:4627221/76‑1 (MQ=255)
                                                tGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTcggc                        <  1:4762969/76‑1 (MQ=255)
                                                tGGCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTcggc                        >  1:3527616/1‑76 (MQ=255)
                                                 ggCACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTcggcg                       <  2:2958761/76‑1 (MQ=255)
                                                   cACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCggcggc                     <  1:438940/76‑1 (MQ=255)
                                                   cACTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCggcggc                     <  1:4430204/76‑1 (MQ=255)
                                                    aCTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCggcggcg                    >  1:5667569/1‑75 (MQ=255)
                                                    aCTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCt                    >  1:73627/1‑76 (MQ=255)
                                                    aCTATTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCt                    <  1:2994983/76‑1 (MQ=255)
                                                     cTATTATGGGGGGTCGGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTg                   <  2:216823/76‑1 (MQ=255)
                                                     cTAGTATGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTg                   <  2:1992561/76‑1 (MQ=255)
                                                      tattatGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTgg                  <  2:3619297/76‑1 (MQ=255)
                                                      tattatGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTgg                  >  2:1678774/1‑76 (MQ=255)
                                                       attatGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGgc                 <  2:5512231/76‑1 (MQ=255)
                                                       attatGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGgc                 <  1:962301/76‑1 (MQ=255)
                                                       attatGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGgc                 <  1:3961296/76‑1 (MQ=255)
                                                        ttatGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGgcg                <  2:2003585/76‑1 (MQ=255)
                                                         tatGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGgcgc               <  1:3449990/76‑1 (MQ=255)
                                                           tGAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCAc             >  1:2341341/1‑76 (MQ=255)
                                                           tGAGCGGTCTGCTGCTGGGGAACTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCAc             <  2:2518689/76‑1 (MQ=255)
                                                            gAGCGGTCTGCTGCTGTGCATATTGCTGGCGCGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCAcc            <  1:5468287/76‑1 (MQ=255)
                                                            gAGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCAcc            >  1:4469604/1‑76 (MQ=255)
                                                             aGCGGTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCCGGCGCaacc           >  2:3046824/1‑73 (MQ=255)
                                                              tgggTATGCTGCTGGAGAGCTTTCTGGCGCGCGCAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCGCCGt          <  1:500028/74‑1 (MQ=255)
                                                                ggTCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCt        <  1:3801807/76‑1 (MQ=255)
                                                                  tCTGCTGCTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCttt      <  1:673403/76‑1 (MQ=255)
                                                                   ctgctgctgGGGATATTGCGGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCtctt     >  2:4620708/1‑76 (MQ=255)
                                                                    tgctgcTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTg    >  1:5317477/1‑76 (MQ=255)
                                                                    tgctgcTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTg    >  1:2828523/1‑76 (MQ=255)
                                                                    tgctgcTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTg    >  2:911705/1‑76 (MQ=255)
                                                                     gctgctGGGGCTCTTGCTGGCGCGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTgg   <  1:5733578/76‑1 (MQ=255)
                                                                      ctgctgGGGATCTTGCTGGCTCGGACAGTTGCCGGATTGCTGGCGCATCTCGGCGGCTGGCGCACCGTCTTTTggg  <  2:1852456/76‑1 (MQ=255)
                                                                         |                                                                        
TCCGCTGGCAGCGACGCTGGCTTCACCGGATAAACGCGGCAAAGTGGTTGGCACTATTATGAGCGGTCTGCTGTTGGGGATCTTGCTGGCACGGACAGTTGCCGGATTGCTGGCGAATCTCGGCGGCTGGCGCACCGTCTTTTGGG  >  NC_000913/2806689‑2806834

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