Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,011,479 T→A 8.4% intergenic (‑417/‑110) PP_0871 ← / → prfC glycine betaine/carnitine/choline ABC transporter permease/peptide chain release factor 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,011,4790TA 8.4% 165.7 / 4.5 59intergenic (‑417/‑110)PP_0871/prfCglycine betaine/carnitine/choline ABC transporter permease/peptide chain release factor 3
Reads supporting (aligned to +/- strand):  ref base T (27/27);  new base A (2/3);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ACTGCGTTTTGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAGCGCC  >  NC_002947/1011340‑1011622
                                                                                                                                           |                                                                                                                                               
aCTGCGTTTTGTTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCTGCCTCTTCGCTGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                        <  5:47107/149‑2 (MQ=255)
    cGTTTTGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAgg                                                                                                                                     <  2:125954/148‑1 (MQ=255)
    cGTTTTGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAg                                                                                                                                      <  5:352058/147‑1 (MQ=255)
     gTTTTGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAgggg                                                                                                                                   >  1:95682/1‑149 (MQ=255)
         tGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCgg                                                                                                                               >  4:75046/1‑149 (MQ=255)
              aTTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAg                                                                                                                          >  8:180285/1‑149 (MQ=255)
              aTTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAg                                                                                                                          >  8:175280/1‑149 (MQ=255)
                   tCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTAACTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGgcgc                                                                                                                      >  1:80773/1‑148 (MQ=255)
                      tGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCcatt                                                                                                                  <  1:355426/149‑1 (MQ=255)
                       gATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCcatt                                                                                                                  >  4:154292/1‑148 (MQ=255)
                          aaaGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCcattcatt                                                                                                              <  7:180285/149‑1 (MQ=255)
                          aaaGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCcattcatt                                                                                                              <  7:175280/149‑1 (MQ=255)
                                 ttGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCgg                                                                                                           >  7:245111/1‑145 (MQ=255)
                                       gTTGGGTTTGTTTTGTTTGTGACGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGAATCCGAACCCTGTGTAGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGTCCGGCACAGTCGCCATTCATTCGGAAGATTTTgg                                                                                                 >  1:109979/1‑149 (MQ=38)
                                          gggTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTc                                                                                               <  3:75046/148‑1 (MQ=255)
                                          gggTTTGTTTTGTTTGGGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACCAAGACTCCACCGAATCCCAACCCTGTGTCGAACCTGTCGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTc                                                                                               <  3:2279/148‑1 (MQ=38)
                                            gtttgtttTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACATGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAg                                                                                            <  4:285599/149‑1 (MQ=255)
                                                 ttttgtttGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGccc                                                                                        >  7:24337/1‑148 (MQ=255)
                                                  tttgtttgTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGTCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCa                                                                                       >  4:185992/1‑148 (MQ=255)
                                                     gtttgtGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCa                                                                                       >  6:19794/1‑145 (MQ=255)
                                                     gtttgtGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCa                                                                                       <  5:19794/145‑1 (MQ=255)
                                                         gtgCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                                       <  7:224523/93‑1 (MQ=255)
                                                         gtgCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGaa                                                                                                                                       >  8:224523/1‑93 (MQ=255)
                                                           gCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTa                                                                             >  7:88518/1‑149 (MQ=255)
                                                            ccGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTAt                                                                            >  1:49175/1‑149 (MQ=255)
                                                                ccgctTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATa                                                                        <  2:80773/146‑1 (MQ=255)
                                                                        cgGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATa                                                                        >  1:193415/1‑141 (MQ=255)
                                                                        cgGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATa                                                                        <  2:193415/141‑1 (MQ=255)
                                                                           gTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGTCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTc                                                             <  3:185992/149‑1 (MQ=255)
                                                                           gTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTc                                                             <  2:95682/149‑1 (MQ=255)
                                                                             aaaCTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCa                                                            >  8:172171/1‑148 (MQ=255)
                                                                                cTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATc                                                        >  7:300360/1‑149 (MQ=255)
                                                                                  cccTCCTACGAAGAATGCACCGCATCCGAACCCTTTGTCGTACCTGTAGGCGCGGTTATACCCGATAAGGTGCCGGCACAGGCGACATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATccc                                                      <  2:109979/147‑1 (MQ=255)
                                                                                   gCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATcccc                                                     <  8:300360/149‑1 (MQ=255)
                                                                                      ccTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTg                                                                                                  >  7:7398/1‑101 (MQ=255)
                                                                                      ccTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTg                                                                                                  <  8:7398/101‑1 (MQ=255)
                                                                                       cTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCa                                                   <  3:154292/147‑1 (MQ=255)
                                                                                       cTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGt                                                 >  4:324291/1‑149 (MQ=255)
                                                                                           gAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAgg                                             <  8:24337/149‑1 (MQ=255)
                                                                                            aaGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCAATCGGCAGATTTTGGTCCAGGCCCAACTTc                                                                                  >  2:319901/1‑111 (MQ=255)
                                                                                                  gCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCACCAGTCAGGCGATTTc                                      >  6:344612/1‑149 (MQ=255)
                                                                                                   cACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTAccc                                                                                                                                            <  7:217543/46‑1 (MQ=255)
                                                                                                   cACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTAccc                                                                                                                                            >  8:217543/1‑46 (MQ=255)
                                                                                                       gCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATg                                 <  3:324291/149‑1 (MQ=255)
                                                                                                             gAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATcc                                                       <  8:93471/121‑1 (MQ=255)
                                                                                                                cccTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAg                        <  7:172171/149‑1 (MQ=255)
                                                                                                                cccTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAg                        >  7:90679/1‑149 (MQ=255)
                                                                                                                cccTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAg                        <  8:88518/149‑1 (MQ=255)
                                                                                                                   tgtgTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGgcc                     >  5:275299/1‑149 (MQ=255)
                                                                                                                    gtgtCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCAt                                  <  6:188876/135‑1 (MQ=255)
                                                                                                                    gtgtCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCAt                                  >  5:188876/1‑135 (MQ=255)
                                                                                                                            ccTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGt             <  2:49175/148‑1 (MQ=255)
                                                                                                                             cTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTcg           <  6:275299/149‑1 (MQ=255)
                                                                                                                               gTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTcgc          >  6:372246/1‑148 (MQ=255)
                                                                                                                               cccttcgcgGGTAAACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTcgc          >  5:68803/6‑148 (MQ=255)
                                                                                                                                  ttcgcgGGTAAACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGac                               >  6:153514/3‑124 (MQ=255)
                                                                                                                                  ttcgcgGGTAAACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGac                               <  5:153514/122‑1 (MQ=255)
                                                                                                                                         gTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAGcgcc  >  7:317718/1‑146 (MQ=255)
                                                                                                                                         gTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAGcgcc  <  8:317718/146‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                               
ACTGCGTTTTGGTTATTTGTCGTGATAAAGCGATTGCCTGTTGGGTTTGTTTTGTTTGTGCCGGCCTCTTCGCGGGTAAACTCGCTCCTACGAAGACTGCACCGCATCCGAACCCTGTGTCGTACCTGTAGGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGGCGCCATTCATTCGGCAGATTTTGGTCCAGGCCCAACTTCAGGTATGATACTGCCCCTTTCAAATCCCCCAGTCAGGCGATTTCCCATGACCAACCAGGCCGCCGAAGTCGCGAAGCGCC  >  NC_002947/1011340‑1011622

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