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	94.3%	88.5 / ‑1.9	36	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (17/16);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01
*	NC_000913	4,293,782	0	C	T	94.3%	90.4 / ‑1.2	36	G80D (GGC→GAT)	nrfG	putative formate‑dependent nitrite reductase complex subunit NrfG
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (16/17);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGGCGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTACG  >  NC_000913/4293541‑4294023                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                   cGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACACCACTCCCGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCCGGCGGTACGTTCGGGGTATCAGCGTCAGCGCGATCACCCACATCAGTTTACCAGCCAGCAAACCCCGGCAGCGCCGCTTCAGGCATTGCAGGAAAAAAACCGCGCTGATCCACCAAAAC‑GCGAACAAGCGGG‑GTTACTGGAAGAGTACT‑At                                                                                                                                                                                                                                          >  2:111308/1‑251 (MQ=255)         cAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGc                                                                                                                                                                                                                                   >  1:37907/1‑250 (MQ=255)            cccAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCTCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTTGCaaaa                                                                                                                                                                                                                               <  1:84148/251‑1 (MQ=255)                 aaTACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGAtt                                                                                                                                                                                                                          <  1:321937/251‑1 (MQ=255)                   tACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTAc                                                                                                                                                                                                                        >  1:362369/1‑251 (MQ=255)                     ccGGTGAAAATGCTGACAACACTCACGATATTGCTGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAACCGATTAc                                                                                                                                                                                                                        >  2:33528/1‑249 (MQ=255)                         tGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCACGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAAt                                                                                                                                                                                                                  >  1:87869/1‑251 (MQ=255)                           aaaaTGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTc                                                                                                                                                                                                                <  1:401470/251‑1 (MQ=255)                             aaTGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGGCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGTGGAGTAACAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGt                                                                                                                                                                                                              <  2:352398/251‑1 (MQ=255)                                    acaacaCTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTggg                                                                                                                                                                                                       >  1:170625/1‑250 (MQ=255)                                     caacaCTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGATAAATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAACGAACA‑GTGGGAGTTACTGGATGAGAACT‑ATCTGTGGCAAAACGATTAAAGCAATTCGTTGCTGGCg                                                                                                                                                                                                      >  1:277511/1‑251 (MQ=255)                                       acacTCACGATATTGATGGAATTTCTCTGTGTCGGCAGTTATCTGTTACGCCCAAAATGGAAGGCGGTACGAGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCACGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACTGCGAACA‑GAGGGCGTTACGGGATGAGTACTTAT‑TGGAGCAAAACGAATACAGCAATACGTTGCAGGCGt                                                                                                                                                                                                     >  2:15487/1‑250 (MQ=255)                                                  atTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGt                                                                                                                                                                                         >  1:435916/1‑251 (MQ=255)                                                         gTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCGGCGTCCGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGACCA‑GTGGGCGTTACTGGATGAGAACT‑ATCTGTGGCAAAAAGATTACCGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACt                                                                                                                                                                                  >  2:253260/1‑251 (MQ=255)                                                                                              aaTGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCg                                                                                                                                             >  1:440102/1‑251 (MQ=255)                                                                                                    aGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGAGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTg                                                                                                                                       >  1:370518/1‑251 (MQ=255)                                                                                                             gtgcgtagagtcgtctacagccCGATCCGCGTCATCAGTGTTCCAGCCAGCCAACCCCGGAAGCGCAGATTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑AACTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCATGATTAc                                                                                                                              <  2:355556/234‑1 (MQ=255)                                                                                                             gTGCGTAGTATCACCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTTCGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTAc                                                                                                                              <  2:260273/251‑1 (MQ=255)                                                                                                               gCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCa                                                                                                                            <  1:253981/251‑1 (MQ=255)                                                                                                                 ggAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAgg                                                                                                                          >  1:165232/1‑251 (MQ=255)                                                                                                                      atctgcgtcagcgCGATCCGCTACATCAGGTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGGGTTACTGGATGAGTACT‑AGCTGTGGCAAAACGATTACAGCAATTCGTGGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGAGAGAACGCTGAACTGTATGCGGCGCTTGCGACGGTGCTTTATTACCAGgccag                                                                                                                      <  2:38550/246‑1 (MQ=255)                                                                                                                                gcgcgtcccGCTCCAGCAGTTGGCCAGCCAGCAAACCCCGGAAGCGCAGCGTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGa                                                                                                           <  2:172871/244‑1 (MQ=255)                                                                                                                                    gATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGc                                                                                                       >  1:354484/1‑251 (MQ=255)                                                                                                                                                      gccagccaGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGa                                                                                      <  2:37907/250‑1 (MQ=255)                                                                                                                                                       ccagccagCAAACCCCGGAAGAGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGTCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATc                                                                                    <  2:179126/251‑1 (MQ=255)                                                                                                                                                              gcAAACCCCGGAAGCGCAGCATCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATGCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAg                                                                             <  2:294443/251‑1 (MQ=255)                                                                                                                                                                  aCCCCTGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTAATGGATGAGTACT‑ATCTAAGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGCGAGGACGCTGAACTGTATGCGGCGCTGGCGACGGCGCTTTATTAGCAGGCCAGCCAGCATATGAGCGCCCAGACTCGCGCCATGATCGCCAAAGCCCt                                                                         >  1:82290/1‑251 (MQ=255)                                                                                                                                                                            cgcAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGAc                                                               <  1:33528/251‑1 (MQ=255)                                                                                                                                                                                        cATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGGCGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGCCAGAAATGAAATc                                                   >  2:191965/1‑251 (MQ=255)                                                                                                                                                                                              aGGATAAAATCCGCTACAATCCAC‑AAAACAGCGGATA‑AGAGTCGTTACTGGATGAGGACT‑ATCTGTGGCAAATCGATTACAGCATTTCGTTGCTGGCGTACCGTCGGGCGTTGCAACTACGTGGTGAGAACGCGGAACAGCATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGAGACGCGCAATGATCGATAAAGCCCACTCGCTGGACAGTAATGAAATCACCGcc                                             <  2:444414/251‑1 (MQ=255)                                                                                                                                                                                                       tCCGCGCTAATCCAC‑AAAACAGCGCACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGAtgctg                                    >  2:233084/1‑251 (MQ=255)                                                                                                                                                                                                           cgcTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTCCTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGACCTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCACTGATGCTGCTgg                                >  1:121627/1‑251 (MQ=255)                                                                                                                                                                                                               aaTCCAC‑AAAACAGCTAACA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGGACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCtt                             <  2:435916/250‑1 (MQ=255)                                                                                                                                                                                                                                aaCA‑GTGGGCGTTACTGGATGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTGCCGATGCGTTTATGCAg            <  1:140688/251‑1 (MQ=255)                                                                                                                                                                                                                                     gTGGGCGTGACTGGTTGAGTTCT‑ATCTGTGGCAAAACGATTACGGCAATGCGTTGCTGGTGTACCGTAAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCTCCGCCCTGATGCTGCTGGCTTTCGATGCGTTTATGCAGGCGa        <  2:364480/251‑1 (MQ=255)                                                                                                                                                                                                                                           gTTACTGGGCGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTCCGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCCATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTAcg  >  1:218430/1‑251 (MQ=255)                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                   CGATGAAACAGCCCAAAATACCGGTGAAAATGCTGACAACACTCACGATATTGATGGTATTTCTCTGTGTCGGCAGTTATCTGTTAAGCCCAAAATGGCAGGCGGTACGTGCGGAGTATCAGCGTCAGCGCGATCCGCTACATCAGTTTGCCAGCCAGCAAACCCCGGAAGCGCAGCTTCAGGCATTGCAGGATAAAATCCGCGCTAATCCAC‑AAAACAGCGAACA‑GTGGGCGTTACTGGGCGAGTACT‑ATCTGTGGCAAAACGATTACAGCAATTCGTTGCTGGCGTACCGTCAGGCGTTGCAACTACGTGGTGAGAACGCTGAACTGTATGCGGCGCTGGCGACGGTGCTTTATTACCAGGCCAGCCAGCATATGACCGCCCAGACTCGCGCAATGATCGACAAAGCCCTCGCGCTGGACAGTAATGAAATCACCGCCCTGATGCTGCTGGCTTCCGATGCGTTTATGCAGGCGAACTACG  >  NC_000913/4293541‑4294023

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