Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,761,329 A→C 26.6% E96A (GAG→GCG)  CJ019_03574 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,761,3290AC26.6% 69.9 / 4.0 60E96A (GAG→GCG) CJ019_03574hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base C (15/3);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 1.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGG  >  ACB122/3761200‑3761462
                                                                                                                                 |                                                                                                                                     
gCATGGTCCAAGGTGAGCCAGACTTCAGCGAGCCGACCCTGTCTGCCGGCGGGGACTGCGTGATCGCCAGCCTGTTCAGCGCCACCCTGCTGGCCGTCTTTTTCGGCCGGCTGACCCTGCGCGCCGGCGAGCCgg                                                                                                                                  <  6:201801/135‑1 (MQ=255)
gCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCgg                                                                                                                                  <  2:245014/135‑1 (MQ=255)
 cATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCCCACTGCTGGCCGTCTTTTTAGGCCTGCTGACACATCTCGCCCGCCCGCCggg                                                                                                                                 >  3:122004/1‑135 (MQ=255)
 cATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGCCCCTGCTCGCCGGCGAGCCggg                                                                                                                                 >  4:253184/1‑135 (MQ=255)
 cATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTAGGCCTGCTGACCCTGCTCGCCCGCGCGCCggg                                                                                                                                 >  5:192704/1‑135 (MQ=255)
 cATGGTCCAAGGTGAGCAAGACTTCAGAGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCTGCCAGCTGACCCTGCTCCCCCGCGAGcccgc                                                                                                                                 >  8:292804/1‑132 (MQ=255)
    ggTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAg                                                                                                                              <  1:123445/135‑1 (MQ=255)
     gTCCAAGGCGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGCTCGCCAGCCTGTTCCGCGCCACACTGCTGGCCGTCTTTTTCGGCCGGCTGACCCTGCTCGCCGGCGAGCCGGGCAGc                                                                                                                             <  7:172852/135‑1 (MQ=255)
        cAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCCGGCcccccc                                                                                                                          >  1:3525/1‑135 (MQ=255)
        aaaGGTGAGGAAGACTCCAGCGATCCCAGCCTGGCTGCCGGCGGTTACTGCGGGATCGCCAGCGTGTTCAGCGCCACACTGCTGGCGGTCTTTTTCGGCCTGCTGACCCTGCTGGCCGGGGAGCCGGGCAGcccc                                                                                                                          <  2:346810/134‑1 (MQ=255)
         aaGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCGTTTTCGGCCTGCTGACCCAGCTCGCCCCCGAGCCCGGGCGccccg                                                                                                                         >  5:361179/1‑134 (MQ=255)
              gAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCGGGGTTTTGGGCCCTCGGGCCCTGCCACCCGGCCCGCCCGGCCGCCCCCTAcccc                                                                                                                   >  3:330172/1‑134 (MQ=255)
                  aaGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCCGCCCCGTACccaac                                                                                                                 >  7:365129/1‑134 (MQ=255)
                  aaGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCCGCCCCGGacccaccc                                                                                                                >  3:311186/1‑135 (MQ=255)
                  aaGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCAGCCCCATACccaacc                                                                                                                >  6:440584/1‑135 (MQ=255)
                  aaGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACccaacc                                                                                                                >  8:215071/1‑135 (MQ=255)
                  aaGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGAGCCCGGCAGCCCCATaaccaacc                                                                                                                >  7:64656/1‑135 (MQ=255)
                   aGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCAGCCCCTTACccaacca                                                                                                               >  2:146696/1‑135 (MQ=255)
                   aGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACccaaccc                                                                                                               >  4:21559/1‑134 (MQ=255)
                   aGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACccaacca                                                                                                               >  7:465237/1‑135 (MQ=255)
                   aGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCACCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTGACccaacca                                                                                                               >  3:375092/1‑135 (MQ=255)
                       ttCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGCGCCCCGCCGCCCCAGACcccacccccca                                                                                                           >  5:403183/1‑130 (MQ=255)
                        tCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCCGCCCCGGACCCCACCCACCTc                                                                                                          >  5:169899/1‑135 (MQ=255)
                         cAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGAACCTGCGACCCCGCGCGCCCGGCCCCCCCCGACCCACCCAACCGcc                                                                                                         >  2:430918/1‑135 (MQ=255)
                            cgtaccGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTGCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGcc                                                                                                      <  8:182421/131‑1 (MQ=255)
                            cgctCCGACCCTGTCTGCCGGCGGTGACTGCGTGAGCGCCAGCCTGTTCCGCGCCACACTGCGGGCCGTCTTTTTCGGCCTGCTGACCCTGCGCGCGGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGcc                                                                                                      <  5:343457/132‑1 (MQ=255)
                              ctCCGACCCTGGCGGCCGGCGGTGACGGCGTGCTCGCCAGCCTGCTCAGCCCCACCCTGCGGGCCGTGTTTTTGGGCGTGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCAc                                                                                                    <  1:55278/134‑1 (MQ=255)
                                 cgccgcTGTCTGCCGGGGGTGCGTGGGTGATCGCCAGCCTGGTCAGCGCCACACTGCGGGCCGTCTTTTTCGGCCGGCTGACCCTGCGCGCCGGCGCGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGcc                                                                                                 <  5:70811/130‑1 (MQ=255)
                                  gACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCCGCCCCGTACCCAACCAACCTCCCCCCCCCCt                                                                                                >  2:10363/1‑135 (MQ=255)
                                      cTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGccc                                                                                            <  2:156098/135‑1 (MQ=255)
                                           ggCCGGCGGTGACTGCGTGAGCGCCAGCCTGTTCCGCGCCACACTGCGGGCCGTCTTTTTCGGCCTGCTGACCCTGCGCCCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAagca                                                                                       <  7:408743/134‑1 (MQ=255)
                                                 cggTGACTGCGTGAGCGCCAGCCTGTTCAGCGCCACACTGCGGGCCGTCTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTg                                                                                 <  4:237348/135‑1 (MQ=255)
                                                             gATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCCGCCCCGTTCCCAACCAACCTACGCCCCCCCGGCCCCAGCCCCGCTGGGCCCTgggcca                                                                     >  8:6655/1‑131 (MQ=255)
                                                             gATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCCGCCCCGTACCCCACCACCCTCCGCCCCGCCTGCCCCAGCAGCGCGGGCCCCTGGTCac                                                                     >  7:160542/1‑134 (MQ=255)
                                                              agcGCCAGCCTGTTCACCGCCACACTGCGGGCCGTCTTTTGCGGCCTGCGGACCCGGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAAt                                                                    <  1:23442/133‑1 (MQ=255)
                                                              aTCGCCAGCCTGCTCAGCGCCACACTGCGGGCCGTCTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAAt                                                                    <  1:250477/135‑1 (MQ=255)
                                                                cGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATcc                                                                  <  2:167535/135‑1 (MQ=255)
                                                                cGCCAGCCTGTGCCGGGCCACTGTGTGGGGGGTGTTTTGCGGCTTGCTGCCGGCGCGCCCCGGCGCGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATcc                                                                  <  3:457720/135‑1 (MQ=255)
                                                                             aGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTcc                                                      <  8:465237/134‑1 (MQ=255)
                                                                              gcgcCACCCTGCGGGCCGTCTTTTTCGGCCTGCGGACCCGGCGCGCCGCCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcg                                                    <  5:314199/135‑1 (MQ=255)
                                                                                  cacaCTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCCCCCTCCGCCCCGCCTGCCCAAGCAGCGCTGGCCACGGGTCAATCCCCGCGCCTCTCCcgcttc                                                >  8:108832/1‑135 (MQ=255)
                                                                                  cacaCTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttc                                                >  4:145274/1‑135 (MQ=255)
                                                                                  cacaCTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCCGCGCGCCGGGCCGCCCCGTACCCCACCCCACTACGCCCCGCCGGCCCCAGCCGCGCGGGCCCCTGGTGCATCCCCCCCCCGCGCCcgcgtc                                                >  1:24141/1‑135 (MQ=255)
                                                                                  cacaCTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCGCGCCGGCGAGCCGGGCCGCCCCGTACCCCACCAACCTCCGCCCCCCCTGCCCAACCCGCGCGGGCCACTGGTCAACCCCCGCCCCGCGcccccgtc                                                >  8:176167/1‑130 (MQ=255)
                                                                                     cggcggggCCGTTTTTTGCGCCCTGCGGCCGCGGCGCGCGGGCGCGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCcgcttcgct                                             <  3:397400/129‑1 (MQ=255)
                                                                                            gCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGCGCCGGGCAGCCCCTTACCCCACCAACCTCCGCCACGCCGGCCCAAGCAGCGCGGGCCCCTGGTCCATCCCCCCGCCTCTCCCGCGTCGCGTCTcgcg                                      >  6:431031/1‑135 (MQ=255)
                                                                                            gCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCATACCCAACCAACCTCCGCCACGCCTGCCCAAGCCGCGCTTGCCCCTGGTCAATACCCGCGCCTCTCCCGCTTAGCTTCTcgcg                                      >  2:409235/1‑135 (MQ=255)
                                                                                            gCCGTCTTTTTCGGCCTGCTGACCCTGCTCCCCGGCGAGCCGGGCAGCCCCCTTCCCCACCAACCACCCCCCCCCCAGCCCAAGCCGCGCAGGCCACAGGTGAAACCCCGCGCCGCGCCCGCGTCGCATcgcgcg                                      >  6:44011/1‑135 (MQ=255)
                                                                                            gCCGTCTTTTTCGGCCTGCTGACCCTGCACGCCGGCGAGCCGGGCAGCCCCGTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTCGTCAATCCCCGCCCCTCGCCCGCATCGCTTCTcgcg                                      >  2:87106/1‑135 (MQ=255)
                                                                                               gTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGTGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTc                                   >  2:21619/1‑135 (MQ=255)
                                                                                                       cGGCCCGCGGCCCCTGCGCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCg                           <  3:65076/135‑1 (MQ=255)
                                                                                                        ggCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCATGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCg                           <  4:41967/134‑1 (MQ=255)
                                                                                                        ggCCTGCGGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCg                           <  6:144000/134‑1 (MQ=255)
                                                                                                            ggctgACCCTGCTGGCTGGGGAGCGGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCTCTGTCCACTGGTCAATCCCAGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACg                      <  4:122004/134‑1 (MQ=255)
                                                                                                                  cccTGCTCGCCGGTGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTcccc                 <  1:21619/134‑1 (MQ=255)
                                                                                                                      gCTCGCCGGCGAGCCGGGCAGCCCCGTCCCCAACCCACCTCCGCCCCGCCTGCCCAAGCCGCGCTGGCCACTGGTCAATCCCCGCGCCGCTCCCGCTTCGCTTCGCGCGGTCTGACATCGCCACGTCCCCGTgcc            >  8:200407/1‑135 (MQ=255)
                                                                                                                       cTCGCCGGCGCGCCGGGCAGCCCCTTACCCCACCCACCTCCGCCCCGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCGCGCGGTCTGACATCGCCACGTCCCCGTgccg           >  4:57235/1‑135 (MQ=255)
                                                                                                                       cTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGCCCCCGTgccg           >  8:81782/1‑135 (MQ=255)
                                                                                                                          gCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCGTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCt        >  1:138531/1‑135 (MQ=255)
                                                                                                                                gAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCgg  <  1:146696/135‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                     
GCATGGTCCAAGGTGAGCAAGACTTCAGCGATCCGACCCTGTCTGCCGGCGGTGACTGCGTGATCGCCAGCCTGTTCAGCGCCACACTGCTGGCCGTCTTTTTCGGCCTGCTGACCCTGCTCGCCGGCGAGCCGGGCAGCCCCTTACCCAACCAACCTCCGCCACGCCTGCCCAAGCAGCGCTGGCCACTGGTCAATCCCCGCGCCTCTCCCGCTTCGCTTCTCGCGTTCTGACATCGCCACGTCCCCGTGCCGCCTGATCGG  >  ACB122/3761200‑3761462

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