Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	714,967	T→A	10.7%	intergenic (+33/+44)	ispH → / ← PP_0607	4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	714,967	0	T	A	10.7%	154.5 / 8.3	56	intergenic (+33/+44)	ispH/PP_0607	4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase/type IV pili biogenesis protein FimT
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base A (3/3);  total (29/27)
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

GATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACA  >  NC_002947/714823‑715111                                                                                                                                                  |                                                                                                                                                  gATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTg                                                                                                                                                 >  7:59196/1‑147 (MQ=255)   tCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGc                                                                                                                                             <  2:242209/149‑1 (MQ=255)    cGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGcc                                                                                                                                            >  8:104843/1‑149 (MQ=255)     gAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGccc                                                                                                                                           <  8:59196/149‑1 (MQ=255)       aCAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCtt                                                                                                                                         >  1:117756/1‑149 (MQ=255)           cTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTgcgc                                                                                                                                     >  1:119332/1‑149 (MQ=255)            tCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTgcgc                                                                                                                                     >  3:193267/1‑148 (MQ=255)              aaGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTgcgc                                                                                                                                     <  4:193267/146‑1 (MQ=255)                ggCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGacac                                                                                                                                <  2:119332/149‑1 (MQ=255)                ggCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGaca                                                                                                                                 >  2:70017/1‑148 (MQ=255)                ggCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAAGGGCTGC‑AAAGCAGCCCCTgcattggcg                                                                                                                                <  2:112808/149‑8 (MQ=255)                   cTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGacac                                                                                                                                >  4:146571/1‑146 (MQ=255)                      gggCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGa                                                                                                                            <  3:30833/147‑1 (MQ=255)                           ccACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGc                                                                                                                                                <  4:35236/122‑1 (MQ=255)                           ccACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGc                                                                                                                                                >  3:35236/1‑122 (MQ=255)                                       aaGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACcg                                                                                                          <  2:89311/148‑1 (MQ=255)                                          agCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACcgctgc                                                                                                      >  3:67982/1‑146 (MQ=255)                                             tGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTag                                                                                                   >  3:57057/1‑148 (MQ=255)                                                aCGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTagcagc                                                                                               <  1:238516/149‑1 (MQ=255)                                                aCGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTagcagc                                                                                               <  5:205502/149‑1 (MQ=255)                                                  gggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGCTTTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCg                                                                                              <  1:67307/149‑1 (MQ=255)                                                   ggCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGt                                                                                            <  4:57057/149‑1 (MQ=255)                                                      gcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCAt                                                                                         <  5:147215/149‑1 (MQ=255)                                                      gcgAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCAt                                                                                         <  8:2071/149‑1 (MQ=255)                                                         aaGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGAAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCt                                                                                                     >  7:53178/1‑134 (MQ=255)                                                          agagAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGc                                                                                      <  4:67982/148‑1 (MQ=255)                                                          agagAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCt                                                                                     <  7:104843/149‑1 (MQ=255)                                                             gaACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGCCACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAg                                                                                  >  2:206074/1‑149 (MQ=255)                                                                   cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcaggagc                                                                                >  5:211172/1‑139 (MQ=39)                                                                   cACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGCTGCTAGC‑GCGGTCATTGcaggagc                                                                                <  6:211172/145‑7 (MQ=39)                                                                      cTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCg                                                                          >  1:156890/1‑148 (MQ=255)                                                                         ttCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTaa                                                                       <  3:146571/148‑1 (MQ=255)                                                                                   aaaGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAc                                                             >  4:74792/1‑148 (MQ=255)                                                                                                ttCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATg                                               >  5:88722/1‑149 (MQ=255)                                                                                                   gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                            >  4:21930/1‑133 (MQ=255)                                                                                                   gctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                            <  3:21930/133‑1 (MQ=255)                                                                                                    ctcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAg                                           <  1:63482/149‑1 (MQ=255)                                                                                                    ctcgctGATCTGACGGGCAGTACCGCCGACAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAg                                           <  3:4910/149‑1 (MQ=255)                                                                                                     tcgctGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGc                                          >  8:217257/1‑149 (MQ=255)                                                                                                            tCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCa                                            >  4:128864/1‑140 (MQ=255)                                                                                                            tCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCa                                            <  3:128864/140‑1 (MQ=255)                                                                                                                  gggCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATc                             <  3:74792/149‑1 (MQ=255)                                                                                                                    gCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCgg                           >  1:226800/1‑149 (MQ=255)                                                                                                                             ccgccAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                            >  2:68789/1‑107 (MQ=255)                                                                                                                             ccgccAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGAct                                                            <  1:68789/107‑1 (MQ=255)                                                                                                                               tctcgcgCAAGGGCTGC‑AAAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagcag                 <  1:114268/142‑1 (MQ=255)                                                                                                                                  aaTTCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacc                                       >  5:35087/1‑123 (MQ=255)                                                                                                                                  aaTTCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGcacc                                       <  6:35087/123‑1 (MQ=255)                                                                                                                                   gcgCAAGGGCTGC‑AGAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAgcggcg            >  7:228188/3‑149 (MQ=255)                                                                                                                                   gcgCAAGGGCTGC‑AAAGCAGCCCCTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGagcag                 >  2:114268/3‑144 (MQ=255)                                                                                                                                     gCAAGGGCTGC‑AAAGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGgcgc           <  8:64877/148‑1 (MQ=255)                                                                                                                                     gCAAGGGCTGC‑AAAGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGcgcg          >  7:64877/1‑149 (MQ=255)                                                                                                                                        ggggCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTc       >  7:143172/1‑149 (MQ=255)                                                                                                                                        agggCTGC‑AAAGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTc       >  3:51162/2‑149 (MQ=255)                                                                                                                                             tgc‑tTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAg                                           >  5:171790/1‑108 (MQ=255)                                                                                                                                             tgc‑tTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGcaca  >  3:215866/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GATCGAACAGCTCAAGGCCTGGGGCGCCACTGGCGCTGAAGAGCTGGACGGGCGCGAAGAGAACATCACCTTTTCAATGCCCAAAGAGCTGCGGGTTCGCTCGCTGATCTGACGGGCAGTACCGCCGCCAATGCAGGGGCTGC‑TTTGCAGCCCTTGCGCGACACAAGACCGCTCCTGCAATGACCGC‑GCTAGCAGCGGTCATTGCTCAGGTCCTCCGTAAGCAACCTGACTCTGCCGCTTGATGCCAGCACCACCTGATATCGGCTGCGAGCAGCGGCGCGTTCGCACA  >  NC_002947/714823‑715111

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