Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,621,090	C→T	R633R (CGC→CGT)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,621,090	0	C	T	98.3%	102.3 / NA	59	R633R (CGC→CGT)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/30);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCG  >  NC_000913/3620843‑3621334                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     gTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGACGGGCTGAGTCAGTGCCGCGCATCCGACAGCCGTGGACAGTTAATTACCGTGAAAGACACGCAGCGACATGAAACGCGGTATGAATACAAAATCGCCGGTGACCTTACCGCTGTCATTGCCCCGGACCGCAGCAGAAACGGGACACAGTACGCCGCGTGGCGAAAGGCCGTCCGTACCACGCCGGGCGGGCTGACGCGTAGt                                                                                                                                                                                                                                                   >  1:297932/1‑251 (MQ=11) gTTATGAACATGCCCGGTTAGGGAAGATGACGGAGGTGCACCGCGAGGAAGGGCTGAGTCAGTAACGCGCATACGACAGCCGTGGACAGCTAATTGCAGTGAAAGACACGCAGGGACATGAAAAGCGGTATGAATACCAAATCGCCGGTGACGCGAACGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACCGTACGATGCGTCGGGAAAGGCCGTCAGTACCAAGCAGGGCGGGCTGCCGAGCAGt                                                                                                                                                                                                                                                   >  2:261341/1‑251 (MQ=11)   tATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTAt                                                                                                                                                                                                                                                 >  2:180397/1‑251 (MQ=14)        ccatgaccGCTTCGGGCAGATGCCGGCGGTGAACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAAt                                                                                                                                                                                                                                            <  1:72851/251‑1 (MQ=14)              ccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATatgatg                                                                                                                                                                                                                                      >  1:191506/1‑251 (MQ=14)              ccGCTTCGGGCAGATGACGGCGGATCACCGCGAGGAACGGCCGAGTCAGTACCGCGCATACGACAGCCGTGGCCAGTTAATTGCCGTGACAGACACGCACGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAACCGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCCGGGCGGCCTGACGCGAAGTATGGAATACGAtg                                                                                                                                                                                                                                      >  2:177913/1‑251 (MQ=11)                    cGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGcc                                                                                                                                                                                                                                 <  1:180397/250‑1 (MQ=14)                     gggCGGATGACGGGGGGGCACCGCGAGTAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTTGACCGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACAGCGCCGGTGCCCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCgg                                                                                                                                                                                                                               <  2:82375/251‑1 (MQ=11)                        cAGATGACGGCGGTGCGCCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACg                                                                                                                                                                                                                            <  1:287458/251‑1 (MQ=14)                        cAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACg                                                                                                                                                                                                                            >  2:191210/1‑251 (MQ=14)                             gACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTc                                                                                                                                                                                                                        <  1:158895/250‑1 (MQ=14)                             gACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCa                                                                                                                                                                                                                       <  1:3651/251‑1 (MQ=14)                                 gcggTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCg                                                                                                                                                                                                                   <  1:202470/251‑1 (MQ=14)                                     tGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTg                                                                                                                                                                                                               >  2:37085/1‑251 (MQ=14)                                       cACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAc                                                                                                                                                                                                             >  1:276222/1‑251 (MQ=14)                                           gcgAGGAAGGGTGGAGTCAGTACCCCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACTGGTCATCCGCCtttccagt                                                                                                                                                                                                         <  2:49752/251‑8 (MQ=11)                                            ggAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGGGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATGGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGGGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGGCATCCGCCTGACCAGTg                                                                                                                                                                                                        <  1:281886/250‑1 (MQ=9)                                                gAAGGGCTGAGTCAGTACCGCGCAGACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGaaaa                                                                                                                                                                                                    <  2:49095/251‑1 (MQ=14)                                                  aGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGCTGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATTCTGCCGGACGGGTCATCCGCCTGACCCGTGAAAACg                                                                                                                                                                                                  >  2:252279/1‑251 (MQ=9)                                                            cattcccGCGCATCCGACGTCCGTGGACAGTTAATTGCCGGAAAAGACACGCAGGGCCAGGAAACCCGGTATGAATACAACATCGCCGTTGCCCGGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGAGGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCacac                                                                                                                                                                                        <  2:12571/246‑1 (MQ=11)                                                                   gcgcATACGACAGCCTGGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTc                                                                                                                                                                                 <  2:276222/251‑1 (MQ=14)                                                                           gACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGa                                                                                                                                                                         <  1:115833/251‑1 (MQ=14)                                                                            aCAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGAt                                                                                                                                                                        >  2:62172/1‑251 (MQ=14)                                                                            aCAGCCGGGGGCAGTTAATTGCCGTGAAAGACAAGAAGGGCCATGAAACGCGGTATGAATAAAACGTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGAt                                                                                                                                                                        <  2:305629/251‑1 (MQ=11)                                                                            aCAGCCGGGGACAGTTAATGGCCGGGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGAt                                                                                                                                                                        <  2:93890/251‑1 (MQ=9)                                                                                ccGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAAAACCCCTTCCGTTCCGATgcac                                                                                                                                                                    >  2:268179/1‑248 (MQ=255)                                                                                            aaTTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCCCACCACCTTCCGTTACGATGTACTTGACCGGCTgt                                                                                                                                                        >  2:236992/1‑250 (MQ=14)                                                                                             aTTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACCCCACCTTCCGTTACGATGTACTTGACCGGCTg                                                                                                                                                         >  2:60481/1‑249 (MQ=14)                                                                                               tGCCGTGAAAGACACGCAAGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGCTACGATGTACTTGACCAGCTgcta                                                                                                                                                      >  2:213421/1‑247 (MQ=9)                                                                                                         gACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCgg                                                                                                                                           >  2:233199/1‑251 (MQ=14)                                                                                                             cGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCttt                                                                                                                                       >  1:293109/1‑251 (MQ=14)                                                                                                                gggtccATGAAACGCGGTAGGAATACAACGTCGCCGGGGCCCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCTCGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCATCCACACCACCTTCCGTTACGATGTACTGGACCGGCTGATACAGGGAACCGGCTTTGAc                                                                                                                                    <  2:284323/247‑1 (MQ=11)                                                                                                                aGGGCCATGACACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGAc                                                                                                                                    >  1:250460/1‑251 (MQ=14)                                                                                                                aGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGACACGGACGGCCGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTggc                                                                                                                                    >  2:173855/1‑249 (MQ=9)                                                                                                                   gCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGAAGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATACGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGCCggg                                                                                                                                 >  1:128003/1‑250 (MQ=9)                                                                                                                     cATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCCGCTGCTACAGGAAACCGGCTTTGACGGCCg                                                                                                                               >  2:161187/1‑251 (MQ=14)                                                                                                                      aTGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACgg                                                                                                                                  >  2:278429/1‑247 (MQ=14)                                                                                                                      aTGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACgg                                                                                                                                  <  1:278429/247‑1 (MQ=14)                                                                                                                           aCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcacaca                                                                                                                         <  2:219341/251‑1 (MQ=14)                                                                                                                                    gAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATca                                                                                                                <  1:296820/251‑1 (MQ=14)                                                                                                                                         caacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACg                                                                                                           <  1:189107/251‑1 (MQ=14)                                                                                                                                          aacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTa                                                                                                                                                                            >  1:177468/1‑185 (MQ=17)                                                                                                                                          aacaTCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTa                                                                                                                                                                            <  2:177468/185‑1 (MQ=17)                                                                                                                                                 ccGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGAcc                                                                                                   <  1:267242/251‑1 (MQ=14)                                                                                                                                                   ggTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCgg                                                                                                 >  1:274401/1‑251 (MQ=14)                                                                                                                                                                    cATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGc                                                                                 <  1:263748/250‑1 (MQ=14)                                                                                                                                                                    cATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGc                                                                                 <  1:234830/250‑1 (MQ=14)                                                                                                                                                                    cATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATAAAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGc                                                                                 <  2:262418/250‑1 (MQ=14)                                                                                                                                                                               aCGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGCAAGCCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGTCAGCCACACCACCTTCCGTTACGCTGTACTTGACCGCCTGCTACAGGAAACCGGCTTTGACGGCCGCCCACAGCGTTTTCCCCCCGACCTGACCGGCCAACTGATCCGCAGCGAGGATGATGGt                                                                     >  2:77908/1‑251 (MQ=11)                                                                                                                                                                                cGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTc                                                                    >  2:54481/1‑251 (MQ=14)                                                                                                                                                                                       agAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCAc                                                             >  1:25161/1‑251 (MQ=14)                                                                                                                                                                                            cGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACt                                                        >  1:161002/1‑251 (MQ=14)                                                                                                                                                                                              ggACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTgg                                                      <  1:62172/251‑1 (MQ=14)                                                                                                                                                                                                          aTGCGGGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGaa                                          <  2:250460/251‑1 (MQ=14)                                                                                                                                                                                                                           gtccgtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgc                          >  2:227141/1‑250 (MQ=14)                                                                                                                                                                                                                                   cACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCGGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGcaccgcaccg                 <  2:217574/251‑1 (MQ=14)                                                                                                                                                                                                                                      gCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACATTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTg               >  2:113089/1‑250 (MQ=14)                                                                                                                                                                                                                                               gcTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGTCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCTAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGaaa     <  2:23253/251‑1 (MQ=14)                                                                                                                                                                                                                                                  gACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGTCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAATCTTATCCGCATCGAGGATGAGGGTCTGGTCACCCACTTGCACTATGACGAATCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCg  <  1:173855/251‑1 (MQ=2)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     GTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCG  >  NC_000913/3620843‑3621334

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