Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,171,255 T→A 15.8% intergenic (+128/+55) PP_4548 → / ← fadD‑I oxidoreductase/long‑chain‑fatty acid‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,171,2550TA15.8% 137.7 / 18.0 57intergenic (+128/+55)PP_4548/fadD‑Ioxidoreductase/long‑chain‑fatty acid‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (21/27);  new base A (3/6);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TTGAGAAGCTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCGGCAGCGAA  >  NC_002947/5171111‑5171400
                                                                                                                                                |                                                                                                                                                 
ttGAGAAGCTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCgg                                                                                                                                              <  1:943011/150‑1 (MQ=255)
    gaAGCTGCGCTGAGGTGCTCCTCCCGAGATTCTCAATTCCTGAAGACGGCGCCGCATCTTATTTTGGGGCCTGGTGTGGCGTAAGTGTCGCAATGGGCCCCCAGGCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCt                                                                                                                                          <  2:222863/150‑1 (MQ=255)
     aaGCTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtg                                                                                                                                         <  2:1124948/150‑2 (MQ=255)
      aGCTGCGCTGAGTTGCGCCTACAGGCACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtgg                                                                                                                                        <  1:518527/150‑3 (MQ=255)
        cTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTcg                                                                                                                                      >  2:672808/1‑150 (MQ=255)
        cTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTcg                                                                                                                                      >  1:121854/1‑150 (MQ=255)
           cgcTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtg                                                                                                                                         <  2:501207/144‑2 (MQ=255)
           cgcTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtg                                                                                                                                         >  1:501207/1‑143 (MQ=255)
                        cTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtgg                                                                                                                                        >  2:280027/1‑130 (MQ=255)
                        cTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtgg                                                                                                                                        <  1:280027/132‑3 (MQ=255)
                           cAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCca                                                                                                                   <  1:960047/150‑1 (MQ=255)
                              ggACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCacac                                                                                                                >  1:455326/1‑150 (MQ=255)
                               gACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACt                                                                                                               <  1:326824/150‑1 (MQ=255)
                                  tCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGAGCTTTGTGTGGCGAAATGGGCGCATAACGGCCCCACCCCTTGATGCATCACCGCCGAAATTCCCGGGCGGCATCGCGGCCCTTTCGCGACACGAGTCCGCCCACACCCCGAg                                                                                                            >  2:576715/1‑150 (MQ=255)
                                    acaATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCca                                                                                                                   >  2:548991/1‑141 (MQ=255)
                                    acaATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCca                                                                                                                   <  1:548991/141‑1 (MQ=255)
                                      aaTGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTg                                                                                                        >  2:241190/1‑150 (MQ=255)
                                          tCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAgg                                                                                                    >  2:974629/1‑150 (MQ=255)
                                          tCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAgg                                                                                                    <  1:1303768/150‑1 (MQ=255)
                                          tCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtgg                                                                                                                                        <  2:48171/114‑3 (MQ=255)
                                          tCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGACGCGGCCCtgg                                                                                                                                        >  1:48171/1‑112 (MQ=255)
                                          tCTGAAGACGGCCCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAgg                                                                                                    <  1:495417/150‑1 (MQ=255)
                                            tGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCt                                                                                                  >  2:464700/1‑150 (MQ=255)
                                            tGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCt                                                                                                  >  2:1049373/1‑150 (MQ=255)
                                            tGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCt                                                                                                  >  1:499470/1‑150 (MQ=255)
                                             gAAGCCGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCAAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCtt                                                                                                 <  1:784531/150‑1 (MQ=255)
                                             gAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCtt                                                                                                 <  1:1049373/150‑1 (MQ=255)
                                                     cgccgcATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAg                                                                                         <  2:232044/150‑1 (MQ=255)
                                                            tCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGAtct                                                                                  >  1:145691/1‑150 (MQ=255)
                                                              ttATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGAtcttc                                                                                <  1:591581/150‑1 (MQ=255)
                                                                       ggCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCca                                                                                                                   >  1:345725/1‑106 (MQ=255)
                                                                       ggCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCca                                                                                                                   <  2:345725/106‑1 (MQ=255)
                                                                        gCCTTGTGTCGCGAAAGGGGCGCAAAGCGCGCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTtgtg                                                                                              >  1:401851/1‑126 (MQ=255)
                                                                        gCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTtgtg                                                                                              <  2:401851/126‑1 (MQ=255)
                                                                         ccTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTg                                                                     <  2:455326/150‑1 (MQ=255)
                                                                              tgtCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATcttctt                                                                              >  1:1387629/1‑136 (MQ=255)
                                                                              tgtCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATcttctt                                                                              <  2:1387629/136‑1 (MQ=255)
                                                                               gtCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGcttctt                                                               <  1:658446/150‑1 (MQ=255)
                                                                               gtCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGcttctt                                                               <  1:241190/150‑1 (MQ=255)
                                                                                  gcgAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGAGTTTTACTGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGAGACCTTACTAAAGACTTTCTTCTGGAg                                                            >  1:88286/1‑150 (MQ=255)
                                                                                      tgggggCGCTTTGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCTCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGGAGGCTTGTGGTCAGGCGATCTTCTTCAAGACTTGCTTCTTTAGctct                                                        >  1:136285/3‑150 (MQ=255)
                                                                                           gcgcAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGc                                                                                                           <  2:972473/94‑1 (MQ=255)
                                                                                           gcgcAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGc                                                                                                           >  1:972473/1‑94 (MQ=255)
                                                                                               tttgcccccccAGCCATTGATGCATCGCCGCCGAAATGCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCTGCCCCCACACTGAGCGTGCAGCCTTGTGGTCATGCCATCTTCTTTAAGCCTTGCTTCTTGTGCTCGTCATCACGc                                               >  2:552220/7‑150 (MQ=255)
                                                                                                aaGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCa                                              <  1:974629/150‑1 (MQ=255)
                                                                                                           gCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGgcgc                                   <  1:156269/150‑1 (MQ=255)
                                                                                                                   tGCATCGCCGCCGAAATCCCACGGCCTCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTATTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCtt                           <  2:1360298/150‑1 (MQ=255)
                                                                                                                     cATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGc                         <  1:979676/150‑1 (MQ=255)
                                                                                                                      aTCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGcc                        <  2:145691/150‑1 (MQ=255)
                                                                                                                      aTCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGcc                        >  2:711674/1‑150 (MQ=255)
                                                                                                                      aTCGCCGCCGAAATCACAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGcc                        >  1:240251/1‑150 (MQ=255)
                                                                                                                               gAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTtggt               <  1:1010179/150‑1 (MQ=255)
                                                                                                                                     ccAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGa                                                             >  1:147707/1‑98 (MQ=255)
                                                                                                                                     ccAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGa                                                             <  2:147707/98‑1 (MQ=255)
                                                                                                                                     ccAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCgg         <  2:1052322/150‑1 (MQ=255)
                                                                                                                                        gggCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCTACACTGAGCGTGCAGGCTTGTGGTCAGCCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCGGCAg      <  1:576715/150‑1 (MQ=255)
                                                                                                                                            cgcgTCGCGGCCCTTTCGCGACACGAGGCAGACCACCCTCTTGCCGGGCAGGCTTGTGGTCAGGCGATCCTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCGGCAGCGaa  <  1:672808/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TTGAGAAGCTGCGCTGAGTTGCGCCTACAGGGACTCACAATGTCTGAAGACGGCGCCGCATCTTATTGTGGGGCCTTGTGTCGCGAAAGGGGCGCAAAGCGCCCCCAGCCATTGATGCATCGCCGCCGAAATCCCAGGGCCGCGTCGCGGCCCTTTCGCGACACGAGGCCGCCCCCACACTGAGCGTGCAGGCTTGTGGTCAGGCGATCTTCTTCAAGCCTTGCTTCTTGAGCTCTTCATCACGCAGTTCGCGGCGCAGGATCTTGCCCACGTTGGTGGTCGGCAGCGAA  >  NC_002947/5171111‑5171400

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