Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	621,495	C→T	S232N (AGT→AAT)	fepD ←	ferric enterobactin ABC transporter membrane subunit FebD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	621,495	0	C	T	98.0%	163.6 / NA	49	S232N (AGT→AAT)	fepD	ferric enterobactin ABC transporter membrane subunit FebD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/18);  minor base G (0/1);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAA  >  NC_000913/621206‑621789                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                       aCGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCNCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCNACTGTGTGCGCGCCACGCGATTGCCCAgcgc                                                                                                                                                                                                                                                                                              <  2:489997/300‑1 (MQ=255)      gCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCAACCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTcg                                                                                                                                                                                                                                                                                         >  1:451204/1‑300 (MQ=255)           aGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGTACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGACAACTATTGCCGTCGCACTACCAAAAAGCACGGCAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGTCGTCTCGgtg                                                                                                                                                                                                                                                                                    >  2:365127/1‑300 (MQ=255)           aGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGACAGCAGACCAATCCACTGTGTGTGCGACACGCGATCGCCCAGCGCCGTCGCGgtg                                                                                                                                                                                                                                                                                    >  2:276419/1‑300 (MQ=255)                gATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCt                                                                                                                                                                                                                                                                               >  2:409216/1‑300 (MQ=255)                 aTCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTg                                                                                                                                                                                                                                                                              >  2:415251/1‑300 (MQ=255)                    aGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCg                                                                                                                                                                                                                                                                           <  1:299994/300‑1 (MQ=255)                     gCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCGCGCGATTGCACAGCGCCGTCGCGGTGTCGCTGcaga                                                                                                                                                                                                                                                                          >  1:89146/1‑297 (MQ=255)                              cACCAATAAACACACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGCGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAAc                                                                                                                                                                                                                                                                 <  2:418183/300‑1 (MQ=255)                                ccAATCAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGCAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAACCTg                                                                                                                                                                                                                                                               >  2:125322/1‑300 (MQ=255)                                ccAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATNGGCACCCACCAGCCAACGCGCCATATNCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTg                                                                                                                                                                                                                                                               >  2:455887/1‑300 (MQ=255)                                              gACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACt                                                                                                                                                                                                                                                 >  1:687719/1‑300 (MQ=255)                                                     gAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAAt                                                                                                                                                                                                                                          >  1:67743/1‑300 (MQ=255)                                                                      gggCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAANTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCNGCCGCGGCTCAAACTGTTCAGCGCGCGACTCAGTAGTAGCGCAGTTGCTCCGGc                                                                                                                                                                                                                         >  2:644099/1‑300 (MQ=255)                                                                           caATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCAAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGATGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGAtta                                                                                                                                                                                                                    >  2:199512/1‑298 (MQ=255)                                                                               cACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCGCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCAc                                                                                                                                                                                                                >  2:652899/1‑300 (MQ=255)                                                                                aCCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCAcc                                                                                                                                                                                                               <  1:389663/300‑1 (MQ=255)                                                                                      ccgATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATc                                                                                                                                                                                                         <  1:296349/300‑1 (MQ=255)                                                                                                      aCAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAgg                                                                                                                                                                                         <  1:191455/300‑1 (MQ=255)                                                                                                       cagcagcagGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGt                                                                                                                                                                                        >  1:439871/1‑300 (MQ=255)                                                                                                               gggcagggGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGTAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGATGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCGGCACAGGGATCAGCACCATTTTTAAGGTATGTAGAt                                                                                                                                                                                >  2:73623/1‑300 (MQ=255)                                                                                                               gggcagggGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGAt                                                                                                                                                                                >  1:225904/1‑300 (MQ=255)                                                                                                                       gTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGATGCCGAGGCTTAAACTGTTCAGCGCGCGACGCAGTAGTAGTGCAGTTGCTCCGGCGGTCAGCACCGAGATCAGCACCACTTTTAGGGTATGTAGATTGCGAata                                                                                                                                                                        >  2:264722/1‑300 (MQ=255)                                                                                                                                               gACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGaaa                                                                                                                                                >  1:690590/1‑300 (MQ=255)                                                                                                                                                    gCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATTAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCaa                                                                                                                                           >  2:103770/1‑300 (MQ=255)                                                                                                                                                              gCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGTGCCAACTATTGCCTTCGCATTACCACAAAGCACGGTGATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCTACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCC‑‑AAACGCAACTGATCGTAGAc                                                                                                                               <  2:266063/300‑1 (MQ=255)                                                                                                                                                              gCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCTAACCGGCTTGCCAGAAACGCAACTGATCGTAg                                                                                                                                 >  1:301939/1‑300 (MQ=255)                                                                                                                                                                 cccaccAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAATACAATCGGGCCAACTATTGCCGTCGCATTACCACAAAGCACGGTAATCGCCAGCAGACCAAGCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCTACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACg                                                                                                                              <  2:654167/300‑1 (MQ=255)                                                                                                                                                                        gCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTAAGGAt                                                                                                                       >  2:621219/1‑300 (MQ=255)                                                                                                                                                                             cgcgCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGc                                                                                                                  <  1:551049/300‑1 (MQ=255)                                                                                                                                                                               cgcCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGAGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGAATCAGTAATAGCGCCGTTGCTCCGGCGATCAGCACCGGGATCAGCACAACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGAAGTCAGTATTAAGCAg                                                                                                                >  1:57076/1‑300 (MQ=255)                                                                                                                                                                                  cATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGc                                                                                                             <  2:581925/300‑1 (MQ=255)                                                                                                                                                                                                    gCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAGAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGNTCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCANGATTAAGCAGGGCGATGCCGCTGGTCAGTcc                                                                                           <  2:500908/300‑1 (MQ=255)                                                                                                                                                                                                        aTAAAGGCAATCGTGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCTGCGATCAGCACCTGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGTCGATGCCGCTGGTCAGTCCTTcc                                                                                       <  2:88842/300‑1 (MQ=255)                                                                                                                                                                                                           aaGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGc                                                                                    <  2:293337/300‑1 (MQ=255)                                                                                                                                                                                                             ggcaaccgtgaaaaCTAGTGCCGTCGCACTACCACCAAGCACGGTAATCGCCAGCAGATCAATCAACTGTGTGCGCGCCACGCGAGTGCCCAGCGCCGTAGCGGGGTCGCTGCCGAGGCTCAAACTGTTCAGCACGCGACTCAGTAATAGCGCAGTTGCTCCTGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCAc                                                                                  <  2:222582/288‑1 (MQ=255)                                                                                                                                                                                                             gacAATCGTGCTAACTATGGCCGTCGCACTACCACAAAGCACGGTGATCGCCATCAGACCAATCACCTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCTGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCAc                                                                                  <  2:225147/298‑1 (MQ=255)                                                                                                                                                                                                                                         cTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGANCCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCAAGGGtt                                                      >  2:601304/1‑300 (MQ=255)                                                                                                                                                                                                                                         cTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGAATCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCCGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGtt                                                      >  2:383680/1‑300 (MQ=255)                                                                                                                                                                                                                                             cacaAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAAc                                                  <  2:666031/300‑1 (MQ=255)                                                                                                                                                                                                                                               caAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCANCACCACTTTTAAGGTATGTAGATTGCNAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGc                                                <  2:468837/300‑1 (MQ=255)                                                                                                                                                                                                                                                      cGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGAc                                         >  1:584514/1‑300 (MQ=255)                                                                                                                                                                                                                                                           aTCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAAc                                    <  2:21773/300‑1 (MQ=255)                                                                                                                                                                                                                                                               ccAGCAGACCAATCAACTGTGTGTGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCACTGCCGCGGTTCAAACTGTTCAGCGCGCGAATCAGTAGTAGCGCAGTTGCTCCGGTGATCAGCACTGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATCTCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGAATTAAGCAGGGCGATGCCGCTGGTCGGTCCTTCCAGCACCGCCGCCAGCGCAACGCCCCCCAGGGTTAAACGCACCGGACTTCACTGcc                                >  2:250005/1‑300 (MQ=255)                                                                                                                                                                                                                                                                       ccAATCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGATGCCGAGGCTCAAACTGTTCAGCGCGCGAGTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAGCGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGACCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGccc                        >  2:266779/1‑300 (MQ=255)                                                                                                                                                                                                                                                                         actCAACTGTGTGCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGATGCCGAGGCTCAAACTGTTCAGCGCGCGGATCAGTAATAGCGCAGTTGCTCCGGCGATCAGCAACGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTCGACGTCAGGATTAAGCAGGGCGGTGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTg                      >  2:260657/3‑300 (MQ=255)                                                                                                                                                                                                                                                                                 gtgtgCGCGCCACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCTCTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCgg              >  1:241401/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                          ccACGCGATTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGcaa     >  2:624498/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                             cgcgATTGCCCAGCACCTTCTCGGTTTCGCTGCCTAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCATAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGcaacaa  <  2:78399/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                                       ACGTCGCACGAGGAAGATCAGCACCGGTGCACCAATAAACGCACTGACCACAGAAACGCGCAGTTCGCCGGGCACAATCACCCGCCCGATGATATCGGCAAACAGCAGCAGGGCAGGGGTAGCAAGTAGCGTGACGGGCAGCGACCAGCGATGATCGGCACCCACCAGCCAACGCGCCATATGCGGCATCATCAGGCCAATAAAGGCAATCGGGCCAACTATTGCCGTCGCACTACCACAAAGCACGGTAATCGCCAGCAGACCAATCAACTGTGTGCGCGCCACGCGACTGCCCAGCGCCGTCGCGGTGTCGCTGCCGAGGCTCAAACTGTTCAGCGCGCGACTCAGTAATAGCGCAGTTGCTCCGGCGATCAGCACCGGGATCAGCACCACTTTTAAGGTATGTAGATTGCGAATATCCAGCGAACCGGCTTGCCAGAAACGCAACTGATCGTAGACGTCAGGATTAAGCAGGGCGATGCCGCTGGTCAGTCCTTCCAGCACCGCCGCCAGCGCCACGCCCGCCAGGGTTAAACGCACCGGACTTAACTGCCCGCCGCCCTGACTGCCGGTAAAGGCAACAA  >  NC_000913/621206‑621789

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