Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,993,531	T→C	R793R (CGT→CGC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,531	0	T	C	94.8%	235.9 / ‑3.3	77	R793R (CGT→CGC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (33/40);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

ACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGC  >  NC_000913/3993290‑3993772                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTAc                                                                                                                                                                                                                                           <  1:699006/251‑1 (MQ=255) aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTAc                                                                                                                                                                                                                                           <  2:323678/251‑1 (MQ=255)        gcgcgggcggTCAGGGCAGGTTGAAACCAAGCACGGCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCCCCCAGTTCTCTGACGAATAAACGCAAGACGAGAATTGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACGGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGGCGCGTCTATTCGTCATCGCTTGTCCGCTTTACCTACGGcgca                                                                                                                                                                                                                                     <  2:78115/241‑4 (MQ=255)                     gTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑Ctt                                                                                                                                                                                                                               <  1:106153/242‑1 (MQ=255)                       gCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACGCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTc                                                                                                                                                                                                                     <  1:639693/250‑1 (MQ=255)                        cAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAACGAAGGGATCATCCCGTTCTTTTTCGAAGAAACGCAAAACGAGAATGGATTTTATATCTACATTCTCGGCGAAAACAAACGGGGTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATAACCGTTTTACCTAACGCCCCAG‑CTTCATCACCttaa                                                                                                                                                                                                                    >  1:80636/1‑248 (MQ=255)                              gAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGGTCATCCAGTTCTTTTTTGAAGGAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATAACACCCCTGCGAAGGCAACAAAGAGGGGCTGGTACGGGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCCTTACCTACGGCTCAAG‑CTTCCTCAACTTCAACCTGcc                                                                                                                                                                                                             >  1:174993/1‑251 (MQ=255)                                 aCCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCa                                                                                                                                                                                                          <  1:311113/251‑1 (MQ=255)                                             aaaTTACCGGCGGTGGGGGACGGCTTTGCCAGCGAAAGGATCATCCAGGTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACCTTCTCGACGAAAGCAACCGGGTTGAGGGATATCACCACTGCGAAGGCAGCAAAAAGGAGCTGGTACGTGACGGCAGTCGCTTCTACTCGTCATCCCATGACCGCTTTTCCTACCGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATc                                                                                                                                                                                                  >  1:294077/1‑247 (MQ=255)                                               aTTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATtg                                                                                                                                                                                            <  1:348376/251‑1 (MQ=255)                                                   ccGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGTTCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACATAAGGCACAAG‑CTTCATCAACTTCAACCTGCAGCAGTTCTATCAGATTGTGa                                                                                                                                                                                         >  2:663997/1‑250 (MQ=255)                                                         gtggtgGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCAACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTg                                                                                                                                                                                   >  1:59842/1‑250 (MQ=255)                                                            gtgtaCGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATgg                                                                                                                                                                               <  2:259092/247‑1 (MQ=255)                                                            gtggACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGAAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTCTCAGATTGTGAAGGTTGATgg                                                                                                                                                                               >  1:31030/1‑251 (MQ=255)                                                                aCGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGt                                                                                                                                                                           <  2:689143/251‑1 (MQ=255)                                                                  ggCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGa                                                                                                                                                                         <  1:681482/251‑1 (MQ=255)                                                                   gCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGAAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGa                                                                                                                                                                         >  1:673368/1‑250 (MQ=255)                                                                        gCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAgg                                                                                                                                                                    <  1:150867/250‑1 (MQ=255)                                                                             cGAAGGGATCAACCAGTTCTTTTTCGAAGAAACGCAAGACGATAATGGCTTTAATATCTACATTCGCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGGCATCGCATGACCGTTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAttc                                                                                                                                                              <  2:346504/251‑1 (MQ=255)                                                                                  ggATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACATACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATTATATTGTGAAGGTTGATGGTCGTGAACAGGTGAttccgttc                                                                                                                                                         >  2:494302/1‑251 (MQ=255)                                                                                   gATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTGTATCAACACTGCGTGGGCAGCAAAGAGGAGCAGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACATACGGCGCAAG‑CTTCATCAAATTCAACCt                                                                                                                                                                                                                >  2:504105/1‑195 (MQ=255)                                                                                   gATCATCCAGTTCTGTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTGTATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACAGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCt                                                                                                                                                                                                                <  1:504105/195‑1 (MQ=255)                                                                                           aGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaa                                                                                                                                                 >  1:127815/1‑250 (MQ=255)                                                                                            gTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAAt                                                                                                                                               >  1:654921/1‑251 (MQ=255)                                                                                                   ttCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCg                                                                                                                                         <  1:421768/250‑1 (MQ=255)                                                                                                     cGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCAACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACGTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGt                                                                                                                                       <  2:59842/250‑1 (MQ=255)                                                                                                       aagaaACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGACGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAAc                                                                                                                                    >  1:372684/1‑251 (MQ=255)                                                                                                         gaaACGCAAGACTAGAATGGCTTTAATATCTACATTCGCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACAt                                                                                                                                  <  1:284811/251‑1 (MQ=255)                                                                                                             cGCAAGACAGGAATGGCCTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGTGGTATATCCCCCATGCGAATGCAGCAAAGAGGAGCTGTGACGTGACGTCAGTCGCTTCTACTCGTCATCTCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAATGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATgctg                                                                                                                              <  2:338479/251‑3 (MQ=255)                                                                                                             cGCAAGAAGAGAATGGCTATAAGATTTACATTCTCGATGAAAGCAACCGGGTTGAGTTATTCCACCCCGGCGAAGGCAGCAAAGGGGAGCGGGTACGTGACGTCAGTCGCTTCTACTCGGCACCGCAGGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTGCTAGCAGATTGTGAAGGTTGATGGTCGTGAACAGGGGAGTCCGTTCCGCACAAAATCTATCGGTAACATgccg                                                                                                                              <  2:420243/251‑1 (MQ=255)                                                                                                                      agaATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAAt                                                                                                                     >  1:593186/1‑251 (MQ=255)                                                                                                                       gaATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATc                                                                                                                    >  1:527296/1‑251 (MQ=255)                                                                                                                        aaTGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTGCAACCTGCCGCAGTTCTATCATATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATc                                                                                                                    <  1:40964/250‑1 (MQ=255)                                                                                                                         aTGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCCCCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGACGCCTGCCAATc                                                                                                                    >  1:694483/1‑249 (MQ=255)                                                                                                                               ttAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGACGCCTGCCAATCAGGAGCAc                                                                                                            >  1:440826/1‑251 (MQ=255)                                                                                                                                tAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCAc                                                                                                            <  1:621699/250‑1 (MQ=255)                                                                                                                                tAATATCTACAGTCTCGACGAAAGCAACCTGGTTGAGGTATATCACCCCTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAG‑CTTAATCAACTTCATCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCAc                                                                                                            <  2:519968/250‑1 (MQ=255)                                                                                                                                       tACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgcc                                                                                                    >  1:321242/1‑251 (MQ=255)                                                                                                                                        aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCTCTGCGAAGGCAGCAAAGAGGCGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCGGCCAATCAGGATCACGATAc                                                                                                       <  1:568072/247‑1 (MQ=255)                                                                                                                                        aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgc                                                                                                     <  1:687677/249‑1 (MQ=255)                                                                                                                                        aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTGCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgcc                                                                                                    <  2:537274/250‑1 (MQ=255)                                                                                                                                        aCATTCGCTACGAAAGCAACCGGGTTGAGGTATATCACCACGTCTAAGTCAGCAAAGAGGAGCTGGTACGTGACTGCAGTCGCTTCTACTCGTCATCGCATGACCGCTGTACCTACGGCTCAAG‑CTGCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTTATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAAAATGCCGCCTGCCAATCAGGATCACGATat                                                                                                       <  2:315730/247‑2 (MQ=255)                                                                                                                                             ctcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTAtt                                                                                              >  1:50534/1‑251 (MQ=255)                                                                                                                                               cgacgaAAGCACCCGGTTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTTGTACGTGACGTCAGTCGCTTCTACTCTTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTTAAGGTTGATAGTCGt                                                                                                                                                                           <  2:636504/172‑1 (MQ=255)                                                                                                                                               cgacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTCCGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGt                                                                                                                                                                           >  1:636504/1‑172 (MQ=255)                                                                                                                                                gacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACa                                                                                           <  1:263573/251‑1 (MQ=255)                                                                                                                                                 acgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGca                                                                                                                                                     <  1:698434/192‑1 (MQ=255)                                                                                                                                                 acgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGca                                                                                                                                                     >  2:698434/1‑192 (MQ=255)                                                                                                                                                   gaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCCCGATACGCCGCTATTACAGc                                                                                         >  2:592873/1‑250 (MQ=255)                                                                                                                                                        cAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAtt                                                                                   >  1:447012/1‑251 (MQ=255)                                                                                                                                                         aaCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAttt                                                                                  <  1:433117/251‑1 (MQ=255)                                                                                                                                                               gTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtg                                                                            >  1:191065/1‑251 (MQ=255)                                                                                                                                                                  gAGGTATATCACCACTGCGAAGGCAGCAAAGAGTAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatg                                                                         <  1:667349/251‑1 (MQ=255)                                                                                                                                                                  gAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgat                                                                          <  1:419244/250‑1 (MQ=255)                                                                                                                                                                   aGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatga                                                                        <  1:494302/251‑1 (MQ=255)                                                                                                                                                                    ggTATATCACCACTGCGAAGGCAGCACAGAGGAGCTGGTACGCGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCACCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCACTATTTTTCGtgatga                                                                        >  1:264343/1‑250 (MQ=255)                                                                                                                                                                    ggTATAACACCCCTGCGAAGGCAGCAAAGGGGAGGTGGTACGTGACGTCAGTCTCTTCTACGCGTCATCGCATGACCGCTTTATCATCGGC‑GAAGTCTTTATCAACTTCAACCTGCCGCAGTTCTATCAGATTGGGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGaa                                                                       <  2:581615/251‑1 (MQ=255)                                                                                                                                                                       atatCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACg                                                                     <  1:663997/250‑1 (MQ=255)                                                                                                                                                                         atCACCACTGCGAAGGCAGCAAAGATGAGCTGGTACTTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGc                                                                  <  1:573882/251‑1 (MQ=255)                                                                                                                                                                           caccacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTTACGGCAGTCGCTTCTACTCGGCATTGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATTGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCg                                                                <  2:384948/251‑1 (MQ=255)                                                                                                                                                                                    gAAGGCAGCAAAGAGGAGCTGGTACGTTACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAg                                                       <  1:467943/251‑1 (MQ=255)                                                                                                                                                                                    gAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCt                                                                                                                                             >  1:292469/1‑165 (MQ=255)                                                                                                                                                                                    gAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCt                                                                                                                                             <  2:292469/165‑1 (MQ=255)                                                                                                                                                                                    gAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCg                                                         >  1:4778/1‑249 (MQ=255)                                                                                                                                                                                          agcaAAGAGGAGCTGGTACGGGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGGTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTGTTCGTGATGAACGTGCCGGAAAGCGAGGCTTAt                                                 <  2:447012/251‑1 (MQ=255)                                                                                                                                                                                                 aggagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCAt                                          <  1:388868/251‑1 (MQ=255)                                                                                                                                                                                                    cgCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTTATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAACGAGGCTTATCCGGCaata                                        >  2:581028/2‑247 (MQ=255)                                                                                                                                                                                                         gTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCtt                                  >  1:11812/1‑251 (MQ=255)                                                                                                                                                                                                                      gTCGCTTCTACTCGTCACCTCATGACCGCTTTACCTACGGCAAAAG‑CGTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGTTGATTCCGTTCCGCACAAAATCTATCGGTAACAGGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactgac                     <  2:191065/251‑1 (MQ=255)                                                                                                                                                                                                                           ttCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTGTCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGCTGAACGTGCCGGACAGCGAAGCTTATCCGGCATGCACCCTTAGCGGAACCTGCCTGttt                >  2:683391/1‑251 (MQ=255)                                                                                                                                                                                                                            tctcctcttcaacgcattcccGTTTTACCTCCGGCTCAAG‑CTTGTTCAACTTCAACCGGCCGCATTTCTATCAGATTGTGAAGTTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAactgactggtg                <  2:628486/231‑4 (MQ=255)                                                                                                                                                                                                                            tCTACTCTTCATCCCATGACCGCTTAACCTACGGCGCAAG‑CTTCATCAACTTCATCCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGtt                 <  2:485826/249‑1 (MQ=255)                                                                                                                                                                                                                            tCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATAtttt                                                                                 >  1:205411/1‑185 (MQ=255)                                                                                                                                                                                                                                cTCGTCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGCTGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAACGCGCGGCTTATCCGGCATGCAATCTTAGCGGAACCTGACTGTTTCAccc           >  2:706316/1‑251 (MQ=255)                                                                                                                                                                                                                                    tCATCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCct       >  2:362718/1‑251 (MQ=255)                                                                                                                                                                                                                                      aTCGCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCctgc     >  2:217209/1‑251 (MQ=255)                                                                                                                                                                                                                                         gCATGACCGCTTTACCTACGGCTCAAG‑CTTCATCCACTTCAACCTGCCGCAGATCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCAGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCtgctgc  >  1:634248/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   ACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAACGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAG‑CTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGC  >  NC_000913/3993290‑3993772

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