Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,216,481	(G)6→7	coding (11/765 nt)	nfeF ←	NADPH‑dependent ferric‑chelate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,216,481	1	.	G	88.1%	102.3 / 11.0	42	coding (11/765 nt)	nfeF	NADPH‑dependent ferric‑chelate reductase
Reads supporting (aligned to +/- strand):  ref base . (1/4);  new base G (13/24);  total (14/28)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

TCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGG‑T‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAAC  >  NC_000913/3216232‑3216726                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       tCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑Cg‑gg‑gggg                                                                                                                                                                                                                                                        >  2:92820/1‑250 (MQ=255)    cgcgACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑Cg‑gg‑ggggtgt                                                                                                                                                                                                                                                     <  2:105901/251‑5 (MQ=255)         cGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTc                                                                                                                                                                                                                                               <  2:39116/251‑1 (MQ=255)             gggCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTg                                                                                                                                                                                                                                           >  2:89186/1‑251 (MQ=255)                cGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGtt                                                                                                                                                                                                                                         >  2:319449/1‑250 (MQ=255)                   ggTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCGCAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGTCGCTGATGCGCTCAACGCGTAACACATTCAGTTCACGGAATCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCg                                                                                                                                                                                                                                     <  1:358046/251‑1 (MQ=255)                              ccAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATg                                                                                                                                                                                                                          <  1:125145/251‑1 (MQ=255)                               cAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATgg                                                                                                                                                                                                                         <  2:273619/251‑1 (MQ=255)                               cAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGG‑T‑GTTATTCATTGTTATCGCCTTCGTGATGGt                                                                                                                                                                                                                        <  1:277810/251‑1 (MQ=255)                                aGACGATGCCCTCTTCCGTTACCGTTGGCGGCACCAAGTGAGCGTCAGGTTGAGGAACGAATAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCCTTGCGAACGCGCTGCGGGTAG‑CG‑GGCGGG‑T‑GTTATTCATTGTTATCTCCTTCGTGATGGt                                                                                                                                                                                                                        >  2:250667/1‑251 (MQ=255)                                 gACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCTTTCGTGATGGTa                                                                                                                                                                                                                       <  1:95552/251‑1 (MQ=255)                                       gCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCGTCGTGATGGTAATCAg                                                                                                                                                                                                                  <  1:319340/250‑1 (MQ=255)                                               cctttACCGTTGGCGGCACAAAGTGAGCGTCAGTTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCTGGTAG‑CGTGG‑GGG‑T‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCt                                                                                                                                                                                                          <  1:40793/247‑1 (MQ=255)                                                                  aaaGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgat                                                                                                                                                                                      <  1:186995/251‑1 (MQ=255)                                                                     gTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                                                                                                                                                                                   <  1:266814/251‑1 (MQ=255)                                                                     gTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAaatgataat                                                                                                                                                                                   <  2:62241/251‑1 (MQ=255)                                                                          cGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTg                                                                                                                                                                              <  2:38997/251‑1 (MQ=255)                                                                                          gtagagGTTGCTGTGATCAGT‑AAAGCTACGCGATAGTAAA‑CCGTCCAGCGCACAGACGCCGAGAACAATGCTCTCAAAACCGGCACTGATGCGCTTAAAGCGTAACACAGTCAGTTTATGGCAGCGCATATCATTGCGAACGCGCTTATGGGAGATG‑GG‑GGG‑T‑GTTAATGATTGGTAACACCTTCGTTATTGTAATCAGATATATCTAAAGAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATg                                                                                                                                                              <  2:188630/249‑1 (MQ=255)                                                                                                        gATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTa                                                                                                                                                 <  1:323747/250‑1 (MQ=255)                                                                                                         aTC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTAt                                                                                                                                                                                                                                       >  1:95482/1‑163 (MQ=255)                                                                                                         aTC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTAt                                                                                                                                                                                                                                       <  2:95482/163‑1 (MQ=255)                                                                                                             gTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACt                                                                                                                                            >  1:52976/1‑251 (MQ=255)                                                                                                                   gCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTTAAACCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGGTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCTTTCGTGATGGTAGTCAGATATATCTAAAGAAAGCTGGCAGTTGGTAATGATTGTTAATCATGAGAAATGCAAGCgcttgtgacagcttatatg                                                                                                                                      <  2:356571/251‑19 (MQ=255)                                                                                                                          gAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAg                                                                                                                               <  1:43357/251‑1 (MQ=255)                                                                                                                                  aaCCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATa                                                                                                                        <  2:356315/251‑1 (MQ=255)                                                                                                                                     cGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATa                                                                                                                        >  1:10236/1‑248 (MQ=255)                                                                                                                                                          ggACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTGCGTGATGGTAATCAGATATATCTAAATAAAACGCGCAAATGATAAGGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATGGTCTAATAAGACGAACCGCCTCTTCTCAGGca                                                                                                <  2:323408/251‑1 (MQ=255)                                                                                                                                                              aaTGCGCTGCAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCt                                                                                                                                                                                                                                  >  2:179609/1‑117 (MQ=255)                                                                                                                                                              aaTGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTTCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCt                                                                                                                                                                                                                                  <  1:179609/117‑1 (MQ=255)                                                                                                                                                                         aaCCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGAAAAATGCAAGCGATTTGTAAAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGcc                                                                                 <  1:89186/251‑1 (MQ=255)                                                                                                                                                                                   gATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAAGCAGATATATCTAAAAAAAACTCGCAAATGATAAGGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGGCTATATTCTGAAAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATgcggc                                                                        <  1:374574/250‑1 (MQ=255)                                                                                                                                                                                       cgcTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGa                                                                    >  2:200029/1‑250 (MQ=255)                                                                                                                                                                                          tCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGc                                                                >  2:11504/1‑251 (MQ=255)                                                                                                                                                                                                       cagtcagtTCACGGATTCGCAGATCATTGCGAACGCGCTGCTTGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTACTCAGATATATCTAAATAAAACTCGCAAATGATAATGTTTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGTGTTCGGCGTGAAGGCCTTATCCGGCCt                                                    >  2:257030/1‑250 (MQ=255)                                                                                                                                                                                                           cagtTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtg                                               <  1:79961/251‑1 (MQ=255)                                                                                                                                                                                                           cagtTAACGGAAGCGCAGATCATGGCGAACGCGCTGGGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATGGTTATCGCTTTCGTGATGGTAATCAGATATATCTAAAAAAAACTCGCAAATGATAATGATTGTTGAACATGAAAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtg                                               <  2:260550/251‑1 (MQ=255)                                                                                                                                                                                                               tCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGtgtg                                             >  2:98045/1‑249 (MQ=255)                                                                                                                                                                                                                cACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTgag                                           <  1:111720/250‑1 (MQ=255)                                                                                                                                                                                                                 aCGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                              <  2:170813/246‑1 (MQ=255)                                                                                                                                                                                                                 aCGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACgtgt                                              >  1:170813/1‑246 (MQ=255)                                                                                                                                                                                                                      aGCGCAGATCTTTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGG‑TTGTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTc                                    <  1:212109/251‑1 (MQ=255)                                                                                                                                                                                                                       gcgcAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTcc                                   >  1:347106/1‑251 (MQ=255)                                                                                                                                                                                                                               cATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTa                            >  1:176575/1‑250 (MQ=255)                                                                                                                                                                                                                                         gcgcTGCGGGTAG‑CG‑GG‑GGGGT‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGAAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACGGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCAtttt                 <  1:200029/251‑1 (MQ=255)                                                                                                                                                                                                                                                        g‑gg‑ggg‑T‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAAc  <  1:6836/251‑1 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       TCACGCGACGGTGGGCGTGGTCCTTCCGGCCAGACGATGCCCTCTTCCGTTACCGTTGGCGGCACAAAGTGAGCGTCAGGTTGAGGAAAGAAGAGTTTGCTGTGATC‑GTCAAAGCCACGCGAT‑GTAAAACCGTCCAGCGCCTCGCCGCCGAGGACAATGCGCTGAAAACCGGCGCTGATGCGCTCAACGCGTAACACAGTCAGTTCACGGAAGCGCAGATCATTGCGAACGCGCTGCGGGTAG‑CG‑GG‑GGG‑T‑GTTATTCATTGTTATCGCCTTCGTGATGGTAATCAGATATATCTAAATAAAACTCGCAAATGATAATGATTGTTAATCATGATAAATGCAAGCGATTTGTAGAACTGATATGTCTATAGTCTGATAAGACGAACCGCCTCTTCTCAGGCATCATTACTCAACGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTGTGAGATGAGTCCACGTTTAAAAATCATTTTTACACTTGCAAGAAC  >  NC_000913/3216232‑3216726

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