Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,389,282 A→C 12.5% E48A (GAG→GCG) ‡ rlmKL → 23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,389,2820AC12.5% 146.6 / 12.9 56E48A (GAG→GCG) ‡rlmKL23S rRNA (guanine(2445)‑N(2))/(guanine(2069)‑N(7))‑ methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/21);  new base C (3/4);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

GGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCG  >  NC_002947/2389134‑2389427
                                                                                                                                                    |                                                                                                                                                 
ggggCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGa                                                                                                                                                   <  7:42931/149‑1 (MQ=255)
          gTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGCCGAGGTGTGCTCGCGCACCTCg                                                                                                                                         >  7:31827/1‑149 (MQ=255)
           tCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCgg                                                                                                                                        <  2:56629/149‑1 (MQ=255)
            cGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGCCGAGGTGTGCTCGCGCACCTCgtc                                                                                                                                       <  8:31827/149‑3 (MQ=255)
            cGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGc                                                                                                                                       >  7:63052/1‑149 (MQ=255)
             gCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGTGCTCGCGCACCTCgtc                                                                                                                                       <  5:47995/148‑3 (MQ=255)
             gCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGTGCTCGCGCACCTCgtc                                                                                                                                       <  7:320656/148‑3 (MQ=255)
              ccccGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCATGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCa                                                                                                                                     >  2:307395/1‑149 (MQ=255)
                       tcaAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCgcggcg                                                                                                                            >  8:267504/1‑149 (MQ=255)
                         aaaTAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGgcgc                                                                                                                           <  4:67019/148‑1 (MQ=255)
                         aaaTAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGcgcc                                                                                                                          >  8:45011/1‑149 (MQ=255)
                         aaaTAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGcgcc                                                                                                                          >  4:79464/1‑149 (MQ=255)
                           aTAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGcgccgc                                                                                                                        >  2:140331/1‑149 (MQ=255)
                               cccAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGa                                                                                                                     >  6:116335/1‑148 (MQ=255)
                                 cAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACAt                                                                                                                  <  8:63052/149‑1 (MQ=255)
                                      cgcTGGAATAGTCATGTCGTACCGTTTCGAACTTTACCTCCCCTGTCCCAAAGGGCTTGACAGCCTGGTCGCCGATGAGGGCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCGATTCGCGGCGCCGCCGACATGGaa                                                                                                              <  8:266426/148‑1 (MQ=255)
                                             aTACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTa                                                                                                      >  2:39415/1‑149 (MQ=255)
                                                 tCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGc                                                                                                  >  1:51487/1‑149 (MQ=255)
                                                        ggACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCtgtgt                                                                                            <  1:126214/148‑1 (MQ=255)
                                                         gACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTgtgtgt                                                                                          <  7:267504/149‑1 (MQ=255)
                                                              tttCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGt                                                                                     >  1:50979/1‑149 (MQ=255)
                                                                   aaCTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGt                                                                                <  1:39415/149‑1 (MQ=255)
                                                                        tACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTgg                                                                            <  5:248518/148‑1 (MQ=255)
                                                                        tACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTgg                                                                            <  3:79464/148‑1 (MQ=255)
                                                                          ccTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGcc                                                                          <  1:140331/148‑1 (MQ=255)
                                                                           cTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACt                                                                                                           <  1:69326/114‑1 (MQ=255)
                                                                           cTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACt                                                                                                           >  2:69326/1‑114 (MQ=255)
                                                                           cTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCaa                                                                        >  1:126834/1‑149 (MQ=255)
                                                                               ccTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACcgc                                                                    <  7:45011/149‑1 (MQ=255)
                                                                               ccTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACcgc                                                                    <  2:51487/149‑1 (MQ=255)
                                                                                     ccAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCt                                                                                                        <  4:114936/107‑1 (MQ=255)
                                                                                     ccAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCt                                                                                                        <  8:137203/107‑1 (MQ=255)
                                                                                     ccAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCt                                                                                                        >  3:114936/1‑107 (MQ=255)
                                                                                     ccAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCt                                                                                                        >  7:137203/1‑107 (MQ=255)
                                                                                      cAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACGTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGc                                                             <  5:188136/149‑1 (MQ=255)
                                                                                          ggCCTTGACAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGTGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGt                                                         <  8:347073/149‑1 (MQ=255)
                                                                                                gAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTc                                                     >  7:213955/1‑147 (MQ=255)
                                                                                                  aaGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAgc                                                 >  7:71918/1‑149 (MQ=255)
                                                                                                       tGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTtc                                            >  1:373863/1‑149 (MQ=255)
                                                                                                          tCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTtctc                                          >  4:8625/1‑148 (MQ=255)
                                                                                                            gCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGAt                                       >  3:279766/1‑149 (MQ=255)
                                                                                                            gCCGAAGAGGCCAAGGGCCTCGGCCCTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGTCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGAt                                       >  2:230512/1‑149 (MQ=255)
                                                                                                               gAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATgaa                                    <  2:126834/149‑1 (MQ=255)
                                                                                                                    ggCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACg                               >  2:258358/1‑149 (MQ=255)
                                                                                                                          gggCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGAGCTTCTCGATGAAGAACGCcgac                          >  2:251266/1‑148 (MQ=255)
                                                                                                                              cTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCGCGATGAAGATCGCCGACGACCt                     >  6:310550/1‑149 (MQ=255)
                                                                                                                                   ccTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTac                 >  3:331800/1‑148 (MQ=255)
                                                                                                                                    cTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTacga               >  1:233408/1‑149 (MQ=255)
                                                                                                                                       gACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACgg            >  2:189184/1‑149 (MQ=255)
                                                                                                                                        aCGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGgt           <  2:50979/149‑1 (MQ=255)
                                                                                                                                         cGAGGTGTGCTCGCGCACCTCGTCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCtgtgt                                                                                            >  3:334667/1‑67 (MQ=25)
                                                                                                                                         cGAGGTGTGCTCGCGCACCTCGTCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCtgtgt                                                                                            <  4:334667/67‑1 (MQ=25)
                                                                                                                                         cGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGgtg          >  4:277777/1‑149 (MQ=255)
                                                                                                                                          gAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGgtgt         <  4:279766/149‑1 (MQ=255)
                                                                                                                                           aGGTGTGCTCGCGCACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGtgtg        >  5:104296/1‑149 (MQ=255)
                                                                                                                                                cgcgAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGTCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCAcgc   <  1:230512/149‑1 (MQ=255)
                                                                                                                                                 gcgAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCAcgcg  >  3:62021/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGGGCGATGCGTCGCCCCGTCAGTCAAATAGCCCAAGGCGCTGGAATACTCATGTCGGACCGTTTCGAACTTTACCTCACCTGTCCCAAAGGCCTTGAAAGCCTGCTCGCCGAAGAGGCCAAGGGCCTCGGCCTTGACGAGGTGCGCGAGCACACCTCGGCCATTCGCGGCGCCGCCGACATGGAAACTGCTTACCGCCTGTGTGTCTGGTCACGTCTGGCCAACCGCGTATTGCTGGTGCTCAAGCGCTTCTCGATGAAGAACGCCGACGACCTTTACGACGGTGTGCACGCG  >  NC_002947/2389134‑2389427

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