Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,163,984	C→T	L116F (CTT→TTT)	btuB →	cobalamin outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,163,984	0	C	T	97.7%	120.2 / NA	43	L116F (CTT→TTT)	btuB	cobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (16/26);  minor base G (0/1);  total (16/27)
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.90e-01

TTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATG  >  NC_000913/4163742‑4164227                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     ttGAACAGCCGCTCATCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑a                                                                                                                                                                                                                                              <  1:20446/250‑1 (MQ=255)  tGAACAGCCGCGCAGCACTGTTCTTGCACCAACCCCCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTTCCGACTTTAGCC‑Ag                                                                                                                                                                                                                                             <  2:52871/250‑1 (MQ=255)  tGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCCGAATCTGGCGGGGGGGAGGGGTTTTTCCGACTTTAGCC‑AGt                                                                                                                                                                                                                                            >  1:494381/1‑251 (MQ=255)       aGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCt                                                                                                                                                                                                                                       <  2:16447/251‑1 (MQ=255)       aGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGGGAGTGGGTCTGCCGGCTTTTTCC‑AGTTCCCt                                                                                                                                                                                                                                       >  1:478892/1‑251 (MQ=255)                            cACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGtt                                                                                                                                                                                                                  <  2:431331/251‑1 (MQ=255)                              ccaaccaCCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGTGTGAGTGGTTCTGTCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGa                                                                                                                                                                                                                <  2:189862/251‑1 (MQ=255)                                aaccaCCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGGCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAat                                                                                                                                                                                                              <  2:51969/251‑1 (MQ=255)                                            gACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGGCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATGTCCGTGGGc                                                                                                                                                                                                  <  2:633871/251‑1 (MQ=255)                                             aCCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGcc                                                                                                                                                                                                 <  1:206453/251‑1 (MQ=255)                                                   cAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCCGGTGAGCTCTCCTTTTTTTTTATTTGCGGTACAAAAGCCAGTCATGGGTTGGGGGTTATTG‑ATTGCGTACCCCTGAATCTGGCGGGGGGGGGTGGTTCTGCCGGCTTTTGCCCGGTTCCCT‑TTGCGCTTGTCCAGCGGGGTGAATATATTCGGGGGCcgcggcc                                                                                                                                                                                           >  2:316146/1‑248 (MQ=255)                                                           cGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGTTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGttt                                                                                                                                                                                   <  1:646675/251‑1 (MQ=255)                                                            gaccgccagaaGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCAGGCGTCGATATCACCCAAAACGGCGGTTCAGGACAGCTCTCATCCCTTTTTATTTTCGGTACCAATGCCAGTCATGTGTTGGTGTTACTTT‑ATGCCGTAAGCCTGAATCTGGCGGGGGTGAGTGGTTCTACCGACTTTACCC‑AGTTCACTATTGCGCTTGTGCAGCGTGTTAGATATAACCGCGGGCGGCGCTCCCCTGtttt                                                                                                                                                                                  >  2:414380/11‑250 (MQ=255)                                                                            ccTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATTGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATa                                                                                                                                                                  <  1:559927/251‑1 (MQ=255)                                                                                 gTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGGACAAATGCCAGTCATGTGTTGGTGTTTATTG‑ATGGCGGACGCCTGAATCTGGCGGGGGTGAGGGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGGTGAATATATCCGTGGGCCGCGCCCCGCTGTTTATGGTTCCGGTGccaaaagcgg                                                                                                                                                             >  1:185653/1‑242 (MQ=255)                                                                                       gaggtgCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAGCGGCGGTTCAGGTCAGTTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCTGggtggt                                                                                                                                                       <  2:466252/248‑1 (MQ=255)                                                                                             ctgcGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTTTTCGCGGTACAAATGCCCGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATTCGTGGGCCGCGCTCCGCTGTTTTTGGTTCCGATGCAATAGGCggggggggggatat                                                                                                                                                 >  1:199675/1‑241 (MQ=255)                                                                                                cgcCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATatcat                                                                                                                                              <  2:199675/251‑1 (MQ=255)                                                                                                 gcCGGCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTTAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATatcatc                                                                                                                                             <  2:500638/251‑1 (MQ=255)                                                                                                                 cGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAAGATGCGGGGTGGTGAATATCATCACGACGCGCGATGAAc                                                                                                                             <  2:379602/251‑1 (MQ=255)                                                                                                                     atCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCcgg                                                                                                                         <  1:278384/251‑1 (MQ=255)                                                                                                                       cACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCcggaa                                                                                                                       <  2:110166/251‑1 (MQ=255)                                                                                                                          ccAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCcggaacgg                                                                                                                    <  2:494381/251‑1 (MQ=255)                                                                                                                                 ggcggTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAAtttt                                                                                                             >  1:544862/1‑250 (MQ=255)                                                                                                                                  gcggTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGTGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTca                                                                                                            <  2:236764/251‑1 (MQ=255)                                                                                                                                    ggTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGTTGTTAATTG‑ATGGCGTGCGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGGCGTTAGCC‑GGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTGCCGATGCAATAGGCGGGGTATTTAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcagc                                                                                                          <  1:48433/251‑1 (MQ=255)                                                                                                                                       tCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTcagcag                                                                                                        <  1:12576/250‑1 (MQ=255)                                                                                                                                                     aTCTATTTTTATTCGAGGTACAAATGCCAGTCATGTGTTGGTGTTAATTT‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCCGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCACTAGGCGGGGTGGGGCAAATCATCACCACGCGCGATGCACCCGGAACGGAAAGTTCCTGAAGGTGGGGAAGCAAt                                                                                          >  2:315377/1‑250 (MQ=255)                                                                                                                                                      tCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAAt                                                                                          >  1:550058/1‑249 (MQ=255)                                                                                                                                                         aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTTCGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCTTGTGCCGCGCTCCGCTGTTTATGGTTCCGAGTCATGAGGCTGGGTGGTTAATATCATCACGATACGCGATTTGCCCGGAACGGAAATTGCAGGAGGTTGATGAAGCAATa                                                                                         <  2:87478/247‑1 (MQ=255)                                                                                                                                                                  cgcggttcaaattgccgttaTGTGTTGGGGTTAATTGGATGCGGCAC‑CCTGAATCCGGCCGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTCCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTa                                                                            <  2:361318/232‑1 (MQ=255)                                                                                                                                                                     ggTACAAATGCCAGTCATGTGTTGGTATTACTTG‑CTGGCGTACGCCTGAATCTGGCGGCGGTGAGTGGTTCTGCCGACTTTAGCC‑CTTTCCCTATTGCGCTTGTCCCGAGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGCTGCAATAGGCGGGGTGGGGAATATCGTCACGACGCGCCATGAACCCGGCACAGAAGTTTCAGCCGGGTGGGGCAGGAATAATTATccgtaccttcc                                                                         >  2:481512/1‑241 (MQ=255)                                                                                                                                                                            aTGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCGGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACTCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGAGAAGCAATAGTTATCAGAACTATGATGTctct                                                                  <  1:438201/251‑1 (MQ=255)                                                                                                                                                                                    cATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGt                                                                                                                                                                                     >  1:606756/1‑128 (MQ=255)                                                                                                                                                                                    cATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGt                                                                                                                                                                                     <  2:606756/128‑1 (MQ=255)                                                                                                                                                                                             gtgtTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑CGTTCCCTATTGCGCTTGTCCCGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGAAGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTggg                                                 >  2:66851/1‑251 (MQ=255)                                                                                                                                                                                                ttAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAAACCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTTTCTACGCAGCAACAACTggggg                                               >  2:14729/1‑250 (MQ=255)                                                                                                                                                                                                  caTTG‑ATGGCGTACGCCTGAATCTGGCGGGGGAGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGCACTATGATGTATCTACGCAGAAACAACTGGGGGAt                                             >  2:382039/2‑250 (MQ=255)                                                                                                                                                                                                   aTTG‑ATGGCGTACGCCTGCATCAGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGCGGTGAAAATCATAACGACGGGCGATGAACCCCGAACGGAAGTTTTAGCAGGGTGGGGAAGCAATAGTTATCAGAAATATGAGTTTTCTACGCAGGAACAACTgggggtta                                            >  2:78553/1‑247 (MQ=255)                                                                                                                                                                                                                cGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGAt                                                                                                                                                                        <  2:282796/114‑1 (MQ=255)                                                                                                                                                                                                                cGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGAt                                                                                                                                                                        >  1:282796/1‑114 (MQ=255)                                                                                                                                                                                                                           ggCGGGGGTGAGTGGTTCTGCTGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg                      <  2:397689/249‑1 (MQ=255)                                                                                                                                                                                                                                              gCCGACTTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTTAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTcgtg  <  1:78553/250‑4 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     TTGAACAGCCGCGCAGCACTGTGCTTGCACCAACCACCGTTGTGACCCGTCAGGATATCGACCGCTGGCAGTCGACCTCGGTCAATGATGTGCTGCGCCGTCTTCCGGGCGTCGATATCACCCAAAACGGCGGTTCAGGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTG‑ATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCC‑AGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATG  >  NC_000913/4163742‑4164227

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