Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,538,387	C→T	R156H (CGC→CAC)	narY ←	nitrate reductase Z subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,538,387	0	C	T	97.9%	144.5 / NA	49	R156H (CGC→CAC)	narY	nitrate reductase Z subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (21/26);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTGCGGCACCACCGGGTTGGCGAAGATTTTGGTTATTACGCCCATCTTGTTGCCCAGACGCGG  >  NC_000913/1538142‑1538628                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  cGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCa                                                                                                                                                                                                                                              <  2:31618/251‑1 (MQ=255)       tCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGttt                                                                                                                                                                                                                                         <  1:14076/250‑1 (MQ=255)        ccAGTTGAAGTAGATTTTTTTGTTCGGACAACCTCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGTGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAGCTTCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGttttt                                                                                                                                                                                                                                       <  1:453344/251‑1 (MQ=255)         cAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGTGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTc                                                                                                                                                                                                                                      <  1:39254/251‑1 (MQ=255)          aGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCGCGAAGTTGCGGTCGTGGGCACGTTTTTCg                                                                                                                                                                                                                                     <  1:132978/251‑1 (MQ=255)                 gTAGATTTTTTTGGTCGGACTTCCGCTGACGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGTCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGCCCTCGAAGTTGCGTTCGTTTGCACGTTTTTCGAACTCAc                                                                                                                                                                                                                              <  2:259816/251‑1 (MQ=255)                    gATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGc                                                                                                                                                                                                                           >  1:328830/1‑251 (MQ=255)                         ttttGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTcagc                                                                                                                                                                                                                      >  2:219541/1‑251 (MQ=255)                              tACGGACAAACGCTTCTGCACAAAAGCCAACCACAACATTTATACTGATCATTAACAGCAC‑‑TGCCGTCTTCTTCGCGTTTGTAGCTAGAGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGAAGTGTTAGCACAGGCGAGGCCGGGACATCATGAAGGTGTTTTCAAACCACCCGTACATCTCCTTTTGCATGGCCTCGCAGTTGCGGTAAAGGGCACGtttatagcacccacctttcagcagatt                                                                                                                                                                                                                 >  2:259946/1‑227 (MQ=255)                                         gCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTg                                                                                                                                                                                                      >  1:42155/1‑251 (MQ=255)                                           ttATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCCGCAGTTCTTCCCAGTTTGGc                                                                                                                                                                                                    >  2:114490/1‑251 (MQ=255)                                              tGCACAAACGCCAGCCACGGCCTTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGCACATCATTAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCTAACTCACCGCTCCTCAGTTCTTCCCCGTTTGGcccc                                                                                                                                                                                                 >  2:196749/1‑251 (MQ=255)                                                       gccagccaCGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTGCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACt                                                                                                                                                                                        <  1:202611/251‑1 (MQ=255)                                                            ccaNGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTc                                                                                                                                                                                   <  2:406780/251‑1 (MQ=255)                                                                ggCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCAACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGAATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCATCGCTCAGCAGTTCTTCTCAGTTTGGCCCCCAGATCAccttttccctc                                                                                                                                                                               >  1:147256/1‑241 (MQ=255)                                                                                      aGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTg                                                                                                                                                            >  1:125014/1‑250 (MQ=255)                                                                                          caATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGCGATTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTGTGTCCATCCGTTTGCCGTCAATCAGTGAAcgc                                                                                                                                                       >  1:244033/1‑251 (MQ=255)                                                                                              gCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTTGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCCTCCGTTTGCCGGCAATCCGTGAACGCGGGc                                                                                                                                                   >  2:224615/1‑251 (MQ=255)                                                                                              gCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCggg                                                                                                                                                    <  1:114490/250‑1 (MQ=255)                                                                                                     tcttcGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTTCCGTCAATCAGTGAACGCGGGCGAGc                                                                                                                                               <  1:66123/248‑1 (MQ=255)                                                                                                              ttGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAatat                                                                                                                                    <  1:18941/250‑1 (MQ=255)                                                                                                               tGTAGATAGCGCTGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATg                                                                                                                                   <  1:419256/250‑1 (MQ=255)                                                                                                                  aGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTg                                                                                                                               <  2:143264/251‑1 (MQ=255)                                                                                                                   gATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGc                                                                                                                              <  1:32148/251‑1 (MQ=255)                                                                                                                         gccgcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGGGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTcct                                                                                                                        >  2:228266/1‑250 (MQ=255)                                                                                                                            gcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTg                                                                                                                     >  2:105918/1‑251 (MQ=255)                                                                                                                            gcTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTg                                                                                                                     >  1:161405/1‑251 (MQ=255)                                                                                                                                gggCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCt                                                                                                                 <  1:293529/251‑1 (MQ=255)                                                                                                                                               cAGCTGGGATTGAGGCCGTGTTCGCACAGGCGCGGCAGGTACATCATGTAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGCTGCGGTCGTGCGCACGTTTTTTGAACTCACCGCTCAGCAGTTACTCCCCTTTTGGCCCCCAGATCACTTTTTCCATCCGTTTGCCGTCAATCAGTCAACGCGGGCGAGCAGTAGGAATCTGTTTGCCTTCCTGTGCGCTATGCAAATGTTCGtt                                                                                                  >  2:212555/1‑250 (MQ=255)                                                                                                                                               cAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTa                                                                                                  <  1:188775/251‑1 (MQ=255)                                                                                                                                                   tGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAata                                                                                                                                     >  1:37921/1‑212 (MQ=255)                                                                                                                                                   tGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAata                                                                                                                                     <  2:37921/212‑1 (MQ=255)                                                                                                                                                      gATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAg                                                                                           <  2:161405/251‑1 (MQ=255)                                                                                                                                                           aGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGaa                                                                                      <  1:224615/251‑1 (MQ=255)                                                                                                                                                                gtgtTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGGCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGtt                                                                                 <  2:158987/251‑1 (MQ=255)                                                                                                                                                                  gtTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGCGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCg                                                                               >  1:292791/1‑251 (MQ=255)                                                                                                                                                                         caGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAat                                                                        <  2:221639/251‑1 (MQ=255)                                                                                                                                                                                  cAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCt                                                               >  1:319185/1‑251 (MQ=255)                                                                                                                                                                                           catGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTGCGGcacca                                                      <  1:311127/251‑1 (MQ=255)                                                                                                                                                                                                           aaCTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTGCGGCACCACCGGGTTGGCGa                                          >  2:31734/1‑247 (MQ=255)                                                                                                                                                                                                           aaCTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCGATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTGCGGCACCACCGGGTTGGCGAAGAt                                      <  2:244033/251‑1 (MQ=255)                                                                                                                                                                                                                 ccGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGCTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCCGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGCTCGTAGTAATCATCAATCTGCGGCACCACCGGGTTGGCGAAGATTTTGGt                                >  1:46111/1‑251 (MQ=255)                                                                                                                                                                                                                  cGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTGCGGCACCACCGGGTTGGCGAAGATTTTGGtt                               >  1:345003/1‑251 (MQ=255)                                                                                                                                                                                                                        tctcCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCg                                                                                              >  2:165652/1‑182 (MQ=255)                                                                                                                                                                                                                        tctcCTTTTGCATGGCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCg                                                                                              <  1:165652/182‑1 (MQ=255)                                                                                                                                                                                                                        tctacTTTTGCATGGCCTCGAAGTTGCGGCCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCTCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACACGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGGAGTCGAAGGTGAAAGGATCGTAGTGATCATCAATCTGCGGCACCACCGGGTTGGCGACGATTTTGGTTATTACg                         >  2:154156/5‑251 (MQ=255)                                                                                                                                                                                                                        tcncCTTTTGCATGGTCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCAACGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCTTAGTAATCATCAATCTGCGGCACCTCCGGGTTGGCGAAGATTTTGGTTATTAc                          <  2:413175/250‑1 (MQ=255)                                                                                                                                                                                                                                      gCCTCGAAGTTGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTGCGGCACCACCGGGTTGGCGAAGATTTTGGTTATTACGCCCATCTTGTTGcc           <  2:345003/251‑1 (MQ=255)                                                                                                                                                                                                                                               ttGCGGTCGTGGGCACGTTTTTCGAACTCACCGCTCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTGCGGCACCACCGGGTTGGCGAAGATTTTGGTTATTACGCCCATCTTGTTGCCCAGACGCgg  >  2:174118/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  CGCTTTTCCAGTTGAAGTAGATTTTTTTGTACGGACAACCGCTTATGCACAAACGCCAGCCACGGCATTTATCCTGATCA‑‑ATCAGCACAATGCCGTCTTCTTCGCGTTTGTAGATAGCGCCGCTTGGGCAGGTCGCCACGCAGCTGGGATTGAGGCAGTGTTCGCACAGGCGCGGCAGGTACATCATGAAGGTGTTTTCAAACTGCCCGTACATCTCCTTTTGCATGGCCTCGAAGTTGCGGTCGCGGGCACGTTTTTCGAACTCACCGCCCAGCAGTTCTTCCCAGTTTGGCCCCCAGATCACTTTGTCCATCCGTTTGCCGTCAATCAGTGAACGCGGGCGAGCAGTAGGAATATGTTTGCCTTCCGGTGCGCTATGCAAATGTTCGTAGTCGAAGGTGAAAGGTTCGTAGTAATCATCAATCTGCGGCACCACCGGGTTGGCGAAGATTTTGGTTATTACGCCCATCTTGTTGCCCAGACGCGG  >  NC_000913/1538142‑1538628

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