Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,849,941	G→A	intergenic (‑92/+33)	hycB ← / ← hycA	formate hydrogenlyase subunit HycB/regulator of the transcriptional regulator FhlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,849,941	0	G	A	86.1%	70.5 / 6.5	39	intergenic (‑92/+33)	hycB/hycA	formate hydrogenlyase subunit HycB/regulator of the transcriptional regulator FhlA
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (16/15);  total (22/17)
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.73e-01

GAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCA  >  NC_000913/2849700‑2850181                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 gagCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAACTGTCATTTTCGACACTCATCGACACACCCCTCCCCAAACAGACGTAACGcc                                                                                                                                                                                                                                         >  2:77318/1‑251 (MQ=255)            cAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCACTTACATAACGTTTCACCAGGCATTCCTCGGGTGATTGTCATTTTTGACGACAACATGCCGTAGAAATGTCATTTTCGACACTCCTAGACACGCCCATCCCCAACCAACCTGAACGCCTGAAAAACggg                                                                                                                                                                                                                              >  2:286384/1‑251 (MQ=255)                            caGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGCCTCTGAACAGGCGGCCTCACAAGTGTGGAAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTAAGGTGATTGTCATTTTTGACGAAACCCTGCCGATGAAATGTCATTTGCGACTCTCCACGACCCGCACATCCCCAAACAGTCATAACACCTGCATAACGGAAAAAGCCtaatcccac                                                                                                                                                                                                              >  2:363694/1‑243 (MQ=255)                                    ctctCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCCGCTCATGCTGCCgg                                                                                                                                                                                                      >  1:41980/1‑251 (MQ=255)                                       tcAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGACTCACAAGTGTGGCAGCCGACACAGAGCGTGGAGTCAGAAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGGAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACCCGCCCTCCCCCAAACAGCCGAAACGCCTGCAAACCGGGCACAGCCTCAGCTCATGATGCCGGGCt                                                                                                                                                                                                   >  2:137085/1‑251 (MQ=255)                                         gggCGCGGCATGGATTGCAGGCCGTGCTGGCGATGCGTCTCTGACCAGGCGGCCTCACAATTTTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCttt                                                                                                                                                                                                 <  1:323282/250‑1 (MQ=255)                                         aGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCttt                                                                                                                                                                                                 <  1:248153/251‑1 (MQ=255)                                          ggCGCGGCATTGATTGCAGACCGTGCTGGCGATGCGTCTCTGAACAGGCAGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTg                                                                                                                                                                                                >  1:135195/1‑251 (MQ=255)                                               ggccttgattgCAGGCCGTGCTGGCGATGCGTCCCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCt                                                                                                                                                                                           <  1:259468/247‑1 (MQ=255)                                               ggCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCt                                                                                                                                                                                           >  1:79207/1‑251 (MQ=255)                                                                        gATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCgg                                                                                                                                                                  >  1:59737/1‑251 (MQ=255)                                                                                        ggcggcCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATa                                                                                                                                                  <  1:147032/251‑1 (MQ=255)                                                                                        ggcggcCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGACCATCCCCAAACAGACGTAACGCCTGAAAAACGGGCAAAGCCTCAGCTCCTGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCTGCACATTGTTTTCAAt                                                                                                                                                   >  2:178532/1‑250 (MQ=255)                                                                                                  caAGTGTGGCAGCCTATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCCCCGGCCGATTGTTTTCAATAGCGGTGtat                                                                                                                                         >  2:39229/1‑250 (MQ=255)                                                                                                       ggtgcatcctaTACCGAGCGTGGAGTCAGCAATTCCAAAACCAATCCCCAGGCATTCCTCAGGTGCTTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCt                                                                                                                                    <  2:38069/240‑1 (MQ=255)                                                                                                          gCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCg                                                                                                                                >  1:124156/1‑251 (MQ=255)                                                                                                          gCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCg                                                                                                                                >  1:113322/1‑251 (MQ=255)                                                                                                               cgaaaCAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAAc                                                                                                                           <  1:145385/247‑1 (MQ=255)                                                                                                               cGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTAACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAAc                                                                                                                           >  1:126404/1‑251 (MQ=255)                                                                                                                 aTACAGAGCGTGGAGTCAACAATTACCAAACGATTAACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACCTTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGc                                                                                                                         >  1:87595/1‑251 (MQ=255)                                                                                                                         cGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAtt                                                                                                                 <  1:88397/251‑1 (MQ=255)                                                                                                                          gTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGAttt                                                                                                                <  1:49729/251‑1 (MQ=255)                                                                                                                                               cGATTCACCAGGCATTCCGCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGaa                                                                                           <  1:265774/251‑1 (MQ=255)                                                                                                                                                     aCCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGTTTTTTTCGAGATAGTGTTCCAGAACCTGtt                                                                                     >  1:217190/1‑251 (MQ=255)                                                                                                                                                         ggCATTCCTCAGGTGATTGTCCTTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTAATCGACACGCCAATCCCAAACCAGACGTAACGCATGCAAAAAGGGCAAATCCTCAGCTCATGCTGCAGGGCTTTGTCCCTTTACCAGTTGGCTTAAAGTCACCGGCACATTGTTTTAAATAGCGGAATATTAGTTGACGTAAACGCCAGCGTTTTTTTATGGATAGCGTTCCTGAACCTGTTcgcg                                                                                 >  2:102708/1‑251 (MQ=255)                                                                                                                                                            aTTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATc                                                                              >  1:380896/1‑251 (MQ=255)                                                                                                                                                              tCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCAGCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATcg                                                                             <  1:39593/250‑1 (MQ=255)                                                                                                                                                                 tCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCgg                                                                          >  1:326266/1‑250 (MQ=255)                                                                                                                                                                                            cATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCTATCGCGGATGTTGACCTTGCCATCAATCATGCCGt                                              <  1:290238/251‑1 (MQ=255)                                                                                                                                                                                               gCCGATGAAATGTCATTTTCGACACTAATCGACACGCCCATCCCCAAAAAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGCGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTCGCCATCAATCCTGCAGTCa                                            >  2:67818/1‑250 (MQ=255)                                                                                                                                                                                                   aTGAAATGTCATTTTCGACACCCACCTACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCACAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCCAATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCt                                       <  2:183387/251‑1 (MQ=255)                                                                                                                                                                                                     gAAATGTCATTTGCGACACTCATCGACACGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCa                                     <  1:371229/251‑1 (MQ=255)                                                                                                                                                                                                            catttcctaccatcctcgccccGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTgcg                               <  2:265596/229‑1 (MQ=255)                                                                                                                                                                                                              ttttCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGCTTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGTTGTTGACCTTGCCATAAATAATGCCGCCAATGCTCAGTTGCGCCt                            >  1:177290/1‑251 (MQ=255)                                                                                                                                                                                                                               cacGCCCATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGGTGCCCTTGCCATCAATCATGCCGCCCATGCTCAGTTCCGCCTCGGCAATCCGCTGtttt           >  2:205862/1‑250 (MQ=255)                                                                                                                                                                                                                                     cATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTc                                                                                    <  2:129296/172‑1 (MQ=255)                                                                                                                                                                                                                                     cATCCCCAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTc                                                                                    >  1:129296/1‑172 (MQ=255)                                                                                                                                                                                                                                       accccAACCAGACGGAACTCCTGCAAAACGGGCAAAGCCTCAGCTCACGCTGCCGGGATTTGGCGCTTTACCATTTGGCTTAAATTCATCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTGTCGAGATAGTGTTCCAGAATCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATACGCTGTTTGTTTTCGcc   <  2:87595/250‑1 (MQ=255)                                                                                                                                                                                                                                        ccccAAACAGACGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCa  <  1:77318/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 GAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATACAGAGCGTGGAGTCAGCAATTACAAAACGATTCACCAGGCATTCCTCAGGTGATTGTCATTTTTGACGAAAACATGCCGATGAAATGTCATTTTCGACACTCATCGACACGCCCATCCCCAAACAGGCGTAACGCCTGCAAAACGGGCAAAGCCTCAGCTCATGCTGCCGGGCTTTGTCCCTTTACCAGTTGGCTTAAATTCACCGGCACATTGTTTTCAATAGCGGTGTATAAGCTGTCGTAAACGCCAGCGATTTTTTCGAGATAGTGTTCCAGAACCTGTTCGCGATCGCGGATGTTGACCTTGCCATCAATCATGCCGTCAATGCTCAGTTGCGCCTGCGCAATCCGCTGTTTGTCTTCGCCA  >  NC_000913/2849700‑2850181

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