Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,220,718	A→T	8.0%	L303H (CTC→CAC)	PP_1960 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,220,718	0	A	T	8.0%	142.3 / 3.5	50	L303H (CTC→CAC)	PP_1960	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base T (2/2);  total (28/22)
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.82e-01

ATGAGGAACATCGTCGGCCACCAGTGTCGAGACTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGTTTTATTACCCTTGC  >  NC_002947/2220571‑2220865                                                                                                                                                        |                                                                                                                                                     aTGAGGAACATCGTCGGCCACCAGTGTCGAGACTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTgag                                                                                                                                                     <  7:53519/149‑1 (MQ=255)    aGGAACATCGTCGGCCACCAGTGTCGAGACTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGTGCGCCTTTGATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTgag                                                                                                                                                     <  1:106584/149‑1 (MQ=255)    aGGAACATCGTCGGCCACCAGTGTCGAGACTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGaa                                                                                                                                                   <  4:287858/148‑1 (MQ=255)          aTCGTCGGCCACCAGTGTCGAGACTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAAc                                                                                                                                                 >  1:59064/1‑144 (MQ=255)                gtccaccaGTGTCGAGTCTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAag                                                                                                                                       <  3:217540/146‑1 (MQ=255)                ggCCACCAGTGTCGAGACTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAag                                                                                                                                       <  2:59064/148‑1 (MQ=255)                       aGTGTCGAGACTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCaa                                                                                                                               >  6:75729/1‑149 (MQ=255)                             gagaCTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATGGAAGGAGGGGAATCGTgg                                                                                                                         >  4:270474/1‑149 (MQ=255)                                  tCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCCTTGCCCCTAGA‑CGCAATTTTCGTCGCGGTAACCGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGa                                                                                                                     >  3:36264/1‑148 (MQ=255)                                          aTTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCt                                                                                                            >  2:262904/1‑149 (MQ=255)                                           ttCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCtt                                                                                                           <  7:202784/149‑1 (MQ=255)                                                ggATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACt                                                                                                      >  3:231551/1‑149 (MQ=255)                                                      ccTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACtttttt                                                                                                 <  1:165703/148‑1 (MQ=255)                                                      ccTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTc                                                                                                >  3:128939/1‑149 (MQ=255)                                                       cTGAAGTTTTCGTCGAAATAAAGT‑‑GCC‑TTCATTTTTTGAGGGTTGCTTCTAGA‑GGCAATTTTGGTCGCTGTTACAGGGAGGCCTCTGCTATTGAGTATCAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCATTGTACttttttag                                                                                               <  8:22669/149‑3 (MQ=255)                                                       cTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTc                                                                                                >  8:49475/1‑148 (MQ=255)                                                        tGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGTAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGa                                                                                              >  1:340928/1‑149 (MQ=255)                                                        tGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGa                                                                                              >  5:144004/1‑149 (MQ=255)                                                        tGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGa                                                                                              >  1:204676/1‑149 (MQ=255)                                                        tGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGa                                                                                              >  3:200513/1‑149 (MQ=255)                                                           aGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGt                                                                                           <  7:49475/149‑1 (MQ=255)                                                               ttCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTcc                                                                                        <  1:12882/148‑1 (MQ=255)                                                                 cGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTc                                                                                      <  7:79473/148‑1 (MQ=255)                                                                 cGTCAAAAGATACT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACCGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTc                                                                                      <  4:36264/148‑1 (MQ=255)                                                                  gTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCa                                                                                     <  6:144004/148‑1 (MQ=255)                                                                   tCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCa                                                                                     >  4:8315/1‑147 (MQ=255)                                                                       aagaaagT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAAc                                                                               <  5:75729/149‑1 (MQ=255)                                                                             gT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTc                                                                                      >  3:79659/1‑136 (MQ=255)                                                                             gT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTc                                                                                      <  4:79659/136‑1 (MQ=255)                                                                                 gCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGa                                                                        >  1:258920/1‑148 (MQ=255)                                                                                     ttCATTTTTGCAGCGTTGCCCCTAGAGCGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGata                                                                      >  1:3693/1‑148 (MQ=255)                                                                                     ttCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGCCGTTCCTCATAGAACGGCTTGatat                                                                     >  3:216428/1‑148 (MQ=255)                                                                                     ttCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGatat                                                                     >  1:285002/1‑148 (MQ=255)                                                                                         tttttGCAGCGTTGCCCCTAGC‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAAccc                                                                >  6:318757/1‑149 (MQ=255)                                                                                                aGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTcc                                                                                        >  1:332321/1‑118 (MQ=255)                                                                                                aGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTcc                                                                                        <  2:332321/118‑1 (MQ=255)                                                                                                      gCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGctct                                                   <  4:216428/148‑1 (MQ=255)                                                                                                                       ttGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGctct                                   <  4:128939/148‑1 (MQ=255)                                                                                                                          gTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGaa                                <  4:200513/148‑1 (MQ=255)                                                                                                                           tCGCGGTAACAGGGACGCCTCTGCTATTGTGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCg                                                                                               <  8:182850/85‑1 (MQ=255)                                                                                                                           tCGCGGTAACAGGGACGCCTCTGCTATTGTGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCg                                                                                               >  7:187822/1‑85 (MQ=255)                                                                                                                           tCGCGGTAACAGGGACGCCTCTGCTATTGTGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCg                                                                                               >  7:182850/1‑85 (MQ=255)                                                                                                                           tCGCGGTAACAGGGACGCCTCTGCTATTGTGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCg                                                                                               <  8:187822/85‑1 (MQ=255)                                                                                                                           tCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAg                               >  5:298413/1‑148 (MQ=255)                                                                                                                               ggTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATaa                          >  1:318094/1‑149 (MQ=255)                                                                                                                                 tAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATc                        <  2:318094/149‑1 (MQ=255)                                                                                                                                           gCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGttt              >  6:217652/1‑149 (MQ=255)                                                                                                                                           gCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGttt              >  2:87044/1‑149 (MQ=255)                                                                                                                                            ccTCTGCTATTGACCAACAAGATCGAAGGAGGGCAATCGTGCCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGtttt             <  2:281765/149‑1 (MQ=255)                                                                                                                                             ctctGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGTTtta            >  5:258692/1‑149 (MQ=255)                                                                                                                                                  cTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGTTTTATTAc        >  2:304671/1‑148 (MQ=255)                                                                                                                                                     ttGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAAAGGCTTGATATAACCCTCTATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGTTTTATTAccc      <  1:262904/148‑1 (MQ=255)                                                                                                                                                       gagaAACAAGATCGAAGGAGGGCGATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGTTTTATTACCCTTGc  >  5:39625/1‑149 (MQ=255)                                                                                                                                                       gagaAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTACCTCATAGAACGGCTTGATATAACACTC‑ATCTAGCTCTCTTAGAACGTCGCTCTGAAGGATAAACGTCCGTGATTTATTACCCATGc  >  1:184207/1‑149 (MQ=255)                                                                                                                                                        agaAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGTTTTATTACCCTTGc  <  1:87044/148‑1 (MQ=255)                                                                                                                                                        agaAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGTTTTATTACCCTTGc  <  2:204676/148‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                     ATGAGGAACATCGTCGGCCACCAGTGTCGAGACTCATAGAAATTCACGGATTCCCTGAAGTTTTCGTCCAAAGAAAGT‑‑GCC‑TTCATTTTTGCAGCGTTGCCCCTAGA‑CGCAATTTTGGTCGCGGTAACAGGGACGCCTCTGCTATTGAGAAACAAGATCGAAGGAGGGCAATCGTGGCCGAATCGCTCCTTGTACTTTTTTCGACGTTCCTCATAGAACGGCTTGATATAACCCTC‑ATCTAGCTCTCTTAGATCGTCGCTCTGAAGGATAATCGTCCGTGTTTTATTACCCTTGC  >  NC_002947/2220571‑2220865

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