Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,856,087	T→A	P80P (CCT→CCA)	USA300HOU_2698 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,856,087	0	T	A	93.7%	203.4 / 6.2	63	P80P (CCT→CCA)	USA300HOU_2698	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base A (26/33);  total (27/36)
Fisher's exact test for biased strand distribution p-value = 6.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

TGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAATTCAAGA  >  CP000730/2855954‑2856227                                                                                                                                      |                                                                                                                                             tGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTc                                                                                                                                       >  1:225108/1‑141 (MQ=255)           tGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGAt                                                                                                                             <  1:632045/141‑1 (MQ=255)            ggTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATa                                                                                                                            <  1:288647/141‑1 (MQ=255)            ggTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATa                                                                                                                            <  1:98413/141‑1 (MQ=255)            gagccTAAAGCTATATCAATATTGAAACAAGATTTATTTCAGCACAAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTTAATTATCAGGAGATGAGTAATGCAATCAAGAGTCAAAACGTCAACAAATGATa                                                                                                                            <  2:780935/138‑1 (MQ=255)               ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                         <  1:861330/141‑1 (MQ=255)               ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                         >  2:285645/1‑141 (MQ=255)               ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                         <  1:464290/141‑1 (MQ=255)               ccTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATCTTAATACTGTAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAt                                                                                                                         <  1:275226/141‑1 (MQ=255)                 tAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGAttt                                                                                                                       <  1:641234/141‑1 (MQ=255)                  aaaGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGCTGTATCACGCCATCATGCGCCAAAACGTCCACAAATGATAGAattt                                                                                                                      >  2:250026/1‑141 (MQ=255)                               tatTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTAcag                                                                                                         <  2:338807/141‑1 (MQ=255)                                  tGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTAcagcag                                                                                                      <  2:410487/141‑1 (MQ=255)                                  tGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTAcagcag                                                                                                      >  1:179733/1‑141 (MQ=255)                                  tGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTAcagcag                                                                                                      >  2:170349/1‑141 (MQ=255)                                    gAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAgcagca                                                                                                    <  2:727479/141‑1 (MQ=255)                                     aaCAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAt                                                                                                   >  1:121239/1‑141 (MQ=255)                                      aCAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCAtt                                                                                                  >  1:160634/1‑141 (MQ=255)                                       cAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTa                                                                                                 <  1:144966/141‑1 (MQ=255)                                       cAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTa                                                                                                 <  1:124943/141‑1 (MQ=255)                                       cAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTa                                                                                                 <  1:24176/141‑1 (MQ=255)                                            tttATTTCATCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGatat                                                                                            <  1:331798/141‑1 (MQ=255)                                                      cACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAAttatta                                                                                  >  2:39924/1‑141 (MQ=255)                                                          aaCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTAAAGCAGCATTAGATTTTTATTTTTTACGtt                                                                              <  2:419832/141‑1 (MQ=255)                                                           aCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTc                                                                             <  2:192987/141‑1 (MQ=255)                                                           aCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTc                                                                             <  2:193063/141‑1 (MQ=255)                                                           aCTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTc                                                                             <  1:158688/141‑1 (MQ=255)                                                                aTTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTgg                                                                        <  1:645243/141‑1 (MQ=255)                                                                  tttAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGtt                                                                      >  2:609714/1‑141 (MQ=255)                                                                          ggAAGTTCAGTCATCATTACCTGTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGa                                                              >  2:399302/1‑141 (MQ=255)                                                                                  aGTCATCATTACCTTGTAAGTTAACAGGTGATGTGTCATGCAATCATGCGCCAAAATGTCAACAAATAATAGATTTTATAGTTGTTAAAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAatttattt                                                      <  1:250026/141‑1 (MQ=255)                                                                                    gcatcatTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTgg                                                    <  2:454405/140‑1 (MQ=255)                                                                                       tcatTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCt                                                 <  2:234754/141‑1 (MQ=255)                                                                                       tcatTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCt                                                 <  1:683865/141‑1 (MQ=255)                                                                                        catTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTg                                                <  1:849781/141‑1 (MQ=255)                                                                                          ttACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTc                                              <  2:497256/141‑1 (MQ=255)                                                                                             ccTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCt                                           <  2:285800/141‑1 (MQ=255)                                                                                             ccTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCt                                           >  2:546299/1‑141 (MQ=255)                                                                                               ttttAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTgg                                         >  1:579036/1‑141 (MQ=255)                                                                                               ttttAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTgg                                         <  2:160634/141‑1 (MQ=255)                                                                                                    aGTTAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCt                                    <  1:358457/141‑1 (MQ=255)                                                                                                      ttAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCttt                                  >  2:161747/1‑141 (MQ=255)                                                                                                       tAACAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTg                                 >  1:868178/1‑141 (MQ=255)                                                                                                        aaCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTATATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGa                                >  1:330851/1‑141 (MQ=255)                                                                                                         aCAGGAGATGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGAc                               <  2:313237/141‑1 (MQ=255)                                                                                                             gagaTGTNTCATGCAATCATGCGCCAAAACGTCAACAAATGTTAGATTTTATAGTTTTTACAGCAGCATTAGATATTGAATTATTGCGTTCACTGGTTACAACAGGATTTATTGGGTCTGTCCCTGGCCGCTTTGACatat                           <  2:651008/141‑1 (MQ=255)                                                                                                               gaTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATAc                         >  2:616747/1‑141 (MQ=255)                                                                                                               gaTGTGTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATAc                         >  1:95090/1‑141 (MQ=255)                                                                                                                aTGTGTCATGCACTCATGCGCCAATAAGTCAACACATGATAGCAATTATAGTAGTTTCAGCAGCACTAGAACTTGCATTATTCCGCTCAATGGTTACAACAGCATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACg                        <  1:736919/141‑1 (MQ=255)                                                                                                                aTGTGTCATGCAATCATGCGCCAAAACATCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGCATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACg                        >  2:177569/1‑141 (MQ=255)                                                                                                                 tgtgTCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACgg                       <  2:598428/141‑1 (MQ=255)                                                                                                                    gtCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGAcc                    >  1:27844/1‑141 (MQ=255)                                                                                                                     tCATGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCt                   >  2:78679/1‑141 (MQ=255)                                                                                                                       aTGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCa                 >  2:490474/1‑141 (MQ=255)                                                                                                                       aTGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCa                 <  1:855332/141‑1 (MQ=255)                                                                                                                       aTGCAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCa                 >  2:722100/1‑141 (MQ=255)                                                                                                                          cAATCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAat              <  2:323893/141‑1 (MQ=255)                                                                                                                            aTCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAatat            >  2:594434/1‑141 (MQ=255)                                                                                                                            aTCATGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAatat            >  1:857881/1‑141 (MQ=255)                                                                                                                                tGCGCCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAt        >  2:176392/1‑141 (MQ=255)                                                                                                                                 gcgcCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAtt       >  2:767595/1‑141 (MQ=255)                                                                                                                                 gcgcCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAAtt       >  2:49741/1‑141 (MQ=255)                                                                                                                                  cgcCAAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAATTc      <  2:364837/141‑1 (MQ=255)                                                                                                                                      tAAACGTCAACAAATGATAGATTTTATAGTTGCTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTGCCTGGACACTTTGACATTTATGGACCTCAAATTTTAATTCAAGa  <  2:225108/141‑1 (MQ=255)                                                                                                                                      |                                                                                                                                             TGCATGGCGTTGGTCCTAAAGCTATATCAATATTGGAACAAGCTTTATTTCAGCACCAACTACATTTTAAAACGGAAGTTCAGTCATCATTACCTTTTAAGTTAACAGGAGATGTGTCATGCAATCATGCGCCTAAACGTCAACAAATGATAGATTTTATAGTTGTTACAGCAGCATTAGATATTGAATTATTACGTTCACTGGTTACAACAGAATTTATTTGGTCTGTCCCTGGCCGCTTTGACATATACGGACCTCAAATATTAATTCAAGA  >  CP000730/2855954‑2856227

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