Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 211,621 G→C 6.8% V5709V (GTG→GTC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947211,6210GC 6.8% 175.8 / 4.5 59V5709V (GTG→GTCPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base C (2/2);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

CGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTAC  >  NC_002947/211474‑211761
                                                                                                                                                    |                                                                                                                                            
cgCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGa                                                                                                                                             <  7:42697/149‑1 (MQ=255)
 gCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTCAc                                                                                                                                            <  2:44984/149‑1 (MQ=255)
 gCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACCGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGAc                                                                                                                                            <  5:184687/149‑1 (MQ=255)
  cTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGAcc                                                                                                                                           <  6:102090/149‑1 (MQ=255)
   tGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGAcc                                                                                                                                           >  8:96528/1‑148 (MQ=255)
    ggCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCg                                                                                                                                          >  8:184918/1‑148 (MQ=255)
        aGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAggt                                                                                                                                     >  8:102578/1‑149 (MQ=255)
        aGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAggt                                                                                                                                     >  6:329834/1‑149 (MQ=255)
           tCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAggtggt                                                                                                                                  <  1:12313/149‑1 (MQ=255)
           tCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAggtggt                                                                                                                                  <  7:146632/149‑1 (MQ=255)
             aGCCTGGGCATCAACTCGGCGGTGGACGCCACTGTTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGATCGCGACGCCTTCGGTGACCGAAGgtggtg                                                                                                                                 <  8:296451/148‑1 (MQ=255)
             aGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTga                                                                                                                                <  1:38602/149‑1 (MQ=255)
             aGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTCACCGAAggcgtcg                                                                                                                                 <  8:171803/148‑6 (MQ=255)
               ccTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAggtggt                                                                                                                                  >  5:237866/1‑145 (MQ=255)
                      aTCAACTCGGAGGTGGAAGCCACTGGTGAGAACTTCGAGAACCTCGAGCTTGGCGGTGCGGATTCGGTTCAAGTGAACGATACAACCGACGAAGTGTTGGCCAAGCTGACCGCGACGCCTTCGGTGAACGAAGGTGGTGAGATAACCTa                                                                                                                       >  6:153172/1‑149 (MQ=255)
                         aaCTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTacac                                                                                                                    >  4:159326/1‑149 (MQ=255)
                                         ccACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACcaac                                                                                                    <  6:344670/149‑1 (MQ=255)
                                               gggCGACCTTCGCTAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCAACACCATCACCCTGACCAACAAAg                                                                                                <  6:237866/145‑1 (MQ=255)
                                                  cGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGt                                                                                           <  1:171629/149‑1 (MQ=255)
                                                   gACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAGAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGCCCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATgg                                                                                            >  1:65255/1‑148 (MQ=255)
                                                     ccTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATg                                                                                             <  6:281151/144‑1 (MQ=255)
                                                     ccTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATg                                                                                             >  5:281151/1‑144 (MQ=255)
                                                      cTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGc                                                                                       >  7:316076/1‑149 (MQ=255)
                                                        tCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGcc                                                                                      >  5:106877/1‑148 (MQ=255)
                                                            gaACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATc                                                                                  >  4:273109/1‑148 (MQ=255)
                                                               ccTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaaca                                                                              <  1:177011/149‑1 (MQ=255)
                                                                cTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATcaacaa                                                                             >  5:198723/1‑149 (MQ=255)
                                                                     gCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACcaca                                                                        <  2:251815/149‑1 (MQ=255)
                                                                        gggCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTg                                                                     >  2:18913/1‑149 (MQ=255)
                                                                         ggCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGt                                                                      <  6:272229/147‑1 (MQ=255)
                                                                           cGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGATCTGCCGATCAACAACCACAGTGCAt                                                                  <  5:329834/149‑1 (MQ=255)
                                                                            ggTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATTACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCAtt                                                                 <  1:74314/149‑1 (MQ=255)
                                                                                cGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAcc                                                             <  3:159326/149‑1 (MQ=255)
                                                                                cGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                              <  1:18913/148‑1 (MQ=255)
                                                                                cGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                              <  5:260595/148‑1 (MQ=255)
                                                                                cGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGAc                                                              <  7:96528/148‑1 (MQ=255)
                                                                                  gCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCtt                                                           <  6:262942/149‑1 (MQ=255)
                                                                                    ttCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGAt                                                                                   >  1:290355/1‑123 (MQ=255)
                                                                                    ttCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGAt                                                                                   <  2:290355/123‑1 (MQ=255)
                                                                                        gTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCt                                                     <  2:65255/149‑1 (MQ=255)
                                                                                            aaGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACACTGAGc                                                 >  5:296604/1‑149 (MQ=255)
                                                                                                        accaccGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGa                                                   <  2:91184/136‑1 (MQ=255)
                                                                                                        accaccGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGa                                                   >  1:91184/1‑136 (MQ=255)
                                                                                                          caccGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCg                                    >  6:306518/1‑149 (MQ=255)
                                                                                                            ccGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGt                                   >  2:254496/1‑148 (MQ=255)
                                                                                                             cgacgaAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGtc                                  >  3:248879/1‑148 (MQ=255)
                                                                                                                   aGTGGTGGCCAAGCTGACAGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACAACATCACACTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACACTGAGCGACGGCAAAACAGTCATCACCg                           >  2:114513/1‑149 (MQ=14)
                                                                                                                        tggCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCATCATCACCGTGCCg                      >  2:264334/1‑149 (MQ=255)
                                                                                                                          gCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGc                    >  1:328154/1‑149 (MQ=255)
                                                                                                                              aGCTGACCGCGACGCCTTCGGTGACCGAATGTGGTTAGATCACCTACACCATCACCCTGAACAACAAAGATTGTCTGCCGATCAACAACCACAGAGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCAACGTGCCGGCTAAt                <  5:153172/149‑1 (MQ=14)
                                                                                                                              aGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAAt                <  8:316076/149‑1 (MQ=255)
                                                                                                                                   aCCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACTTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCAc           <  1:254496/149‑1 (MQ=255)
                                                                                                                                       cgACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTg       <  7:102578/149‑1 (MQ=255)
                                                                                                                                       caccaccTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTg       >  4:154827/6‑149 (MQ=255)
                                                                                                                                        gACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGttg      <  6:296604/149‑3 (MQ=255)
                                                                                                                                        gACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTgg      <  5:306518/149‑1 (MQ=255)
                                                                                                                                        accaccTTCGGTCACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTgg      >  8:104104/5‑149 (MQ=255)
                                                                                                                                        accaccTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAACATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTgg      >  5:240916/5‑149 (MQ=255)
                                                                                                                                            ccTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTAc  <  4:154345/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CGCTGGCCAGATCAGCCTGGGCATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCGAGAACCTCGAGCTGGGCGGTGCGGCTTCGGTTCAAGTGACCGATACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACGCCTTCGGTGACCGAAGGTGGTGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACCCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCTAATGGCACCGTGGGTAC  >  NC_002947/211474‑211761

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