Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,293,781	2 bp→AT	coding (239‑240/597 nt)	nrfG →	putative formate‑dependent nitrite reductase complex subunit NrfG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,293,781	0	G	A	100.0%	129.3 / NA	38	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/23);  total (15/23)
*	NC_000913	4,293,782	0	C	T	94.7%	120.2 / ‑5.8	38	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (15/21);  minor base G (0/2);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTC  >  NC_000913/4293533‑4293997                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                        gggCGAAGAGATGAAACAGCCCAAAATACCGGTGATAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGACCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGCTCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGATTCAGGCATTGCAGGATAAAATCCGCGCTACTCCACAAAACAGCGAACAGTGGGCGATACTgggtg                                                                                                                                                                                                                        >  2:26349/1‑249 (MQ=255)                    cccAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTGCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCaaaa                                                                                                                                                                                                     <  1:253625/251‑1 (MQ=255)                        aaaTACCGGTGAAAATTCTTACAACACTCCCGATATTGATTGTATTTGTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCTCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCCGGCATTGCAGGATAAAATCCGCGCTATTCCACAAAACAGCGAACAGTGGTCGTTACTGGAGGAGTACTATCTGTTGCAAAACg                                                                                                                                                                                                   <  2:304906/249‑1 (MQ=255)                            aCCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTAc                                                                                                                                                                                              <  1:26349/250‑1 (MQ=255)                                gtaaaagtGCTGACAACACTCACGATATTGCTGGTATTTCTCTGAGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGGTTGCCAGCCAGCAAAACCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTCATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGAACTATCTGTGGCAAAACGATTACatca                                                                                                                                                                                          >  2:12605/8‑247 (MQ=255)                                   aaaaTGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTc                                                                                                                                                                                      <  2:215714/251‑1 (MQ=255)                                      aTGCTTACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGtt                                                                                                                                                                                   <  2:221656/251‑1 (MQ=255)                                        gCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCGGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGc                                                                                                                                                                                 <  1:78448/251‑1 (MQ=255)                                          tGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTg                                                                                                                                                                               >  1:270074/1‑251 (MQ=255)                                                     aCGATATTGATGGTAGTTCTTTGTGTCGGCAGGTATCTGTGAAGCCCAAAATGGCAGGCGGTACGGGCGGAGTATCAGCGTCATCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAAGCCGCGCTAATCCACAGAACAGCTACCATTTGGCGTTACTGGATGAGTCCTATCTGTGGCAAATCTATTACAGCAGTTCTTAGCTGGCGTACCGTCa                                                                                                                                                                    <  2:309861/251‑1 (MQ=255)                                                         tatTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGt                                                                                                                                                                                    <  1:277892/231‑1 (MQ=255)                                                         tatTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGACAGCCAGCAAACCCCGGCAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGt                                                                                                                                                                                    >  2:277892/1‑231 (MQ=255)                                                              aTGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCa                                                                                                                                                           <  1:42043/251‑1 (MQ=255)                                                              aTGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCa                                                                                                                                                           <  1:298857/251‑1 (MQ=255)                                                                   aTTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAACATGGCAGACGGTACGTACGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCCAACCCCGGAAGCGCAGCTTCAGACATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCCGGCGTTGCAACTAc                                                                                                                                                      >  1:176020/1‑251 (MQ=255)                                                                                              aaGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATgcgg                                                                                                                           >  1:275249/1‑251 (MQ=255)                                                                                              aaGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCAGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGAGAGAACGCTGAACTGTATgcgg                                                                                                                           >  2:220631/1‑251 (MQ=255)                                                                                                      aaTGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCg                                                                                                                   >  1:29465/1‑251 (MQ=255)                                                                                                       aTGGCAGGCGGTACGTGCGGAGTATCAGCTCCCGCGCGATCCGCTCCATCAGTTTGC‑AGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGTCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGAc                                                                                                                 <  1:305441/251‑1 (MQ=255)                                                                                                                    cGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTttatta                                                                                                     <  1:237483/251‑1 (MQ=255)                                                                                                                     gTGCGGAGTATCAGCGTCAGCGCGATCCGCTACCTCAGTTCCCCAGCCAGCAAAGCCCGGAAGCGCCGCTTCAGGCATTGAAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGCCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACCACGTGGTGAGAACGATGACCTGTATGCGGCGCTGGCGACGGTGCTTTATTAc                                                                                                    >  1:148527/1‑251 (MQ=255)                                                                                                                     gTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGAGTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGTTTATCTGTATGCGGCGCTGGCGACGGTGCTTTATTAc                                                                                                    <  1:264932/251‑1 (MQ=255)                                                                                                                     gTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCAGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTAGCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTAc                                                                                                    <  2:10020/251‑1 (MQ=255)                                                                                                                                 agcgtcagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGccagccag                                                                                        <  1:297861/251‑1 (MQ=255)                                                                                                                                                      aTCAGTCTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTcgc                                                                   >  1:91387/1‑251 (MQ=255)                                                                                                                                                              gccagccaGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATACAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGAt                                                           >  1:296664/1‑251 (MQ=255)                                                                                                                                                              gccagccaGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGa                                                            <  1:163038/250‑1 (MQ=255)                                                                                                                                                                 agccaccAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCg                                                         >  1:254153/1‑250 (MQ=255)                                                                                                                                                                  gccagcAAACCCCGGGAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGAc                                                       <  2:293533/251‑1 (MQ=255)                                                                                                                                                                      gcAAGCCCCGGAAGCGCATCGTCAGGCAATGCATGATAAAATCCTCGCTAATCCACAAAACATCGAACAGGGTGCGTTACTGTATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCTTACCGTCATGCGTGGCAACCGCGTGGTGAGAACGCGGAACTGTATGCGGCGCTGGCGACGGTGCTTGCTTACCAGTCCAGCCAGCATATGTCCGCCCAGACTCGCGCAAtcagcggcaaag                                                   <  2:148527/251‑12 (MQ=255)                                                                                                                                                                                   gcgcAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGa                                      <  2:91387/251‑1 (MQ=255)                                                                                                                                                                                       aGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCtt                                                                                                          >  1:271652/1‑179 (MQ=255)                                                                                                                                                                                       aGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATGACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTTAGAACGCTGAACTGTATGCGGCGCTGGCGACTGTGCtt                                                                                                          <  2:271652/179‑1 (MQ=255)                                                                                                                                                                                         cTTCAGGCATTGCAGGATAAAATACGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTac                                >  2:284976/1‑250 (MQ=255)                                                                                                                                                                                           tCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCACTGATCGACAAAGCCCTCGCGCTGGACAGTAATg                              >  2:56433/1‑251 (MQ=255)                                                                                                                                                                                            aaGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGAGGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGAGCTGGACAGTAATGa                             <  1:37158/250‑1 (MQ=255)                                                                                                                                                                                                                       aaTCCACAAAACAGCGAACAGGGTGCGGTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGGCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACGGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCTCCCTGATGCTACTGGCTTc  <  2:270074/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCGTTATTACCAGGCCAGCCAGCATATGACGGCCCAGACTCGCGCAATTATCGACa                                                      <  2:285608/179‑1 (MQ=255)                                                                                                                                                                                                                                           gTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGACGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACa                                                      >  1:285608/1‑179 (MQ=255)                                                                                                                                                                                                                                                         ||                                                                                                                                                                                                                        GGGCGAAGCGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTC  >  NC_000913/4293533‑4293997

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: ‑