Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,796,627	C→T	11.6%	intergenic (‑185/+33)	pvdH ← / ← PP_4224	diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,796,627	0	C	T	11.6%	132.6 / 7.5	56	intergenic (‑185/+33)	pvdH/PP_4224	diaminobutyrate‑2‑oxoglutarate transaminase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (21/26);  new base T (3/3);  total (27/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 = 2.91e-01

GACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGCCC  >  NC_002947/4796480‑4796773                                                                                                                                                    |                                                                                                                                                    gACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCa                                                                                                                                                    <  2:268594/149‑1 (MQ=255)  aCATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTAc                                                                                                                                                   <  3:341712/149‑1 (MQ=38)  aCATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGCCTGTGCTGGCCTCTTCGGAGGTGAACCCGCTCCTACGGGCGGGTTTAc                                                                                                                                                   <  1:217733/149‑1 (MQ=37)    aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAccc                                                                                                                                                 >  7:108124/1‑149 (MQ=255)    aTGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAccc                                                                                                                                                 <  3:256467/149‑1 (MQ=255)      gggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTGACCcgc                                                                                                                                               >  4:72417/1‑149 (MQ=255)       ggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                               >  1:19112/1‑148 (MQ=255)       ggTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCcgc                                                                                                                                               <  2:53631/148‑1 (MQ=255)             tCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAgg                                                                                                                                        >  7:283188/1‑149 (MQ=255)             tCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAgg                                                                                                                                        >  3:76213/1‑149 (MQ=255)                       cAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTGGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAgg                                                                                                                              <  6:127852/149‑1 (MQ=255)                               tGGGCGGTTATCGCGTTTGAAACGCTTGACGGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCATCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑Ac                                                                                                                     <  7:344504/149‑1 (MQ=38)                                   cGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCa                                                                                                                 <  4:362282/149‑1 (MQ=255)                                    ggTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCaa                                                                                                                <  7:1057/149‑1 (MQ=255)                                           gCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCa                                                                                                          <  8:108124/148‑1 (MQ=255)                                                ggAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGAcgc                                                                                                    >  3:295543/1‑149 (MQ=255)                                                                ggggAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGcaac                                                                                     >  4:236667/1‑148 (MQ=255)                                                                   gaggaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACaaaa                                                                                 <  4:295543/149‑1 (MQ=255)                                                                    aggaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAc                                                                                                                                                   <  2:359006/83‑1 (MQ=255)                                                                    aggaTTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCAc                                                                                                                                                   >  1:359006/1‑83 (MQ=255)                                                                      gaTTTAGCGGTTGTTCGGCGATTTAATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGTATCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAAccc                                                                              <  7:133188/149‑1 (MQ=38)                                                                          tatcGGTTGTTAAGCGAATTTATTTTTGTCTGTGCTGGCCTCATCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGATGCCAGCACAGGCGATAATACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGc                                                                            <  1:259023/145‑1 (MQ=38)                                                                            gCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGc                                                                        >  5:8408/1‑149 (MQ=255)                                                                              ggTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCg                                                                      >  1:331290/1‑149 (MQ=255)                                                                              ggTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCg                                                                      <  4:263586/149‑1 (MQ=255)                                                                               gttgttACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTAACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGt                                                                     >  7:341564/1‑149 (MQ=255)                                                                                       gcgATTTTTTTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACgg                                                             <  5:272261/149‑1 (MQ=255)                                                                                       gcgATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAGCCCGCCCGTAGGAGCGGGTTCACCCGGGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCATCACTACgg                                                             >  7:123485/1‑149 (MQ=255)                                                                                         gATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCt                                                           >  3:275197/1‑149 (MQ=255)                                                                                         gATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCt                                                           <  4:76213/149‑1 (MQ=255)                                                                                         gATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAACAGGCCAGCACAGGCGATAA‑ACATCAAAGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAACCCGACACTACgg                                                             >  4:91341/1‑147 (MQ=255)                                                                                            tttattttTGTCTGTGCTGGCCTCATCGGGGGTAAACCCGCCGGTAGGCGCGGGTTGGCCCGGGAAGGGGCCCGCACAGGGGATTA‑ACATGAAGGCGCAGGCCGGGCCGGGACCTGCGAGAAGACCCGCAAGCCGTGACTACgggggt                                                         >  6:202680/1‑144 (MQ=255)                                                                                             ttattttTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTc                                                        <  2:61298/148‑1 (MQ=255)                                                                                                     gccTGTGCTGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAg                                               >  4:66503/3‑149 (MQ=25)                                                                                                      tCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAgg                                              >  1:194964/1‑149 (MQ=255)                                                                                                      tCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAgg                                              <  2:19112/149‑1 (MQ=255)                                                                                                      tCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAg                                               <  8:283188/148‑1 (MQ=255)                                                                                                          tgCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGc                                             <  4:275197/146‑1 (MQ=255)                                                                                                             tGGCCTCTTCGCGGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCgg                                       >  8:33615/1‑149 (MQ=37)                                                                                                               gCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCgg                                       <  3:236667/147‑1 (MQ=255)                                                                                                                        cgGGTGAACCCGCTCCTACGGGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTg                             >  5:305151/1‑148 (MQ=37)                                                                                                                        cgGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGc                            <  2:194964/149‑1 (MQ=255)                                                                                                                         gggTGAACCCGCTCCTACGGGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCa                           >  6:55399/1‑149 (MQ=38)                                                                                                                           gTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCAc                          >  3:361876/1‑148 (MQ=255)                                                                                                                              aaCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc                      <  5:202680/149‑1 (MQ=255)                                                                                                                              aaCCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATc                      >  3:55315/1‑149 (MQ=255)                                                                                                                                cccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGc                    >  6:26445/1‑149 (MQ=255)                                                                                                                                 ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                   >  3:218211/1‑149 (MQ=255)                                                                                                                                 ccgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCAGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGcc                   >  8:189244/1‑149 (MQ=255)                                                                                                                                  cgcccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCt                  <  6:8408/149‑1 (MQ=255)                                                                                                                                     ccgTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGca               >  6:194666/1‑149 (MQ=255)                                                                                                                                          ggAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCATCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATc          <  8:123485/149‑1 (MQ=255)                                                                                                                                          cgggCGGGTTTACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATc          <  7:33615/146‑1 (MQ=255)                                                                                                                                           gAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCa                           >  8:1690/1‑131 (MQ=255)                                                                                                                                           gAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCa                           <  7:1690/131‑1 (MQ=255)                                                                                                                                                  ttCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGccc  <  3:72417/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GACATGGGTAAATCCTCGCGAGCAAGCGGATGGGCGGTTATCGCTTTGGAAACGCTTGAGTGGGGGGAGGATTTAGCGGTTGTTACGCGATTTTATTTTTGTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCCCGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCACAGGCGATAA‑ACATCAATGCCCAGGACGCGCCGGTACCTGCAACAAAACCCGCAAGCCGTCACTACGGCTGTCAAAGCGCAGGCTACCGGCACACCGCTGCACGATCGCCTGCACAATCGCCAGCCC  >  NC_002947/4796480‑4796773

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: ‑