Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	722,091	C→T	G775E (GGG→GAG)	kdpD ←	sensor histidine kinase KdpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	722,091	0	C	T	97.8%	122.7 / NA	45	G775E (GGG→GAG)	kdpD	sensor histidine kinase KdpD
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (10/34);  total (11/34)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

AAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCCCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGATCTCGCTGGCCTGGCGGGCGTGGGGTGATCCTTCGCTTGCC  >  NC_000913/721861‑722336                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         aaaTTGCAAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGCAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCACCGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCCAAGAGTGGCTCGTCAACGTggctggaagac                                                                                                                                                                                                                                      >  2:370470/1‑242 (MQ=255)      gCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGTCTCGTCAACGTGGATTAAGGT‑CAGCg                                                                                                                                                                                                                                 <  1:332643/250‑1 (MQ=255)       ccAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGt                                                                                                                                                                                                                               <  1:395996/251‑1 (MQ=255)          gTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCt                                                                                                                                                                                                                            <  1:22044/251‑1 (MQ=255)                       ccccGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGGCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGAGTAATCAGCCCCCGTTCAAAGAGTGGCTCTTCAGCGTGGATTAAGGG‑CAGTGGTTCTGGCAGAGAAAGAt                                                                                                                                                                                                               <  2:360836/251‑1 (MQ=255)                         ccGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCTCCGCGTTCTCCAGCAGATTAATCATCGCCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTg                                                                                                                                                                                                             <  1:406969/251‑1 (MQ=255)                                                ccGCGAGCAAACGTATCAAATATCGTTTGCTCCTGGCCTGGCGGAAGACCGGGGCCGGTATCCCAGACATCCAGTTGTAGATTTTCGCCCGCAACGTGGGCATCGCTCCCAATTTCGGCCTGCGCACCCGCATAGTTCACCGCGTTCTCCATCAGATTAAGCAGCACTCGTTGAAAGAGTGGGTCGTCAACGTGGATTAAGGT‑CAGCGGTGCTGGCAGAGAAATATTGATGGGCGTTGATAAACCCGGTTc                                                                                                                                                                                      <  2:49094/251‑1 (MQ=255)                                                    gAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGc                                                                                                                                                                                  <  1:34102/251‑1 (MQ=255)                                                     aGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAATCCCGGTTCCAGCa                                                                                                                                                                                 <  1:319297/251‑1 (MQ=255)                                                                       cGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTCTCCCAGACATCCCGTTGTAGATTTTCGCCCTCACCGTCGGCATCGCTACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCCACGTGGGTTAAAGTACA‑CGGTTCTGGGCGAGAAAGATTGGAGGGGGGCGATATACCCCGTTccggccatttcagggggcggccgc                                                                                                                                                               >  2:408942/1‑229 (MQ=255)                                                                            gctactggcctggcGGAAGACCGGTGCCGTTATCCCAGACCTCCAGTCGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAAGTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGCT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCATCGCGCTGCCGact‑act                                                                                                                                                         <  2:189241/247‑1 (MQ=255)                                                                            gCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGact‑ac                                                                                                                                                          <  1:104796/250‑1 (MQ=255)                                                                               cctggcctggCGGAAGACCGGGGCCGTTATCCCAGGCCTCCAGTTGTAGATTTTCGCTCGCAACGTGGGCATCGATACCAAGTTCGGCCTGCGCACCCGCATATGTCGCCGCGTGCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCACCGTGGATTAGGGT‑CAGCGGTTCTTGCAGAGAAATATTGATGGGCGACGATAAACCCGGTGCCAGGATTTGCATCGCGCTGCCGACT‑Actt                                                                                                                                                        <  2:412385/249‑1 (MQ=255)                                                                                       ggcGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGGTTAAGGG‑CAGCGGGTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCCTTTGCAGAGAGCTGCCGACT‑ACTTCTTCCAGCGt                                                                                                                                              >  1:271249/1‑251 (MQ=255)                                                                                               aCCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCAct                                                                                                                                      <  1:155320/251‑1 (MQ=255)                                                                                                  ggggCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACtctt                                                                                                                                   <  1:62674/251‑1 (MQ=255)                                                                                                      ccGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATCTTTCACAGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGGGGGTTTAGGG‑CCGCGGCTCTGGCAGGGAAAGATTGGTGGGCGAAGGAAAACCCGGGTCCAACATTTGCAACGCGCGTCCGCCTCACTT‑TTCCCCCGTTAACCACcctttctt                                                                                                                               >  2:273638/1‑251 (MQ=255)                                                                                                      ccGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACtctttctt                                                                                                                               <  2:253514/251‑1 (MQ=255)                                                                                                          tATCCCAGACATCCAGTTGTAGATTGTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCaa                                                                                                                            <  2:357550/250‑1 (MQ=255)                                                                                                                      ccAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTGCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAgccgcc                                                                                                                <  1:408942/250‑1 (MQ=255)                                                                                                                       cAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCgg                                                                                                              <  1:406160/251‑1 (MQ=255)                                                                                                                            gTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAATGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACACTTTCTTCAAATTAAAGCCGCCGGACTGa                                                                                                         <  2:118671/251‑1 (MQ=255)                                                                                                                              aGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGaa                                                                                                        <  2:26475/250‑1 (MQ=255)                                                                                                                                 ttttCGCCCTCAACGTGGGCATCGATAACAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAAGT‑CAACGGCTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAACGTTAACCACTCTTTCTTCAAATTAAAACCGCCGGACTGAATTcg                                                                                                    >  1:336257/1‑251 (MQ=255)                                                                                                                                                 gggCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCCACGTGGATTTAGGT‑CAGCGGTTCTGGCAGAGAAAAATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGGCT‑ACTTTTTCCCGCGTTAACCACCCTTTTTTTAAATTAAAGCCGCCGGACTGAAATCGCGCCCTATCCAGtcg                                                                                     >  2:102766/1‑248 (MQ=255)                                                                                                                                                 gggCATCGATACCAATTGCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGTCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGCGCCATATCCAGTAg                                                                                     <  2:157618/250‑1 (MQ=255)                                                                                                                                                  ggCATCGATACCAATTTCGGCCTGCGCACCCGCATATTGCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGGGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGCGCCATATCCAGTAGat                                                                                   <  2:271249/251‑1 (MQ=255)                                                                                                                                                     aTCGATACCAATTTCGTCCGGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGGGACGATAAACCCGGTGCCAGCATTTGCAGCGCTCTTCCGACT‑TCTTTTTCCAGCGTTAACCACTTTATCTTGAAATCAAATCCGCCGGACTGAATTCGCGCCATATCCAGTAGattat                                                                                <  1:182533/251‑1 (MQ=255)                                                                                                                                                                          cgcACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCGCTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCTCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCa                                                          <  1:262782/251‑1 (MQ=255)                                                                                                                                                                                         cACCGCGGTCTCCAGCAGATTAAGCAGCACCCGTGCAAAGAGGGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAGCCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTGCCAGCGTTAACCTCTCTTTCTTCAGGTTAAAGCCGCCGGACTGAATTCGCGGGATATCCTGTAGATTATTCGCC‑AGGCGGGTAGTTTGCAGCACATGCTGACGGa                                           <  2:427563/251‑1 (MQ=255)                                                                                                                                                                                            cgcgTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGAGGCCAGCCAACGTGGGTTAAGGC‑CCGCGGGTCTGGCAGAGGACGATTGCTGGGCGGCGCTAAACCCGGTTCCAGCAATTGCAGCCCGCTGCCAACT‑ACTTCTTCCACCGGTAACCACTCTTTTTTCAAATTAAAGACGCCGGCCTGAATTCGCGCCCTATCCAGTAGATTATTCCCC‑AGTCGGGTAGTGTTCAGCACACGCTGACGGAtct                                        >  1:48563/1‑251 (MQ=255)                                                                                                                                                                                             gcgTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGAtctc                                       <  1:102766/251‑1 (MQ=255)                                                                                                                                                                                                 tctcCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGATCTCGct                                    <  1:235888/250‑1 (MQ=255)                                                                                                                                                                                                    ccAGCATATGAATCAGCACCCGTTCAAAGGGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGG‑CGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGATCTCGctggcct                               <  2:235684/251‑1 (MQ=255)                                                                                                                                                                                                    ccAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTACAAATTAAAGCCGCCGGACTGAAGTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGATCTCGctggcc                                <  1:183626/251‑1 (MQ=255)                                                                                                                                                                                                    ccAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCTCCAGAAAGTCGGTTCAGAT‑TAGCGGCTCTCGGAGGGAAATAATGACGGGCGCCGAAACACCCGGTTTCAACAATTTCAACGGGCTGCCGACT‑ACTTCTTTCAACGTTTACAACTCTTTATTAAAATTAAAGCCGCCCGACTGAACTAGAGCCCTATTCAGTAGATTATTCACCCAG‑CGGGGAGGGTTCACCACATGCTTGCGGGTCTCCctggcc                                >  1:381862/1‑251 (MQ=255)                                                                                                                                                                                                             tAATCAGCACCCGTTCAAAGAGTGGCTCGTCCACGGGGATTAAGGT‑CAGCGGGTCTGGCAGAGAAAAATTGATGGGGGACGATAAACCCCGTTCCAGCATTTGCAGCGCGCTGCCCAAT‑ACTTCTTCCCGCGTTAACCCCCCTTTCTTCAAATTAAAGCCGCCGGACTGAATTTGCGCCATATCCAGGAGATTATTCACC‑CGTCGGGTAGTGTTCAGCACATTCTGACGGATCTCGCTTGCCCggcggggg                       >  1:369624/1‑249 (MQ=255)                                                                                                                                                                                                             tAATCAGCACCCGTTCAAAGAGTGGCTCCTCAACGTGGATTAAGGT‑CAACGGTTCTTGCAGAGAAAGATTGATGGGCGAAGATAAAACCGGTTACAACAATTGCAGCGCGCTTCCGAAC‑AATTCTTCCAACGTTAACCACCCTTTCTTCAAATTAAAGCCGCCCGAATGAATTCGCGCCATATCCAGTAGATTATTCACA‑AGGCGGGGCGTGTTCAGAACATACCGACGGATATCGCTGGCCTggcgggcg                       >  1:369652/1‑251 (MQ=255)                                                                                                                                                                                                               aTCAGGACCCGTTCAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CATCCGTACTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCATCATTTGCAGCGCGCTGCCGACT‑ACTGCTTCCAGCGTTAACCCCTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGTGCCATACCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGATCGCGCTGGCCTGGCGGGCGTg                     <  2:130958/251‑1 (MQ=255)                                                                                                                                                                                                                       ccGTTCAAAGAGTGGCTCGTGAACGCGGATGAAGGT‑CAGCGGTTCAGGCAGAGAAAGATTGAGGGGCGACGATAAACCCGGTTCCAGCTTTTGCAGCGCGCTGCCGACT‑ACTTCGTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAAGTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCATCACATGCTGACGGATCTCGCTGGCCTGGCGGGCGTGGGGTGATc             <  2:19017/251‑1 (MQ=255)                                                                                                                                                                                                                            cAAAGATTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTTATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCCCGCTGCCGACT‑CCTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGCCTGAATTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGATCTCGCTGGCCTGGCGGGCGTGGGGTGATCCTTc         <  1:78237/250‑1 (MQ=255)                                                                                                                                                                                                                            cAAAGAGTGGCTCGTCCACGTGGATTAAGGG‑CAGCGGGTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACtctttctt                                                                                                                               >  1:147128/1‑133 (MQ=255)                                                                                                                                                                                                                            cAAAGAGTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACtctttctt                                                                                                                               <  2:147128/133‑1 (MQ=255)                                                                                                                                                                                                                                gagTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGAGTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACACTTTCGTCAAAGTAAAGCCGCCGGACTGAATTCGCGCCATAGCCAGGAGATTATTCACC‑AGGCGGGTAGGGTTCAGCACATGCTGACGGATCTCGCTGGCCTGGCGTGCGTGGGGTGATCCTGCGCTTg    <  2:381862/251‑1 (MQ=255)                                                                                                                                                                                                                                  gTGGCTCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGATCTCGCTGGCCTGGCGGGCGTGGGGTGATCCTTCGCTTGcc  <  1:104300/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         AAATTGCCAGTCCAAGCCCTACCCCCGGTACTGCCGACTCTTTATTCCCGCGAGCAAACTTATCAAATATCGTCTGCTCCTGGCCTGGCGGAAGACCGGGGCCGTTATCCCAGACATCCAGTTGTAGATTTTCGCCCTCAACGTGGGCATCGATACCAATTTCGGCCTGCGCACCCGCATATTTCACCGCGTTCTCCAGCAGATTAATCAGCACCCGTTCAAAGAGTGGCCCGTCAACGTGGATTAAGGT‑CAGCGGTTCTGGCAGAGAAAGATTGATGGGCGACGATAAACCCGGTTCCAGCATTTGCAGCGCGCTGCCGACT‑ACTTCTTCCAGCGTTAACCACTCTTTCTTCAAATTAAAGCCGCCGGACTGAATTCGCGCCATATCCAGTAGATTATTCACC‑AGTCGGGTAGTGTTCAGCACATGCTGACGGATCTCGCTGGCCTGGCGGGCGTGGGGTGATCCTTCGCTTGCC  >  NC_000913/721861‑722336

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