Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,212,520 C→T 100% intergenic (‑54/+134) PP_1060 ← / ← PP_1061 glutamate synthase large subunit/ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,212,5200CT98.2% 181.8 / ‑4.4 56intergenic (‑54/+134)PP_1060/PP_1061glutamate synthase large subunit/ATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/23);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACCATTCGTGGTGGCTGAAGGGGAGGCTCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGACCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCGTTGCCGCCGCTACCGTA  >  NC_002947/1212376‑1212661
                                                                                                                                                |                                                                                                                                             
acaaTTCGTGGTGGCTGAAGGGGAGGCTCGCCAGGGTGAACAGCACACAGGCTGCGAAGAAGGCATAACGACTGAGAAGTGACAGTCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATa                                                                                                                                           <  6:329418/146‑1 (MQ=255)
aCCATTCGTGGTGGCTGAAGGGGAGGCTCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATa                                                                                                                                           <  2:409082/149‑1 (MQ=255)
       gtggtgGCTGAAGGGGAGGCTCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGt                                                                                                                                    >  1:168679/1‑149 (MQ=255)
           tggCTGAAGGGGAGGCTCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGcc                                                                                                                                >  7:354470/1‑149 (MQ=255)
           atcttGACGGGGAGGCTCGCCAGGGGGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGcc                                                                                                                                <  3:304814/145‑1 (MQ=255)
                 aaGGGGAGGCTCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCAc                                                                                                                          >  2:16255/1‑149 (MQ=255)
                      gAGGCTCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGCTCATAGCCAAGTTGCCGGGCACTGCa                                                                                                                      >  1:97396/1‑148 (MQ=255)
                       aGGCTCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCaaa                                                                                                                    <  8:354470/149‑1 (MQ=255)
                           tCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATa                                                                                                                                           >  1:356909/1‑122 (MQ=255)
                           tCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATa                                                                                                                                           <  2:356909/122‑1 (MQ=255)
                             gCCAGGGTTAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCaaaataaaa                                                                                                              <  3:262297/149‑1 (MQ=255)
                                  ggTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCggg                                                                                                         >  3:237451/1‑149 (MQ=255)
                                           cacacaGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAACTTGCCGGGCACTGCAAAATAAAATCGCGGTAGCGaa                                                                                                 >  7:291403/1‑148 (MQ=255)
                                                aGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCAt                                                                                           >  5:80643/1‑149 (MQ=255)
                                                    ggcgAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTa                                                                                                      >  3:226953/1‑134 (MQ=255)
                                                    ggcgAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTa                                                                                                      <  4:226953/134‑1 (MQ=255)
                                                    ggcgAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATtctc                                                                                       <  5:156495/149‑1 (MQ=255)
                                                      cgAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTACTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGTCAAACCACGATCATAGCAAAGTTGCAGGGCACTGCAAAATACAATCGGGGTATCGAAAATCATtatc                                                                                       >  3:98000/1‑147 (MQ=255)
                                                        aagaagGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGc                                                                                   <  2:43591/149‑1 (MQ=255)
                                                           aagGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCg                                                                                  >  6:213188/1‑147 (MQ=255)
                                                            agGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTa                                                                                <  4:237451/148‑1 (MQ=255)
                                                               cATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAg                                                                            >  7:75379/1‑149 (MQ=255)
                                                               cATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAg                                                                            >  5:301232/1‑149 (MQ=255)
                                                                 tAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGa                                                                                     >  6:61309/1‑138 (MQ=255)
                                                                 tAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGa                                                                                     <  5:61309/138‑1 (MQ=255)
                                                                 tAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATAATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCg                                                                          >  1:39701/1‑149 (MQ=255)
                                                                   aCGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCgg                                                                         <  5:213188/148‑1 (MQ=255)
                                                                           gaAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGATTCACCcg                                                                >  2:268609/1‑149 (MQ=255)
                                                                               tGACAGGCTCATACAGGGTCCTTGGGATGGCGCGGTCAGGCGCAGGGCATGGGGGCAAACCCCGATAATAGCCAAGTTGCAGGGCCGTGCTAATTTAAATCGGTGTAGCGAAAATCATTATCGAGCGGAGGCGCGGTGTCACCCTCGa                                                             >  6:327309/1‑148 (MQ=255)
                                                                               tGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGa                                                             >  2:380026/1‑148 (MQ=255)
                                                                                 aCAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGaaa                                                           <  3:196751/148‑1 (MQ=255)
                                                                                 aCAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGa                                                             >  4:196751/1‑146 (MQ=255)
                                                                                 aCAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAg                                                          >  6:237534/1‑149 (MQ=255)
                                                                                 aCAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAg                                                          >  6:400071/1‑149 (MQ=255)
                                                                                    ggCTCATACAGTGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAgggc                                                       <  6:275868/149‑1 (MQ=255)
                                                                                      cTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCATGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAgggccg                                                     <  6:318224/149‑1 (MQ=255)
                                                                                      cTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTAACCCGCGAAAgggccg                                                     >  5:318224/1‑149 (MQ=255)
                                                                                          tACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCCCCCGCGAAAgggccgggcc                                                 >  8:371114/1‑149 (MQ=255)
                                                                                               ggTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCctgctg                                            <  2:97396/149‑1 (MQ=255)
                                                                                                    ttGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAgg                                                         >  7:27869/1‑131 (MQ=255)
                                                                                                    ttGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAgg                                                         <  8:27869/131‑1 (MQ=255)
                                                                                                       cGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAAtc                                     <  3:364917/148‑1 (MQ=255)
                                                                                                          tGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCa                                 <  1:16255/149‑1 (MQ=255)
                                                                                                                    aGGCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTCCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTAACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACccgccgag                        >  8:295291/1‑146 (MQ=255)
                                                                                                                      gCTCAGGGCCTGGGGGCAAGCCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCAc                     >  4:144605/1‑149 (MQ=255)
                                                                                                                      gCTCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCGTTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCAc                     >  5:213119/1‑149 (MQ=255)
                                                                                                                      gCTCAGGGCCTGGGGGCAAACCCCGACCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCAc                     >  5:393043/1‑149 (MQ=255)
                                                                                                                       cTCAGGGCCTGGGGTCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCAcc                    <  5:37227/149‑1 (MQ=255)
                                                                                                                        tCAGGGCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAAGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCg                   >  1:183038/1‑149 (MQ=255)
                                                                                                                          atggCCTGTGGGCAAACCCCGATCTTAGCGAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCAGCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCGtt                 <  7:295291/147‑1 (MQ=37)
                                                                                                                            ggCCTGGGGGCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCGTTgc               >  7:244732/1‑149 (MQ=255)
                                                                                                                                  ggggCAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCGTTGCcgccgc         >  6:335234/1‑149 (MQ=255)
                                                                                                                                      cAAACCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCAAATCTCCACCCGCCGTGCACCGTTGCCGCCGCTAcc     <  2:39701/149‑1 (MQ=255)
                                                                                                                                        aaCCCCGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCGTTGCCGCCGCTACCg    >  5:92473/1‑148 (MQ=255)
                                                                                                                                            ccGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCGTTGCCGCCGCTACCGTa  >  2:102622/1‑146 (MQ=255)
                                                                                                                                            ccGATCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCGTTGCCGCCGCTACCGTa  <  1:102622/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
ACCATTCGTGGTGGCTGAAGGGGAGGCTCGCCAGGGTGAACAGCACACAGGCGGCGAAGAAGGCATAACGACTGAGAAGTGACAGGCTCATACAGGGTCCTTGCGATGGCGCGGTCAGGCTCAGGGCCTGGGGGCAAACCCCGACCATAGCCAAGTTGCCGGGCACTGCAAAATAAAATCGGGGTAGCGAAAATCATTCTCGAGCGTAGGAGCGGGTTCACCCGCGAAAGGGCCGGGCCTGCTGGCCAATCTCCACCCGCCGTGCACCGTTGCCGCCGCTACCGTA  >  NC_002947/1212376‑1212661

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: