Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,164,123	C→G	A162G (GCA→GGA)	btuB →	cobalamin outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,164,123	0	C	G	97.7%	153.4 / NA	44	A162G (GCA→GGA)	btuB	cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (21/22);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑CAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCT  >  NC_000913/4163882‑4164347                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 cAGCTCTCATCTATTTTGATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCGGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTGCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTgggg                                                                                                                                                                                                                         <  2:375148/251‑1 (MQ=255)  aGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGGACGGAAATTTCAG‑GAGGGTgggg                                                                                                                                                                                                                         >  1:216727/1‑250 (MQ=255)   gCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGCACCCGGAACGGGAATTTCAG‑GAGGGTGGGGa                                                                                                                                                                                                                        >  1:340284/1‑250 (MQ=255)            tATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGAGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCCCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCCTCACGACGCGCGATGAAACCGGAACGGAGATTTCAG‑GAGGGTGGGGAAGCAATatgt                                                                                                                                                                                                              >  1:311547/1‑248 (MQ=255)             aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAGTTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATTTTTCCGTGGGGCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGATCGGAAAGTTCAG‑GAGGGTGGGGAAGCAATa                                                                                                                                                                                                                 <  2:213556/247‑1 (MQ=255)              tttttATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGCAGCAATAGTTAt                                                                                                                                                                                                            >  1:71389/1‑251 (MQ=255)                        cgGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTGCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                                                                                                                                                                  <  2:213688/251‑1 (MQ=255)                         ggTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTatg                                                                                                                                                                                                  >  1:26564/1‑250 (MQ=255)                          gTACAAATGCCAGCCATGTGTTGGAGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCAAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATACGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATCTCATCACGACGCGCGATGAACCCGGAACGGGAATTTCCG‑GAGGGGGGGGAAGCAATAGTTAACAGAACTatga                                                                                                                                                                                                 >  2:355912/1‑250 (MQ=255)                               aaTGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACCACGCGCGATGAACCCGGAACGGAAATTTCAGGCAGG‑TGGGGAAGCAATAGGTATCAGAACTAAGGTGtctc                                                                                                                                                                                           >  2:259742/1‑251 (MQ=255)                                aTGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTctct                                                                                                                                                                                          >  1:99559/1‑251 (MQ=255)                                      gtcattttttGGTGTTAATGGAGGGCGGACGCTTGAATCTGGCGGGGGTGAGTGGTTCTGCCGTCCTTATCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTGTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GGGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACgcag                                                                                                                                                                                    <  2:362416/243‑1 (MQ=255)                                                         ggATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTGATCAGATCTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGa                                                                                                                                                                 <  2:331738/250‑1 (MQ=255)                                                          gATGGCGAACGCCTGAATCTGGCGGGGGGGACTGGTTCTGCCGACCTTAGCCAGTTCCCTATCGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATaaaac                                                                                                                                                                >  2:256414/1‑248 (MQ=255)                                                                gTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCCGCAACAACTGGGGGATAAGACACGGGt                                                                                                                                                          >  1:82918/1‑251 (MQ=255)                                                                    gCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACg                                                                                                                                                      <  2:378800/251‑1 (MQ=255)                                                                              ggCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGAGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GGGGGGGGGGAAGCCATAGGTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCCGATGGGCg                                                                                                                                            >  1:122997/1‑251 (MQ=255)                                                                                 gggggTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAtt                                                                                                                                         >  1:124365/1‑251 (MQ=255)                                                                                    ggTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATg                                                                                                                                      <  2:216727/251‑1 (MQ=255)                                                                                                 gCCGACCTTAGCCAGTTCCCTATTGCTCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCa                                                                                                                            <  1:206665/248‑1 (MQ=255)                                                                                                                 tCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTACGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGACTATCAGCACGACGCGCGATGAACCCGGAACGGAAATTTCTG‑GAGGGGGGTGAAGCAATAGTTATCAGAACTATGATGTATCTAAGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGGTGcc                                                                                                         >  2:355479/1‑251 (MQ=255)                                                                                                                  cgcTATTGCGCTTGTCCAGCGGGGTGAATAGCTCCGTGGGCCGCTCGCCGCTTGTTATTGTTTCGATGCAATAGGCTGGGTGGTGAATATCATCACGACGCGAGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACGATGATGTCTCTACGCAGCAACAACTGGGGGATAATACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGTTTATGATTTTTTTGcc                                                                                                         <  2:245991/248‑1 (MQ=255)                                                                                                                        tGCGCTTGTCCAGCGTGTTGAATATATCCGTGGACCGCGCTCCGCTGTTTATGGTTCCGATGCAATAAGCGGGGTGGTGAATATCCTCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCACCAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATCCTCATGGTTATGATGTTGTTGCCTATGGt                                                                                                   >  1:22847/1‑250 (MQ=255)                                                                                                                               gTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGATCCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGTGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCg                                                                                            <  1:320369/250‑1 (MQ=255)                                                                                                                                        gtTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAgc                                                                                  <  2:124365/251‑1 (MQ=255)                                                                                                                                          tGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAgcg                                                                                 <  1:354743/250‑1 (MQ=255)                                                                                                                                            aaTATATCCGTGGGCCGCGCTCCGCTTTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATCCCGGAACGCAAGCGcag                                                                              <  1:341549/251‑1 (MQ=255)                                                                                                                                                     gTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAAc                                                                     >  1:255935/1‑251 (MQ=255)                                                                                                                                                        ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCCG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGAt                                                                  >  1:400959/1‑251 (MQ=255)                                                                                                                                                        ggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGAt                                                                  <  2:255935/251‑1 (MQ=255)                                                                                                                                                               cTCCGCTGTTTATGGTTCCGCTGCAATAGGCGGGGTGGTGAATCTCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GCGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGCTAAGACACGGGAAACCCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACTGGAACGCAAGCGCAGCCAGATAACGCTGGttttt                                                           >  2:155924/1‑251 (MQ=255)                                                                                                                                                                 ccGCTGTTTATGGTTCCGATGCAAGAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACTGATAACGATGGtttttt                                                          <  2:270294/250‑1 (MQ=255)                                                                                                                                                                  ggCTGTTTATGGTTCCGATGCAGTAGGCGGGGCGGTGAATATCATCACGACGTGCGATGAACCCGGAACGGCAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTa                                                         >  2:40918/2‑250 (MQ=255)                                                                                                                                                                        ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg                                                                                                                                                                       >  2:309955/1‑134 (MQ=255)                                                                                                                                                                        ttATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggggg                                                                                                                                                                       <  1:309955/134‑1 (MQ=255)                                                                                                                                                                          aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACg                                                <  2:19479/251‑1 (MQ=255)                                                                                                                                                                            ggTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGTTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCt                                              <  1:152061/251‑1 (MQ=255)                                                                                                                                                                                 cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATg                                         <  1:92379/251‑1 (MQ=255)                                                                                                                                                                                      cAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACTCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGgcgcg                                    <  1:259742/251‑1 (MQ=255)                                                                                                                                                                                                                  gacgCGCGATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGCTGCCTATGGTAATACCGGAACGCACGCGCAGACAGGTAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTgg         >  1:99171/1‑250 (MQ=255)                                                                                                                                                                                                                    cgcgcgATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGc      <  1:185820/251‑1 (MQ=255)                                                                                                                                                                                                                       gcgATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGc   >  2:239191/1‑251 (MQ=255)                                                                                                                                                                                                                       gcgATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGc   <  1:71136/251‑1 (MQ=255)                                                                                                                                                                                                                        cgATGAACCCGGAACGGAAATTTCAG‑GAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCt  >  1:122264/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 CAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAG‑CAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCT  >  NC_000913/4163882‑4164347

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