Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,647 G→T 13.9% intergenic (+59/‑237) PP_2661 → / → PP_2662 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,048,6470GT13.9% 137.6 / 10.6 57intergenic (+59/‑237)PP_2661/PP_2662hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base T (3/5);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACT  >  NC_002947/3048504‑3048796
                                                                                                                                               |                                                                                                                                                     
ggATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCcg                                                                                                                                                 >  2:354730/1‑150 (MQ=255)
 gATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCcgc                                                                                                                                                <  2:828915/150‑1 (MQ=255)
  aTGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCcgcg                                                                                                                                               <  1:371839/150‑1 (MQ=255)
    gtgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCGGCCCTTCGCGGGCTCGCCCGCGaa                                                                                                                                             <  1:1106078/150‑1 (MQ=255)
       cAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGCCCGCGAAggg                                                                                                                                          <  2:478941/150‑1 (MQ=255)
         ggCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGcc                                                                                                                                        <  2:665249/150‑1 (MQ=255)
         ggCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAAATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGcc                                                                                                                                        <  1:28717/150‑1 (MQ=255)
          gCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCg                                                                                                                                       >  2:305128/1‑150 (MQ=255)
              ggCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAg                                                                                                                                   >  1:103473/1‑150 (MQ=255)
              ggCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCaaa                                                                                                                                   >  1:378972/1‑149 (MQ=38)
                ctggcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggc                                                                                                                                 <  1:955749/150‑1 (MQ=255)
                  ggcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCCGCCTTGCGGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCaaagcag                                                                                                                               <  2:378972/150‑6 (MQ=37)
                      tggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGcccc                                                                                                                           >  2:74618/1‑150 (MQ=255)
                      tggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGcc                                                                                                                             <  1:918705/148‑1 (MQ=255)
                      tggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGcc                                                                                                                             >  2:918705/1‑148 (MQ=255)
                                       ttATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcg                                                                                                                                <  2:122358/128‑1 (MQ=255)
                                       ttATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcg                                                                                                                                >  1:122358/1‑128 (MQ=255)
                                        tATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGAGGCTGCTTTGCAGCCCTTAATATATTACTCCGTGAAGTGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                                                                                                         <  1:467027/150‑1 (MQ=255)
                                            aaaTGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATa                                                                                                     >  2:331604/1‑150 (MQ=255)
                                             aaTGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATAc                                                                                                    <  1:354730/150‑1 (MQ=255)
                                              aTGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGTTGTGCAGCCCGCGAAGGGCCCCAAGGCGGCCCCGGTTAATTCAATCGTATTGATAcc                                                                                                   >  1:11694/1‑150 (MQ=255)
                                                   cGCTTGCCAGCGTCCGCGGGCACTGGCTTCCTGAAGCCAGGAGATGCCGCTCAGATAAGGTAATTGGGGCTGCTTAGCACCCCTTGGCGGGCGAGCTGGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCt                                                                                              <  2:239721/150‑1 (MQ=255)
                                                          cAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGaa                                                                                                                                             <  2:937448/96‑1 (MQ=255)
                                                          cAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGaa                                                                                                                                             >  1:937448/1‑96 (MQ=255)
                                                           aGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGc                                                                                                >  1:624986/1‑140 (MQ=255)
                                                           aGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcc                                                                                      >  2:1083628/1‑150 (MQ=255)
                                                           aGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCCGCCCCGGTTCATTCAATCGTATTGATACCCGc                                                                                                <  2:624986/140‑1 (MQ=255)
                                                            gCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGttcatt                                                                                                                   <  2:128073/120‑1 (MQ=255)
                                                            gCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGttcatt                                                                                                                   >  1:128073/1‑120 (MQ=255)
                                                                      gCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTa                                                                                                           <  1:719771/118‑1 (MQ=255)
                                                                      gCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTa                                                                                                           >  2:719771/1‑118 (MQ=255)
                                                                            gCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGtc                                                                     >  1:1124670/1‑150 (MQ=255)
                                                                             cTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGc                                                                                                >  1:645636/1‑122 (MQ=255)
                                                                             cTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGc                                                                                                <  2:645636/122‑1 (MQ=255)
                                                                             cTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGtct                                                                    <  1:1057809/150‑1 (MQ=255)
                                                                               tCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCg                                                                  <  2:1124670/150‑1 (MQ=255)
                                                                                 cTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGaa                                                                >  1:277235/1‑150 (MQ=255)
                                                                                        cAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTGGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAAGCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGAc                                                         <  1:305128/150‑1 (MQ=255)
                                                                                            aaatgccGCTCAGATAATGTAAGTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGAGGCTGCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCaa                                                     <  1:74618/145‑1 (MQ=255)
                                                                                                tCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGCGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAg                                                 >  2:709724/1‑150 (MQ=255)
                                                                                                  cGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGt                                               >  2:1004518/1‑150 (MQ=255)
                                                                                                         aTAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGa                                        >  2:729958/1‑150 (MQ=255)
                                                                                                                      ggCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCggg                           >  1:106670/1‑150 (MQ=255)
                                                                                                                         tgctTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATaa                                                                        >  2:857264/1‑102 (MQ=255)
                                                                                                                         tgctTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATaa                                                                        <  1:857264/102‑1 (MQ=255)
                                                                                                                                 aGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAgg                                                           <  2:1074783/107‑1 (MQ=255)
                                                                                                                                 aGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAgg                                                           >  1:1074734/1‑107 (MQ=255)
                                                                                                                                 aGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAgg                                                           >  1:1074783/1‑107 (MQ=255)
                                                                                                                                 aGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAgg                                                           <  2:1074734/107‑1 (MQ=255)
                                                                                                                                   cgcTTCGCGGGCTAGCCCAGGAAGGGCCGCAACGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCg              >  2:592785/3‑150 (MQ=255)
                                                                                                                                   cccTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCg              <  2:103473/150‑1 (MQ=255)
                                                                                                                                    ccTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCAc                               <  1:827398/132‑1 (MQ=38)
                                                                                                                                    ccTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCAc                               >  2:827398/1‑132 (MQ=38)
                                                                                                                                     cTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGcc            >  1:342428/1‑150 (MQ=255)
                                                                                                                                          gggtcGAGCCCGGGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCaaa       <  2:252684/146‑1 (MQ=255)
                                                                                                                                          cgGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCTGCCCCaaa       <  2:106670/150‑1 (MQ=255)
                                                                                                                                             gCGAGCGCGAGAAGGGCCCCAAGGCGGCCCCGGTACATTCACGCTTCTTGACACCCTCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCACGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCa    <  2:68278/150‑1 (MQ=255)
                                                                                                                                               gAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACt  >  1:128966/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACT  >  NC_002947/3048504‑3048796

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