Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,706,587	T→A	12.6%	L202Q (CTG→CAG)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,706,587	0	T	A	12.6%	153.4 / 5.5	59	L202Q (CTG→CAG)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/23);  new base A (5/3);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTC  >  NC_002947/5706441‑5706732                                                                                                                                                    |                                                                                                                                                   gCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTgg                                                                                                                                                  <  5:147360/149‑1 (MQ=255)        gaAGACAGCCAAGTCTCGCGCATGCTCTGGCTTGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGa                                                                                                                                           <  3:69069/149‑1 (MQ=255)        gaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGATGTc                                                                                                                                             >  8:112004/1‑147 (MQ=255)         aaGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAg                                                                                                                                          <  2:199988/149‑1 (MQ=255)           gACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                        <  4:285298/149‑1 (MQ=255)           gACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                        <  5:7268/149‑1 (MQ=255)                 cAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATAc                                                                                                                                   >  4:103383/1‑148 (MQ=255)                  aaGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTc                                                                                                                                             <  7:41199/137‑1 (MQ=255)                  aaGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTc                                                                                                                                             >  8:41199/1‑137 (MQ=255)                     tctcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCAGGCTGTc                                                                                                                                             >  3:132548/1‑134 (MQ=255)                     tctcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCAGGCTGTc                                                                                                                                             <  4:132548/134‑1 (MQ=255)                         gcgcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGAAGCAAAAGAAGCAATGTACGCCCTGCTCGAC‑GGCGTACACGACAGCCACGCTGTCGAGGCGATACAAGTCACg                                                                                                                           >  5:270892/1‑148 (MQ=255)                                 tctGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCAGGCTGTCgt                                                                                                                                           >  1:295076/1‑123 (MQ=255)                                 tctGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCAGGCTGTCgt                                                                                                                                           <  2:295076/124‑2 (MQ=255)                                         gAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgg                                                                                                          >  3:158476/1‑149 (MQ=255)                                           gCACACGATCATGGCTGCATTTCTCGAGAACCGGCCTGTCCTAGTGCTTGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACTCGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgggg                                                                                                        <  6:270892/149‑1 (MQ=255)                                           gCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgggg                                                                                                        >  6:29273/1‑149 (MQ=255)                                                  aTCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCacacac                                                                                                 >  8:13265/1‑149 (MQ=255)                                                  aTCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCacacac                                                                                                 <  1:173270/149‑1 (MQ=255)                                                    cAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCg                                                                                               >  8:65082/1‑149 (MQ=255)                                                            cATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTa                                                                                       >  4:79159/1‑149 (MQ=255)                                                              ttGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTAcc                                                                                     <  1:58816/149‑1 (MQ=255)                                                                    gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATc                                                                                >  6:315462/1‑148 (MQ=255)                                                                    gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATc                                                                                <  8:26518/148‑1 (MQ=255)                                                                           tccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                         >  4:113654/1‑148 (MQ=255)                                                                            ccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtca                                                                       >  2:165002/1‑149 (MQ=255)                                                                             cctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGAt                                                                                 >  5:279933/1‑138 (MQ=255)                                                                             cctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGAt                                                                                 <  6:279933/138‑1 (MQ=255)                                                                                         cTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTgg                                                         <  3:79159/149‑1 (MQ=255)                                                                                              aGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGgcaacc                                                    <  4:150584/149‑1 (MQ=255)                                                                                                     gCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCt                                             >  4:22506/1‑149 (MQ=255)                                                                                                     gCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGAAACCGCAACCt                                             >  5:108917/1‑149 (MQ=255)                                                                                                      cAAAAGCTGCAATGTACGCCCTGCTCGACGGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCACCCGCAACCt                                             >  1:43848/1‑149 (MQ=255)                                                                                                       aaaaGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCacacac                                                                                                 <  8:180514/96‑1 (MQ=255)                                                                                                       aaaaGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCacacac                                                                                                 >  7:180514/1‑96 (MQ=255)                                                                                                           gCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCt                                       <  2:43848/149‑1 (MQ=255)                                                                                                           gCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCt                                       <  5:301436/149‑1 (MQ=255)                                                                                                             tGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACATGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTAATC‑ACCACAAGTGCCAACCGCAACCTCGTTCtg                                      >  6:212192/1‑148 (MQ=255)                                                                                                               cAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCtgcaa                                   <  7:112004/149‑1 (MQ=255)                                                                                                               cAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCtgcaa                                   >  5:47853/1‑149 (MQ=255)                                                                                                               cAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCtgcaa                                   >  3:305432/1‑149 (MQ=255)                                                                                                                aaTGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCtgcaat                                  >  3:138557/1‑149 (MQ=255)                                                                                                                 aTGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCtgcaatg                                 <  3:22506/149‑1 (MQ=255)                                                                                                                         ccTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCt                         >  2:282596/1‑149 (MQ=255)                                                                                                                           tGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGc                       >  5:36176/1‑149 (MQ=255)                                                                                                                             cTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcgg                     <  1:266052/149‑1 (MQ=255)                                                                                                                              tCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggt                    >  8:156817/1‑149 (MQ=255)                                                                                                                                 aC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttcg                 >  6:33656/1‑149 (MQ=255)                                                                                                                                    ggggGACACGACAGCCTGGCTGTCGTGGCGATACAACTCACGCTCGGAACGGATGGCGGGGCACAAACCGTTACCCTACCGGATCTATTGCTCATCAACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttcgg                <  1:244183/149‑1 (MQ=255)                                                                                                                                       gTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGgcaacc                                                    >  3:71116/1‑109 (MQ=255)                                                                                                                                       gTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGgcaacc                                                    <  4:71116/109‑1 (MQ=255)                                                                                                                                         cctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCt                                       >  2:197498/4‑120 (MQ=255)                                                                                                                                         cctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCt                                       <  1:197498/117‑1 (MQ=255)                                                                                                                                          ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTTCGGTTCGGTtcgtag         >  5:323593/3‑147 (MQ=255)                                                                                                                                           acGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGc        >  3:192962/1‑149 (MQ=255)                                                                                                                                              aCAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCtgcaat                                  >  6:29197/1‑120 (MQ=255)                                                                                                                                              aCAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCtgcaat                                  <  5:29197/120‑1 (MQ=255)                                                                                                                                              aCAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCAt     <  6:108917/149‑1 (MQ=255)                                                                                                                                                  ccTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTc  >  5:188644/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATC‑ACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTC  >  NC_002947/5706441‑5706732

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