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%	144.9 / NA	42	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 (16/26);  total (16/26)
*	NC_000913	4,293,782	0	C	T	100.0%	143.2 / NA	42	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/26);  total (16/26)

GATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTAC  >  NC_000913/4293542‑4294022                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                                 gATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATc                                                                                                                                                                                                                                        <  2:79197/251‑1 (MQ=255)          gCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCaaa                                                                                                                                                                                                                               <  2:109836/251‑1 (MQ=255)            ccAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAAc                                                                                                                                                                                                                             <  1:257041/251‑1 (MQ=255)                                  gACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTgg                                                                                                                                                                                                       >  1:265283/1‑251 (MQ=255)                                           cACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTc                                                                                                                                                                                              <  1:266637/251‑1 (MQ=255)                                                   tGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTg                                                                                                                                                                                      >  1:39658/1‑251 (MQ=255)                                                     aTGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGc                                                                                                                                                                                     >  1:113116/1‑250 (MQ=255)                                                       ggTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAAc                                                                                                                                                                                  <  1:74819/251‑1 (MQ=255)                                                         tATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACt                                                                                                                                                                                 >  2:105772/1‑250 (MQ=255)                                                          aTTTCTCTGTGTCGTCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCTGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCTCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGTCGTACCGTCAGGCGTTGCAACTa                                                                                                                                                                                <  1:331087/250‑1 (MQ=255)                                                                  gtgtCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTga                                                                                                                                                                        <  2:223755/250‑1 (MQ=255)                                                                    gtCGGCCGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTgaga                                                                                                                                                                      >  2:339050/1‑250 (MQ=255)                                                                                  gTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTAt                                                                                                                                                        <  1:271177/250‑1 (MQ=255)                                                                                    tAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTgg                                                                                                                                                                                                       <  2:367965/201‑1 (MQ=255)                                                                                    tAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTgg                                                                                                                                                                                                       >  1:367965/1‑201 (MQ=255)                                                                                     aaGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATgcg                                                                                                                                                     <  1:45524/250‑1 (MQ=255)                                                                                          cAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTg                                                                                                                                               >  1:112669/1‑251 (MQ=255)                                                                                          cAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTg                                                                                                                                               <  1:121479/251‑1 (MQ=255)                                                                                            aaaTGGCAGGCGGTACGTGCGGAGGATCATCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCATGCGTTGCAATTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTgg                                                                                                                                              <  1:359350/250‑1 (MQ=255)                                                                                              aTGGCAGGCTGTACGTGCGGAGTATCAGCGGCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAAGCCACAAAACATCGAACAGTGGGCGTTACTGGATGAGTTCTATCTGTGGCAAAACGATTACAGCAATTCGGTGCTGGCGGACCGTCAGGCGGTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGa                                                                                                                                           <  1:170201/251‑1 (MQ=255)                                                                                                    ggCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGG‑CGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCt                                                                                                                                    <  1:306510/251‑1 (MQ=255)                                                                                                     gCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCt                                                                                                                                    <  1:78402/251‑1 (MQ=255)                                                                                                                        agcgtcagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGgccagcc                                                                                                                   >  1:6295/1‑249 (MQ=255)                                                                                                                             cagcgCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCata                                                                                                             <  2:112669/250‑1 (MQ=255)                                                                                                                                      ccGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCGGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGccc                                                                                                    <  2:131641/250‑1 (MQ=255)                                                                                                                                           aCATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTc                                                                                              >  2:120156/1‑251 (MQ=255)                                                                                                                                           aCATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTAAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTc                                                                                              >  1:255819/1‑251 (MQ=255)                                                                                                                                               cAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATGACCAGGCCAGCCAGCATATGACCGCCCAGACTcgcg                                                                                           <  2:103783/250‑1 (MQ=255)                                                                                                                                                                                cTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTCGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTa                                                          >  2:258837/1‑250 (MQ=255)                                                                                                                                                                                     ggCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGaaa                                                    <  2:68197/251‑1 (MQ=255)                                                                                                                                                                                       cATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATc                                                  <  1:152148/251‑1 (MQ=255)                                                                                                                                                                                           gCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGTGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCg                                              <  2:239995/251‑1 (MQ=255)                                                                                                                                                                                               gATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCt                                          <  2:39658/251‑1 (MQ=255)                                                                                                                                                                                                            cTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTCTCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGCACTGTCTGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCt                             >  1:164902/1‑251 (MQ=255)                                                                                                                                                                                                            cTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTACTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCt                             >  2:178607/1‑251 (MQ=255)                                                                                                                                                                                                                  cacaAAACAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTa                                                                                                                                                         >  1:332381/1‑121 (MQ=255)                                                                                                                                                                                                                  cacaAAACAGCGAACAGTGGGCGTTACTGGATGAGTACGATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTa                                                                                                                                                         <  2:332381/121‑1 (MQ=255)                                                                                                                                                                                                                     aaaaCAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGGCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCg                    <  2:375726/251‑1 (MQ=255)                                                                                                                                                                                                                       aaCAGCGAACAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGtt                  >  1:88180/1‑251 (MQ=255)                                                                                                                                                                                                                               aCAGTGGGCGTTACTGGATGAGTACTATCTGTGGCAAAAAGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCACGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAgg          >  1:39961/1‑251 (MQ=255)                                                                                                                                                                                                                                 ggTGGGCGTGACTGGATGAGTACTATCTGTTGCAAAACGATTACAGCAAGTCTGTGATGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCGGGCGACGGTGCTTTATTACCAGGCCAGCCCGCATATGACCGCCCAGACTCGCGCAATGAGCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCTCCTCCCTTATGCTGCTGGCTTCCTATGCGTTTATTCAGGCg        <  2:174369/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggTTACTGGATGAGTACGAGCTGTGGCAAAACGATTACAGCAATTCGTTGCGGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTAGGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAAGCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTAc  <  2:85989/250‑1 (MQ=255)                                                                                                                                                                                                                                                ||                                                                                                                                                                                                                                                 GATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCACAAAACAGCGAACAGTGGGCGTTACTGGGCGAGTACTATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTAC  >  NC_000913/4293542‑4294022

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