Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,891,085	G→A	P272S (CCG→TCG)	cysJ ←	sulfite reductase, flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,891,085	0	G	A	98.5%	183.4 / NA	66	P272S (CCG→TCG)	cysJ	sulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/33);  minor base C (1/0);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGT  >  NC_000913/2890842‑2891330                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      gCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCATTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTa                                                                                                                                                                                                                                                 <  1:430635/250‑1 (MQ=255)  cTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAAc                                                                                                                                                                                                                                               >  1:59295/1‑251 (MQ=255)  cTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGAGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAAc                                                                                                                                                                                                                                               >  1:453894/1‑251 (MQ=255)      aaCTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCAt                                                                                                                                                                                                                                           <  2:71329/251‑1 (MQ=255)         ttCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCAGATACCAGACGACCAGCGCGTCACCCGACTGGTAACGCATGcc                                                                                                                                                                                                                                        >  1:17097/1‑251 (MQ=255)              tttATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGAAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGt                                                                                                                                                                                                                                   >  1:52941/1‑251 (MQ=255)                tatcgccgtgcatcGGCAGTAGTGTTGCACTGCGGGTACGCGTGTCGTTTTTCTCAAGAAGGTGGGCGGTGTTGACGGTCAGTTCGACGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGGCGGTGACAGGTTCATCGCCGTTCAGCCACAGGAGTTTGACAAGTTCTTTCACCAGTGCCGGAGCGTTCTGATACCAGACGCTCAGCGCGTCACCCGACTGGTAACGCATGCACGAGTc                                                                                                                                                                                                                                  <  1:175847/237‑1 (MQ=255)                   cGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGGTCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGGTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCAcc                                                                                                                                                                                                                               <  1:318798/250‑1 (MQ=255)                      ccaccaGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAg                                                                                                                                                                                                                           <  1:328915/251‑1 (MQ=255)                           cgCGGCAGCAGTGTGTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGTCGGTGTTGACGGTCATTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCGCCCGACTGTTAACTCATGCCCGAGTCACCTAAGTCAAt                                                                                                                                                                                                                      <  2:369919/250‑1 (MQ=255)                           aGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCTTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAAt                                                                                                                                                                                                                      <  2:14852/251‑1 (MQ=255)                               gcagcaGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGATCAACCCCTTTCAGCCAAAGCAGTTCGACAAGTTCTTTCACCAGTTCCGGATCGTTCTGATACCAGACACCCATCGCGTCACCCGACTGGTCACGCATGCCCGAGTCACCTAAGTCAAttt                                                                                                                                                                                                                    >  2:252598/1‑249 (MQ=255)                                cagcagTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCg                                                                                                                                                                                                                  >  1:346830/1‑250 (MQ=255)                                   cagggtTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCCGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGCTACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGCGTCACCTAAGTCAATTTCAatat                                                                                                                                                                                                              >  2:227441/5‑251 (MQ=255)                                    agtgttgcatttcggttAAGCTTGGGGTAATTCTCAACAATGGTGGGGGTGTTGAAGAGCAGTTCGAGGTGCCACTGTTGCGCTTCGTTCAGATGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCCCAGCAGTTCGACAAGTTCTTTCACCAGGGCCGGATCATTCTCATACCAGACGACGAGCGCGTCACCCGACTGGTAACGCATGCCCGCGTCAACTAAGTCAATTTGGatatt                                                                                                                                                                                                             <  2:230584/235‑2 (MQ=255)                                    agTGTTTCACTGCTGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTTTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATg                                                                                                                                                                                                             <  2:447851/251‑1 (MQ=255)                                       gtTTCACTGCGGGTAAGCGTGGCTTAATTCGCAACAATGTTGGCGGGGTTGACGGGCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCg                                                                                                                                                                                                          <  2:265901/251‑1 (MQ=255)                                         ttCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGGCAGTCAGTTCGAAGTACCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCATCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCCCCAGTGCAGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGACACCTAAGTCAATTTCGATATGGCGa                                                                                                                                                                                                         >  1:289375/1‑250 (MQ=255)                                           cACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGGTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCTCAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGAGATGGCGAAc                                                                                                                                                                                                       <  2:317821/250‑1 (MQ=255)                                                    tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGACGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCCGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTc                                                                                                                                                                                              >  1:372044/1‑250 (MQ=255)                                                    tAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTc                                                                                                                                                                                              >  1:345334/1‑250 (MQ=255)                                                    tAAGCGAGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACCGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTc                                                                                                                                                                                              >  1:367644/1‑250 (MQ=255)                                                     gaGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTc                                                                                                                                                                                              >  2:109603/2‑249 (MQ=255)                                                     aaGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCa                                                                                                                                                                                             <  1:105426/250‑1 (MQ=255)                                                      agcttgtcGTAATGCTAAACAATGTTGGCGGTGGTGACGGTCAGTTCGAAGGGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTGGGCCCGCGACGGTGACAGTTTAATCGCCTTGCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGGGGCGGATCGTTCTGATTCCAGACTCCCAGCGTGTCGCCCGACTTGGAACGCATGCTCGAGTTACCTAAGGCGATTGCGATATGGCGAAGTTCCTTTTCaga                                                                                                                                                                                           <  2:221543/244‑1 (MQ=255)                                                      aGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTc                                                                                                                                                                                              >  1:400822/1‑248 (MQ=255)                                                      aGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTc                                                                                                                                                                                              >  1:444505/1‑248 (MQ=255)                                                       gtgTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACTGTCAGGTCGAAGTGCCAATGTAGCGCTTCGTTCATAGGCACCGTTTTTCCCTCGACGGTGAAAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCTTCACCCGACTGGTAACGCATGCCCGAGGCACCTAAGACAATTTCGATATGGCGAACGTCTTTTTCagag                                                                                                                                                                                          <  2:361071/249‑1 (MQ=255)                                                          tGGCGGAATTCGCAACAAGGTTGTCGGTGTTGCCGGGCAGTTCGAAGGGCCCCTGTAGCGCTTCTGTCAGATGCAACGTTTTGCCCTCTACGGTGTTAGGTTCATCGCCTTTCAGCCACAGCAATGCGTCAAGTTCTTTCTCCAGTGCCGGATCGTTCTGATACCAGACGTGCAGCGCTTCACTCGTGGGGTAGCGCAGTCCTGAGTCGCCTAAGGCAATTGGGATATTGCTGACGTCTTGTTCAGAGtt                                                                                                                                                                                        <  2:410476/250‑1 (MQ=255)                                                           ggCGTAATTCTCAACAATGTTGGTGGTGTTGACGGTCAGTTCGAAGTGCTACTGTAGCGCTTCGTTCAGAGTCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCATTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTGCGATATGGCGAACGTCTTTTTCAGAGTTAc                                                                                                                                                                                      <  2:177800/251‑1 (MQ=255)                                                           ggCGTAATTCGCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTAc                                                                                                                                                                                      >  1:68957/1‑251 (MQ=255)                                                                aaTTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCg                                                                                                                                                                                 <  1:217564/251‑1 (MQ=255)                                                                    ctcAACAATGGTGGCGGTGTTGACGGTCAGTTAGAAGTGCCCCTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGATTCATCGACTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGATCTGAAACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGACACCTAAGTCACTTTCGCTATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTa                                                                                                                                                                              >  1:372419/1‑250 (MQ=255)                                                                       aacaaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTAAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCACCTTTCAGCAACAGCAGTTCCACAAGTTCTTTCACCCGTGCCGGATCGTTCTGATCCCAGACGCCCAACGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTCAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTACttt                                                                                                                                                                          >  1:288744/1‑251 (MQ=255)                                                                         caaTGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTc                                                                                                                                                                        >  1:102201/1‑251 (MQ=255)                                                                                        ggACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAgg                                                                                                                                                         <  2:408865/250‑1 (MQ=255)                                                                                         gACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGGTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCGCCCGACTGGGAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTGTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGc                                                                                                                                                        <  2:360706/251‑1 (MQ=255)                                                                                                  ttCGAAGTGCCACCGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACCGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGCGAGGCTAGccacc                                                                                                                                                >  2:286489/1‑250 (MQ=255)                                                                                                      aaGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAgcgg                                                                                                                                           <  1:100944/251‑1 (MQ=255)                                                                                                          gCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGcgcg                                                                                                                                       <  1:92264/251‑1 (MQ=255)                                                                                                             aCTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCt                                                                                                                                    >  1:11375/1‑251 (MQ=255)                                                                                                                    cgcTTCGTTCAGAGGCAACTTTTTGCCCGCGGCGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTGCGGCAATTTCGTTCACCCGTGCCGGGTCGTTCTGATTCCAGACGCCCAGCGCGTCACCCGACTGGGAACGCATGCCCGAGGCCCCTAAGTCAATTTCGAAATGGCGAACGGCTTGTTCAGAGTTACGCCCGGTAATTTTCTGGGTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGgtgt                                                                                                                             <  2:37340/251‑1 (MQ=255)                                                                                                                     gcTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGt                                                                                                                             >  1:400578/1‑250 (MQ=255)                                                                                                                      cTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGAGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTAc                                                                                                                           >  1:221827/1‑251 (MQ=255)                                                                                                                          gTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACgtgc                                                                                                                       >  1:35584/1‑248 (MQ=255)                                                                                                                                  gCAACTTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCTGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCGTTTTCAGAGTTACGCCCGGTAATTTTCTGGGTAACAGAGAGGCTAGCCACCATCGGCGCGTCTTTGCTGTACGGGCTGGTGTgg                                                                                                               <  1:447352/251‑1 (MQ=255)                                                                                                                                       gTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTc                                                                                                          >  1:423524/1‑251 (MQ=255)                                                                                                                                         tttGCCCTCGACGGCGACAGGTTCATCGCCTTTCAGCCACAACAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGACTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGCGAGGCTAGCCCCCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGGTTTCAt                                                                                                        >  2:70715/1‑251 (MQ=255)                                                                                                                                                   aCGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGcc                                                                                              <  1:312627/251‑1 (MQ=255)                                                                                                                                                        gACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGTGCTGGTGTGGATTTCATTTACCGCGCCAGTa                                                                                          <  1:412572/250‑1 (MQ=255)                                                                                                                                                                cATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGAt                                                                                 >  2:144049/1‑251 (MQ=255)                                                                                                                                                                           aGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGc                                                                      <  1:204042/251‑1 (MQ=255)                                                                                                                                                                                  gcagTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACa                                                               <  1:404400/251‑1 (MQ=255)                                                                                                                                                                                       gcGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATGCCAGACGCCCAGCTCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTTAGTTTCTGGTTAACAGAGAGGCTAGCCACCATCGGCGCGTCTTTGCTGTATGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGTGACGGATTGCGAAGGTGCCGCGACAGgcg                                                           <  2:453894/249‑1 (MQ=255)                                                                                                                                                                                          acacgTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACTCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTTCGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCAcg                                                       <  2:102201/247‑1 (MQ=255)                                                                                                                                                                                                tCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGAttt                                                 >  1:350329/1‑251 (MQ=255)                                                                                                                                                                                                       ccAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAgagc                                           >  1:149319/1‑247 (MQ=255)                                                                                                                                                                                                         aGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTGCATTTACCGCGCCAGTAGCGGCGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATc                                        <  1:413002/251‑1 (MQ=255)                                                                                                                                                                                                          gTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTCATTTTCTGGGTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGAGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGTTTTAAGCGCATc                                        >  1:51783/1‑250 (MQ=255)                                                                                                                                                                                                              cGGATCGTTCTGATACCAGACGCCCAGCGCGTCCCCCGACTGTTAACGCATGCCCGAGGCACCTAAGTCAATTGCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTGCTGGTTAACAGAGAGGCTAGCCACCAGCGGCTCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCTCGACAGGCTCACGCTATTTAAGCGCATCAACCa                                   <  2:201623/251‑1 (MQ=255)                                                                                                                                                                                                              cGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGCAGGTGCCGCGACAGGCGCACGCGGTTTAAgcgcctccaccc                                   >  2:368018/1‑243 (MQ=255)                                                                                                                                                                                                               ggATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGTGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCGAGCCTCCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCAc                                  <  1:428608/251‑1 (MQ=255)                                                                                                                                                                                                                         gATACCAGACGCCCAGCGCGTCACCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGgcg                           <  1:429661/248‑1 (MQ=255)                                                                                                                                                                                                                              cAGACGCCCAGCTCGTCCCCCGACTGGTAACGCATGCCCGAGTCACCTAAGTCAAGTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCACCCACGCGGGCGCGCCACTc                   <  2:393444/251‑1 (MQ=255)                                                                                                                                                                                                                              cAGACGCCCAGCGCGTCACACGACTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGCGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGAGTGGATTTCATTTACCGCGCCAGTAGCGCCGGATTGCGAAGAGGCCGCGACAGGCGCCCGCGATTTAAGCGCAACAAACACGCGGGCGCGCCCCTc                   >  2:175333/1‑251 (MQ=255)                                                                                                                                                                                                                                               aCCCGCCTGGTAACGCTTGCCCGAGTCACCTAAGTCCATTTCAATATGGCGAACGTCTTTTTCAGCGTTCCGCCCGGAAATTTTCTGGATCACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGGTGTACGGGCTGGTGTGGATTTCATTTACCGCCCCAGTCGCGGCGGCTTGCGGACGTGGCGCGACACGAGCCCGCGATTTTAGCGCATCCACCACGCAGGCGCGCCACTCGCAGgcagccacctggg  >  2:144366/1‑243 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      GCTGTAACTTCGCTTTATCGCCCACCAGCGGCAGCAGTGTTTCACTGCGGGTAAGCGTGGCGTAATTCTCAACAATGTTGGCGGTGTTGACGGTCAGTTCGAAGTGCCACTGTAGCGCTTCGTTCAGAGGCAACGTTTTGCCCTCGACGGTGACAGGTTCATCGCCTTTCAGCCACAGCAGTTCGACAAGTTCTTTCACCAGTGCCGGATCGTTCTGATACCAGACGCCCAGCGCGTCACCCGGCTGGTAACGCATGCCCGAGTCACCTAAGTCAATTTCGATATGGCGAACGTCTTTTTCAGAGTTACGCCCGGTAATTTTCTGGTTAACAGAGAGGCTAGCCACCAGCGGCGCGTCTTTGCTGTACGGGCTGGTGTGGATTTCATTTACCGCGCCAGTAGCGACGGATTGCGAAGGTGCCGCGACAGGCGCACGCGATTTAAGCGCATCAACCACGCGGGCGCGCCACTCGCTGGCAGCAGCCTGGT  >  NC_000913/2890842‑2891330

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