Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,286,965	A→G	15.5%	intergenic (+439/+101)	narI → / ← rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,286,965	0	A	G	15.5%	86.8 / 5.0	41	intergenic (+439/+101)	narI/rttR	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (14/20);  new base G (2/5);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGT  >  NC_000913/1286678‑1287250                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                              gAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCCTGTCTCTTCGCCCTTCGTGTCGTCGCCTGAGGCTACGCTCTCTCGCTGAc                                                                                                                                                                                                                                                                                              >  2:541359/1‑289 (MQ=17)     aaCAGGTTTGGAGTATGTAAGGAGGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                                                                                                                                                                                           <  2:32068/288‑1 (MQ=255)           ttttGATTATTTAAGGAAGCTATGGTGGGGGGGGAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                                                                                                                                                                     <  2:903945/288‑1 (MQ=255)           ttttGAGTATTGAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCTCCCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGTTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGcc                                                                                                                                                                                                                                                                                     <  2:234929/288‑1 (MQ=255)            tttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTTGGGTCGTTGCTTGCGTCAACGCCCTCTCGCTCACGCTCGCGCCGa                                                                                                                                                                                                                                                                                   >  2:686115/1‑289 (MQ=255)            tttGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGa                                                                                                                                                                                                                                                                                   >  2:544845/1‑289 (MQ=255)              tGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGGGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                                                                                                                                                                                                                                                  >  2:730669/1‑288 (MQ=255)                    tgtAAGGAAGCGATGGTTGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                                                                                                                                                                                           <  2:321657/287‑1 (MQ=255)                       aaGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                                                                                                                                        <  2:822824/289‑1 (MQ=255)                       aaGGAAGAGGTGGGGGGGGGGGAAGGATTACGCAGCGCTGCGCGTCTCGGCCTTCGGGTAGGTGCCTGCGGGAACGCTCTCTCGCTGGCGCTCAAGGCAAACCTGGGGCGAATCTTATCATCCTTCCCCGCTTGGGCAGAATTTTTGATTGCGGATTCGTTTGAGATTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                                                                                                                                        <  1:817192/289‑1 (MQ=14)                            aGCGATTGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCTAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGAGTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCt                                                                                                                                                                                                                                                                    <  2:865905/288‑1 (MQ=255)                              cGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCACTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATCGCGGATTCGTTCGAGAATTCCGCTGCTTACGAACGTTTTGGTGGTGCAGTAAGTATTACGCAGCGCATGGCGCTTCGCTCTTCGGCACGATGCCTGAGTCACCGCCTTTCCGTTTACGCTCAAGccacacccacttcaacaccttg                                                                                                                                                                                                                                                                 >  1:352591/1‑269 (MQ=9)                                       tggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTc                                                                                                                                                                                                                                                         <  2:41871/288‑1 (MQ=255)                                                 aTTGCTCATCGGTGCTGGCTGCGCCCTCCGGTTGATTGGCTGCGGCAACGTCGTCTGGCGGGTGATAAGGTCAGACCTTGTTGGTAGTTTCTCCCCCTTCCCCGTTTTGGCATAGTATTTAATTGCGTAGTCGTTTGAGATTTCCGTGGCTTATGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATggg                                                                                                                                                                                                                                               <  1:581968/288‑1 (MQ=9)                                                  ttACTCAGCGCTGGGCGCTTCGCCCTTCGGGTCGTGGCCTGCGGAAACGCTCTTTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAAGTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                                                                                                                             <  2:352591/289‑1 (MQ=255)                                                                   gggcgcccgtgggtgcGTGGCCTGCGGAAACGTTTTTTCGGTGGCGCGAAAGTCGAACCTTAGTTGAAGATTGTCATCCGTGTCGGCTGGGTCAGAAGATTTTATGGGGGAGGCGTTTCAGATTTCGGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                                                                                                                                                                            <  1:544845/274‑1 (MQ=14)                                                                    ttCGCCCTTCGGCTTGTTGCCTGCGGCATCGCTCTCTCGCTTGCGCTCTAGTCAAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTTGGCAGAATATTTGATTGCGTATTCGTTTGAGAATTCCGGTTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGAt                                                                                                                                                                                                                           <  2:719098/289‑1 (MQ=14)                                                                           ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                                                                                                                                                                                             <  2:822517/264‑1 (MQ=255)                                                                            tccggTCGTTGCCTTCGGAAACGCCCTCTCGCTGGCGCTCAAGTCGAACCTTGTTCGAAGCTTCTCATCCTTCCCCGCTTGGTCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCgttggg                                                                                                                                                                                                                   <  2:437184/286‑6 (MQ=11)                                                                                     tgcgtgcgGCAACGCATTCTCGCTGGCGATCAAGTCGAACCGTGGTGGTAGCTTGTCATCCTTCCCCGCGTGGGCAGAAGATTTGATTGCGGAGTCGTTTGAGAATTCCGGGGCTTTTGAAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                                                                                                                                          <  1:114545/289‑1 (MQ=17)                                                                                     tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGCTTCGTTTGAGAATTCCGGGGCTGTTGAAAGTGATGGTGGTGGGGGAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGAAACGCTATCTCGTTGACGCTCCAGCCGAAACTTAGTAGAAGCTTCCTATCCTGTCCCGCAGGGGCAGAACATTTggtcgggctctgcttcacactagaat                                                                                                                                                                                                           >  1:822517/1‑263 (MQ=14)                                                                                      gCCTGCGGCACCGCTCTCTGTCTGTCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAggg                                                                                                                                                                                                         <  2:199643/289‑1 (MQ=255)                                                                                       ccTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCgttggga                                                                                                                                                                                                                  <  2:386209/279‑7 (MQ=16)                                                                                                 cGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAACCTTCTCATCCTTCCCCTCATGGGCAGAATATTTAATGTCCGATTCGTTGGGAAGTTCAGGGAcctgtcccct                                                                                                                                                                                              >  1:386209/1‑280 (MQ=11)                                                                                                                 ggcAAGGCGAACCTTGGTGGTAGCTTCTGAGTCTTCCGCGGTTGGGTAGAATTGTTGATAGGGGATTCGTTTGAGATTTTCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                               <  1:730669/286‑1 (MQ=14)                                                                                                                                                        gCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                                                                                                                                        >  1:657561/1‑160 (MQ=255)                                                                                                                                                        gCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                                                                                                                                                                                        <  2:657561/160‑1 (MQ=255)                                                                                                                                                                            tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                               >  2:534950/1‑229 (MQ=255)                                                                                                                                                                            tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTggggg                                                                                                                                                                               <  1:534950/229‑1 (MQ=255)                                                                                                                                                                                gATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGCACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCTGGGCCTTATACTTGTTTGACGCTTCCCTCGAACGTCCAGGGACTTCTCTAGGTTCTTGTGGTGTCGGTAGGATTACTAAGCGCTGCGCGCTCGGACCTTGCGGTCGTTgtgtcccgcaaacctctgccgg                                                                                                                >  1:5886/1‑267 (MQ=2)                                                                                                                                                                                                                tGATGGTGGTGGTGGAAGGATTACACGGCGGTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTTTCATCGGTCCCCGCATGGGCAGAATATGTGATGGCCGATGCGCTTGAGAGTTCAGTGTCTTTTGCAGGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTgc                                                                                                                                                            <  1:116009/212‑1 (MQ=14)                                                                                                                                                                                                                   ttgtggtggGGGAAGGATTACTCGGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGGGCCGAAATTTAGTCGAAGCTTCTCATCCTTCCCCGTATGGGCAGAAGATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAACTTCTCa                                                                             <  2:5886/286‑1 (MQ=255)                                                                                                                                                                                                                          gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGcc                                                                                                                                                  >  2:116009/1‑212 (MQ=255)                                                                                                                                                                                                                          gggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTCGGGCCACGCTCTCTCGCTGGCGCTCGTGTCGGAACTTGGACGAAGACTCCCATCCtttc                                                                      >  2:790857/1‑286 (MQ=14)                                                                                                                                                                                                                                  gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCtt                                                                                                                                                      <  1:669461/200‑1 (MQ=255)                                                                                                                                                                                                                                    ttACTCAGCGGTTCGCGTTTCGCCCTTCGGGTCGTATCCGGCGGCAACCCTCTCTCGCTGGCGTTCGAGCTGAACCTTAGTGAAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                           <  2:323592/289‑1 (MQ=14)                                                                                                                                                                                                                                            cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCgg                                                                                                                                            >  2:669461/1‑200 (MQ=255)                                                                                                                                                                                                                                                                    cGTTGCCCGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGAGTCGAACCTTGGTAGAAGCTTCTCATCCTTCCCCGCATGGGCCGACTCTCTAATTGAGGATTCGTTGGGCAGtt                            >  2:12932/1‑288 (MQ=17)                                                                                                                                                                                                                                                                             cGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACttt                  >  2:692925/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                            gCTGGCGCTCGAGTCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGAGTGGGAGAAGGAATACTCAGCGCTGCGCGCTTCGGTCTTCGGGGCGATGCCTAAGAGAAAGCGCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtgg   <  2:489073/289‑1 (MQ=14)                                                                                                                                                                                                                                                                                             cTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGTGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGTTTCGTTTGGAAGTTCAGGGACTTTTGAAAGTGAtggtggt  >  2:91912/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                              GAGTAACAGGTTTTGATTATTTAAGGAAGCGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGT  >  NC_000913/1286678‑1287250

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑