Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,563,403 T→A 7.5% K141* (AAG→TAG)  PP_3146 ← oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,563,4030TA 7.5% 151.3 / 2.6 53K141* (AAG→TAG) PP_3146oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base A (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

TCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGCACACCGTCCTGGCC  >  NC_002947/3563258‑3563550
                                                                                                                                                 |                                                                                                                                                   
tCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGaa                                                                                                                                                  <  6:101537/149‑1 (MQ=255)
  cAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCGCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAg                                                                                                                                                 <  1:279566/148‑1 (MQ=255)
   aGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAggc                                                                                                                                               <  3:70416/149‑1 (MQ=255)
              gTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTgcg                                                                                                                                    >  4:19620/1‑149 (MQ=255)
                ggCCCGAGCCCATGTCGACCAGGCCACCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGcgcc                                                                                                                                  >  5:88015/1‑149 (MQ=255)
                         ccATGTCGACCAGGCCACCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTg                                                                                                                         <  6:88015/149‑1 (MQ=255)
                              gcGACCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCa                                                                                                                     <  5:204588/147‑1 (MQ=255)
                                 aCCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGgcggcg                                                                                                                                           >  5:115048/1‑123 (MQ=255)
                                 aCCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGgcggcg                                                                                                                                           <  6:115048/123‑1 (MQ=255)
                                 aCCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTc                                                                                                                 <  2:108063/149‑1 (MQ=255)
                                            aGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGccc                                                                                                       >  2:327371/1‑148 (MQ=255)
                                             gCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGccc                                                                                                       >  1:304875/1‑147 (MQ=255)
                                              cGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAg                                                                                                    <  6:201933/149‑1 (MQ=255)
                                                       cGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCg                                                                                           >  1:289669/1‑149 (MQ=255)
                                                         cTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCa                                                                                         <  3:19620/149‑1 (MQ=255)
                                                              gacgaTTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGaa                                                                                                                                                  <  8:191267/87‑1 (MQ=255)
                                                                 gaTTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGc                                                                                 <  7:213200/149‑1 (MQ=255)
                                                                  aTTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGcc                                                                                <  2:219142/149‑1 (MQ=255)
                                                                       gTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGc                                                                           >  2:172644/1‑149 (MQ=255)
                                                                        ttCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGCACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAAGGCCGTAGCg                                                                          >  6:199598/1‑149 (MQ=255)
                                                                         tCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTAGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCgg                                                                         >  1:277848/1‑149 (MQ=255)
                                                                         tCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTAGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCgg                                                                         >  1:274296/1‑149 (MQ=255)
                                                                           aaCTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCggcc                                                                       >  8:292228/1‑147 (MQ=255)
                                                                           aaCTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGc                                                                       >  3:300868/1‑149 (MQ=255)
                                                                           aaCTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGc                                                                       >  6:57324/1‑149 (MQ=255)
                                                                           aaCTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGc                                                                       >  6:84641/1‑149 (MQ=255)
                                                                                cccGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGAt                                                                    <  1:327371/147‑1 (MQ=255)
                                                                                   ggCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCt                                                               <  4:139385/149‑1 (MQ=255)
                                                                                          aCCAGGCGTGTTTCATGCCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTCCGCCGCGGTGTTGGCCACTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAACGGGCGATGCGCTCGGacac                                                        >  6:23927/1‑149 (MQ=255)
                                                                                                  tgtTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCa                                                <  5:57324/149‑1 (MQ=255)
                                                                                                       cATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGa                                            <  5:84641/148‑1 (MQ=255)
                                                                                                         tGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTc                                         <  1:172644/149‑1 (MQ=255)
                                                                                                              aTAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGc                                    <  3:1166/149‑1 (MQ=255)
                                                                                                                 gCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCg                                  <  2:304875/148‑1 (MQ=255)
                                                                                                                  cccAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGTCGAACTTCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATACCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGc                                <  2:289669/149‑1 (MQ=255)
                                                                                                                  cccAGGGCGGCCAGGTCGTCCTGTTCCTCCTAGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGc                                <  2:274296/149‑1 (MQ=255)
                                                                                                                  cccAGGGCGGCCAGGTCGTCCTGTTCCTCCTAGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGc                                <  2:277848/149‑1 (MQ=255)
                                                                                                                   ccAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGcc                               >  2:302654/1‑149 (MQ=255)
                                                                                                                   ccAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCCCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGcc                               >  6:148283/1‑149 (MQ=255)
                                                                                                                      gggCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGc                            <  7:292228/149‑1 (MQ=255)
                                                                                                                        gcggcCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACt                                              <  7:254807/129‑1 (MQ=255)
                                                                                                                        gcggcCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACt                                              >  8:254807/1‑129 (MQ=255)
                                                                                                                          ggcCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCtt                          >  5:285183/1‑147 (MQ=255)
                                                                                                                              aGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAgg                      >  1:204984/1‑147 (MQ=255)
                                                                                                                                gtcgtcCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAgc                  >  3:243100/1‑149 (MQ=255)
                                                                                                                                    tcCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGca               >  1:153188/1‑148 (MQ=255)
                                                                                                                                        gTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGCACACCg          >  6:119557/1‑149 (MQ=255)
                                                                                                                                            ctcctTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGCACACCGTCCt      >  5:282581/1‑149 (MQ=255)
                                                                                                                                             tcctTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCt                              <  4:239205/124‑1 (MQ=255)
                                                                                                                                             tcctTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCt                              >  3:239205/1‑124 (MQ=255)
                                                                                                                                             tcctTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGCACACCGTCCTg     >  1:299284/1‑149 (MQ=255)
                                                                                                                                              cctTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGCACACCGTCCTgg    >  4:130595/1‑149 (MQ=255)
                                                                                                                                               ctTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGCACACCGTCCTg     <  6:285183/147‑1 (MQ=255)
                                                                                                                                                ttGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGTCGGCCTGCTTCAGGTAGCGCACACCGTCCTGGcc  >  6:248691/1‑149 (MQ=255)
                                                                                                                                                 tGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGCACACCGTCCTGGcc  >  6:307816/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCCAGCGGGTGCAGGTGGCCCGAGCCCATGTCGACCAGGCCGCCAGCGTACTGGTCGCTGGCGACGATTTCGTTCAACTGCCCGGCGCTGACCAGGCGTGTTTCATGGCGATAGCCCAGGGCGGCCAGGTCGTCCTGTTCCTCCTTGAAGGCGGCGAACTGCGCCGCGGTGTTGGCCAGTTCGCAGAAGCCCCAGCGCAGATCGCAGGCAATGCCGTAGCGGGCGATGCGCTCGGCCACCAGTTCCACTGAGTCGATGCCGGCCTGCTTCAGGTAGCGCACACCGTCCTGGCC  >  NC_002947/3563258‑3563550

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: