Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,149,437 T→G 10.0% intergenic (‑423/+25) PP_mr51 ← / ← PP_4535 C4 ncRNA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,149,4370TG10.0% 172.1 / 11.4 60intergenic (‑423/+25)PP_mr51/PP_4535C4 ncRNA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base G (3/3);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

CGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGG  >  NC_002947/5149298‑5149583
                                                                                                                                           |                                                                                                                                                  
cgcgTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATccccccc                                                                                                                                           >  1:269635/1‑149 (MQ=255)
     cTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATcc                                                                                                                                                <  1:156263/139‑1 (MQ=255)
     cTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATcc                                                                                                                                                >  2:156263/1‑139 (MQ=255)
       tGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTg                                                                                                                                     >  4:2971/1‑148 (MQ=255)
        gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGtt                                                                                                                                      <  4:43743/146‑1 (MQ=255)
        gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGtt                                                                                                                                      <  5:103407/146‑1 (MQ=255)
        gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTgg                                                                                                                                    <  5:200044/148‑1 (MQ=255)
        gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTgg                                                                                                                                    >  3:278660/1‑148 (MQ=255)
              cTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTAt                                                                                                                                                  <  8:207129/128‑1 (MQ=255)
              cTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTAt                                                                                                                                                  >  7:207129/1‑128 (MQ=255)
                 tgCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGTTg                                                                                                                                     <  7:106799/138‑1 (MQ=39)
                 tgCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGTTg                                                                                                                                     >  8:106799/1‑138 (MQ=39)
                     ttCTGAATGCCCCCGAGCAATCCGCCCCATAAGACATAATGTGGCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGGGGCCCAGCCACTCCTCATCCACGTTCCAACCAGGAGGCTAACCCCCCCAGGATGGACCGTCAGTGAcgc                                                                                                                      >  7:304907/1‑147 (MQ=255)
                     ttCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGtt                                                                                                                                      <  1:144109/133‑1 (MQ=39)
                     ttCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGtt                                                                                                                                      >  2:144109/1‑133 (MQ=39)
                        tGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAACCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCATGGTTGGAACTTCAGTTACTCgg                                                                                                                    >  3:107248/1‑148 (MQ=255)
                          aaTGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGc                                                                                                                   <  8:13479/147‑1 (MQ=255)
                            tGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATc                                                                                                               <  3:2971/149‑1 (MQ=255)
                              aCCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAAt                                                                                                                >  6:322474/1‑146 (MQ=255)
                                cACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGAc                                                                                                           >  3:151038/1‑149 (MQ=255)
                                 aCGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGAcc                                                                                                          <  5:322474/149‑1 (MQ=255)
                                   gAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCAAAGCCAATCCTCATCCAAGTTCCAACCAGGAGCCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACcgc                                                                                                        >  7:5205/1‑149 (MQ=255)
                                       aaTCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGt                                                                                                    >  1:165419/1‑149 (MQ=255)
                                        aTCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGtt                                                                                                   <  2:165419/149‑1 (MQ=255)
                                         tCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGt                                                                                                                           >  1:181034/1‑124 (MQ=255)
                                         tCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGt                                                                                                                           <  2:181034/124‑1 (MQ=255)
                                            gACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTtctc                                                                                                >  7:248798/1‑148 (MQ=255)
                                                                tCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTc                                                                            >  3:5523/1‑148 (MQ=255)
                                                                tCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTc                                                                            <  4:5523/148‑1 (MQ=255)
                                                                         cccgccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACa                                                                  >  4:181153/1‑149 (MQ=255)
                                                                                tatCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTAGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGc                                                           >  5:128391/1‑149 (MQ=255)
                                                                                 atCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAAGCTGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATATTTGcc                                                          <  8:52370/149‑1 (MQ=255)
                                                                                    ttGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAAGCTGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTcc                                                                           <  6:195366/129‑1 (MQ=255)
                                                                                    ttGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCAAAGTTCCAACCAGGAGGCTATACCCCCCTGGTTGGAACTTAAGTTACTAGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTcc                                                                           >  5:195366/1‑129 (MQ=255)
                                                                                            aaCTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCa                                               <  8:5205/149‑1 (MQ=255)
                                                                                            aaCTGCAGTGTCCCAGCCAATCCTCATCCAACTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCa                                               <  4:107248/149‑1 (MQ=255)
                                                                                                        cAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTcc                                                                           <  1:287826/109‑1 (MQ=255)
                                                                                                        cAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTcc                                                                           >  2:287826/1‑109 (MQ=255)
                                                                                                        cAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                   >  3:54245/1‑149 (MQ=255)
                                                                                                        cAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                   >  3:58593/1‑149 (MQ=255)
                                                                                                         aaGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTg                                  >  8:218710/1‑149 (MQ=255)
                                                                                                           gCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGAt                                <  4:151038/149‑1 (MQ=255)
                                                                                                             cAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAAGTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCaa                                                    <  2:250279/127‑1 (MQ=255)
                                                                                                             cAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCGTCGAAGACATAGTTGCCCTGCaa                                                    >  1:250279/1‑127 (MQ=255)
                                                                                                              aaTCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAg                                                               >  7:36423/1‑115 (MQ=255)
                                                                                                              aaTCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAg                                                               <  8:36423/115‑1 (MQ=255)
                                                                                                              aaTCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCg                              >  6:30582/1‑148 (MQ=255)
                                                                                                                     aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATa                                                                >  3:272584/1‑107 (MQ=255)
                                                                                                                     aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATa                                                                <  4:272584/107‑1 (MQ=255)
                                                                                                                     aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAg                                              <  3:222669/125‑1 (MQ=255)
                                                                                                                     aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAg                                              >  4:222669/1‑125 (MQ=255)
                                                                                                                     aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGcc                                           >  3:147678/1‑128 (MQ=255)
                                                                                                                     aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGcc                                           <  4:147678/128‑1 (MQ=255)
                                                                                                                     aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACaa                       <  4:278660/148‑1 (MQ=255)
                                                                                                                            ttCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTAGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGc               <  6:128391/149‑1 (MQ=255)
                                                                                                                              ccAACCAGGAGGCTATCCCCCCCTGGTTGAAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCgg             >  4:112932/1‑149 (MQ=255)
                                                                                                                                aaCCAGGGGGGGATAGCCTCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt           >  2:123493/1‑149 (MQ=255)
                                                                                                                                aaCCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt           <  8:248798/149‑1 (MQ=255)
                                                                                                                                         gCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCgg  <  4:58593/149‑1 (MQ=255)
                                                                                                                                         gCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCgg  <  4:54245/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
CGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGG  >  NC_002947/5149298‑5149583

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: