Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,168,350	Δ1 bp	coding (483/633 nt)	ycfQ ←	repressor for bhsA(ycfR)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,168,348	0	G	.	98.6%	274.0 / ‑2.6	69	coding (485/633 nt)	ycfQ	repressor for bhsA(ycfR)
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (34/34);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.11e-01

CGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGCGGGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATCATGAAG  >  NC_000913/1168233‑1168482                                                                                                                    |                                                                                                                                       cGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAGTTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGtt                                                                                                                >  2:999219/1‑139 (MQ=255)  gTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGCGGGATTTCCCCGCGCGCTTGTCGtt                                                                                                                <  1:1375799/139‑1 (MQ=255)              tGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAActg                                                                                                   >  2:195975/1‑139 (MQ=255)                     gTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCaaa                                                                                            <  2:259645/139‑1 (MQ=255)                      tttttCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCa                                                                                              <  2:580399/136‑1 (MQ=255)                      tttttCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCa                                                                                              >  1:580399/1‑136 (MQ=255)                        tttCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTg                                                                                         <  2:807747/139‑1 (MQ=255)                                gATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATaaa                                                                                                        <  2:267400/116‑1 (MQ=255)                                gATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATaaa                                                                                                        >  1:267400/1‑116 (MQ=255)                                gATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTg                                                                                 >  2:1041702/1‑139 (MQ=255)                                      ccTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCGCTTGCATCGc                                                                           <  1:1269097/139‑1 (MQ=255)                                         tCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTCCAAAGTGCGCTCTTGc                                                                                >  1:1266341/1‑131 (MQ=255)                                         tCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTCCAAAGTGCGCTCTTGc                                                                                <  2:1266341/131‑1 (MQ=255)                                           gcgTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTccc                                                                                                                               <  1:1478379/82‑1 (MQ=255)                                           gcgTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTccc                                                                                                                               >  2:1478379/1‑82 (MQ=255)                                             gTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGt                                                                    <  1:1041702/139‑1 (MQ=255)                                              tGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTg                                                                                         >  1:201008/1‑117 (MQ=255)                                              tGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTg                                                                                         <  2:201008/117‑1 (MQ=255)                                                         aCATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTa                                                        <  2:838888/139‑1 (MQ=255)                                                                 tGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGat                                                <  1:182648/139‑1 (MQ=255)                                                                  gAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                               <  1:999219/139‑1 (MQ=255)                                                                  gAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                               <  1:385187/139‑1 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATaaaaa                                                                                                      >  2:618293/1‑78 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATaaaaa                                                                                                      <  1:618293/78‑1 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGTGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                         >  1:407005/1‑139 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCTTGGCGTGATTTCAACGTATTGGCGATATCTCCt                                         >  1:489713/1‑139 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCg                                                                     <  1:1236814/111‑1 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCg                                                                     >  2:1236814/1‑111 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGat                                                <  1:1590702/132‑1 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGat                                                >  2:1590702/1‑132 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                         >  2:629323/1‑139 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                         >  1:1619703/1‑139 (MQ=255)                                                                        taCAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCt                                         >  1:51634/1‑139 (MQ=255)                                                                           aaTTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                               >  1:1070453/1‑130 (MQ=255)                                                                           aaTTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                               <  2:1070453/130‑1 (MQ=255)                                                                            aTTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAAt                                     <  1:629323/139‑1 (MQ=255)                                                                            aTTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAAt                                     <  2:489713/139‑1 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGctc                                                                                    <  1:19165/91‑1 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGctc                                                                                    >  2:19165/1‑91 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGc                                                                                >  1:1346008/1‑95 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGc                                                                                <  2:1346008/95‑1 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTg                                                     >  1:201513/1‑122 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTg                                                     <  2:201513/122‑1 (MQ=255)                                                                             ttAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATc                                    >  1:1383608/1‑139 (MQ=255)                                                                               aaGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCcg                                                                                                                             <  2:413929/48‑1 (MQ=255)                                                                               aaGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCcg                                                                                                                             >  1:413929/1‑48 (MQ=255)                                                                                   aaTTCTGCCAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGga                              <  1:195975/139‑1 (MQ=255)                                                                                          ccAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTg                                        >  2:146070/1‑122 (MQ=255)                                                                                          ccAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTg                                        <  1:146070/122‑1 (MQ=255)                                                                                          ccAGATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGggtg                       >  1:1740938/1‑139 (MQ=255)                                                                                             gATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGTGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGgtggtg                    <  2:407005/139‑1 (MQ=255)                                                                                             gATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                               <  1:514446/112‑1 (MQ=255)                                                                                             gATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                               >  2:514446/1‑112 (MQ=255)                                                                                             gATGTGTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGgtggtg                    <  2:1740938/139‑1 (MQ=255)                                                                                               tgtgTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                               >  2:1432287/1‑110 (MQ=255)                                                                                               tgtgTCACGTCACAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGata                                               <  1:1432287/110‑1 (MQ=255)                                                                                                        tcacAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAActgct                                                                                                 <  2:1698146/51‑1 (MQ=255)                                                                                                        tcacAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAActgct                                                                                                 >  1:1698146/1‑51 (MQ=255)                                                                                                        tcacAGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGAt         >  1:516132/1‑139 (MQ=255)                                                                                                         cacaGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTg                                        >  2:1231297/1‑107 (MQ=255)                                                                                                         cacaGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTg                                        <  1:1231297/107‑1 (MQ=255)                                                                                                         cacaGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATc        >  2:410196/1‑139 (MQ=255)                                                                                                         cacaGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATc        >  2:1170829/1‑139 (MQ=255)                                                                                                           caGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCg                                                                     <  1:1497849/76‑1 (MQ=255)                                                                                                           caGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCg                                                                     >  2:1497849/1‑76 (MQ=255)                                                                                                            aGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAAc                                                           <  2:1691737/85‑1 (MQ=255)                                                                                                            aGTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAAc                                                           >  1:1691737/1‑85 (MQ=255)                                                                                                             gTGGGGC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATCATGa    >  1:1698764/1‑139 (MQ=255)                                                                                                               ggggC‑GGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATCATGAag  <  1:1864312/139‑1 (MQ=255)                                                                                                                    |                                                                                                                                       CGTTCCGGCAATCTGCATCAGTTTTTCCAGCGATGCACCTTCGCGTGCGCTGATCGACATCCCCTGAATAATACAATTAAGGAATTCTGCCAGATGTGTCACGTCACAGTGGGGCGGGATTTCCCCGCGCGCTTGTCGTTGACATAAAAACTGCTGCAAAGTGCGCTCTTGCATCGCATGGCGTGATTTCAACGTATTGGCGATATCTCCTGAATCTGGGGAGAGGGTGGTGCAGTTGTTGATCATGAAG  >  NC_000913/1168233‑1168482

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