Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,591,418	C→T	E253E (GAG→GAA)	hipA ←	serine/threonine‑protein kinase toxin HipA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,591,418	0	C	T	100.0%	188.7 / NA	56	E253E (GAG→GAA)	hipA	serine/threonine‑protein kinase toxin HipA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/28);  total (28/28)

CTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCAATCGGTAATTTAATGATGTGCGTCGTCGGCGTTATTC  >  NC_000913/1591171‑1591662                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     cTGAACTGATGTCGTAAAACGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCCATCAACCACTGGAAGAACTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCATCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATAtcttcc                                                                                                                                                                                                                                                   >  2:386789/1‑251 (MQ=255)         atgTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCa                                                                                                                                                                                                                                            >  1:177708/1‑250 (MQ=255)                  aaTGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTTCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTTACACATATCTTCCTGTGGCAAGCGaaata                                                                                                                                                                                                                                   >  1:257349/1‑247 (MQ=255)                     ggCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTACATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGCTGCCTGGCCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAAc                                                                                                                                                                                                                               >  1:597196/1‑251 (MQ=255)                      gagTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGACAATCAACCACTGGAAGACCTGGAATTTCAAAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCCAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCCTCTGATTCATATTTCACCGATGAAGGTATACCGACTGTCTGACACATATCTTCCTGTGGCAAGCgccgtacaacc                                                                                                                                                                                                                              >  2:379157/3‑241 (MQ=255)                        gTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGt                                                                                                                                                                                                                            <  1:362964/251‑1 (MQ=255)                          gAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTc                                                                                                                                                                                                                          <  1:460305/251‑1 (MQ=255)                           agtagATAACGGCCGCAAGCCTGAATAAATACGGAGCAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCAGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGCAGGTAAACCGACTGTCTGACACATATCTTCCTGTGACAAGCGAAGTCAAATCGTTCg                                                                                                                                                                                                                         >  2:31564/5‑251 (MQ=255)                           aGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGCAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATCTCTTCCTGTGGCAAGCGAAGTAAAACCGTTCg                                                                                                                                                                                                                         >  2:53796/1‑251 (MQ=255)                            gTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGc                                                                                                                                                                                                                        <  2:40862/251‑1 (MQ=255)                              cGATAAGTGCCGGCAGCCTGAATAACTACGGAGAAGTTTTTTGCATGACCGTCCGTTGAGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCACGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCATCGATGAAGGTAAACCGAATGTCTGACACATATCTTCATGTGGCGAGCGCAGTAAAACCGTTCGCTc                                                                                                                                                                                                                      >  2:484186/1‑251 (MQ=255)                                         gccAGCCTGAATAAATACGGAGAAGTTTGTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAAc                                                                                                                                                                                                           <  1:259111/251‑1 (MQ=255)                                           cAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCa                                                                                                                                                                                                                                            <  2:192357/216‑1 (MQ=255)                                           cAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCa                                                                                                                                                                                                                                            >  1:192357/1‑216 (MQ=255)                                                   ataaataCGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGAGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCaa                                                                                                                                                                                                  >  1:203578/1‑250 (MQ=255)                                                             agaAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGAcc                                                                                                                                                                                       <  1:528423/251‑1 (MQ=255)                                                                    ttttGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCATCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAAc                                                                                                                                                                                 <  1:26351/250‑1 (MQ=255)                                                                           tGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTaaaa                                                                                                                                                                                                                                >  2:149703/1‑196 (MQ=255)                                                                           tGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTaaaa                                                                                                                                                                                                                                <  1:149703/196‑1 (MQ=255)                                                                              ccgtccgtTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACAtt                                                                                                                                                                      <  1:162495/251‑1 (MQ=255)                                                                                 tccgtTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTcc                                                                                                                                                                   >  1:108773/1‑251 (MQ=255)                                                                                           cAATCAACCACTGGACGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATg                                                                                                                                                         >  2:322341/1‑251 (MQ=255)                                                                                               cAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCGCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGAttt                                                                                                                                                     <  2:240948/251‑1 (MQ=255)                                                                                                aaCCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTc                                                                                                                                                    <  2:177708/251‑1 (MQ=255)                                                                                                  ccACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTg                                                                                                                                                  >  1:73092/1‑251 (MQ=255)                                                                                                       ggAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCTAATGTCTGACACAGATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTc                                                                                                                                              <  1:386789/250‑1 (MQ=255)                                                                                                                ggAATTTCATAAAATCATAGCGATCTTTCAGCGCCGCGCTGGAACCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAATCCGTTATCTCAGCATTCCAACGCCTGTCAAAACGTTCTACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACAtt                                                                                                                                     <  2:495449/250‑1 (MQ=255)                                                                                                                            aaTCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATg                                                                                                                                                         <  1:244359/218‑1 (MQ=255)                                                                                                                            aaTCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATg                                                                                                                                                         >  2:244359/1‑218 (MQ=255)                                                                                                                               cATAGCGATCTTTCAGCGACTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTc                                                                                                                     <  2:520457/251‑1 (MQ=255)                                                                                                                                  aGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGc                                                                                                                   <  2:282900/250‑1 (MQ=255)                                                                                                                                                 ccTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTaa                                                                                                   <  1:290322/251‑1 (MQ=255)                                                                                                                                                    cGCTGTACCCCATCAAAAAAGCTATGATCCGCGCGATGCCTGGGCCTCCCTCTGATTCATATTTCACCGATGAAGGTAAACCGAATATCTGACACATATCTTCGTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAGCGCCTGTCAAAACCTTCGACCGCTAACGCGCGCACATTTCCCGCTCTAATGATTTCTGCTTACGGAACATTCAACCCAAGTTCTTTCCCCAGCAGCAGACAGTAATa                                                                                                 <  2:341763/250‑1 (MQ=255)                                                                                                                                                        ggACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCAt                                                                                            >  1:344225/1‑251 (MQ=255)                                                                                                                                                           ccccATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTAt                                                                                         <  2:74583/251‑1 (MQ=255)                                                                                                                                                             ccATCAAAAAAGCCATGATCCGCGCGATGCCTGGGACTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCa                                                                                       <  2:480222/251‑1 (MQ=255)                                                                                                                                                               aTCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTAAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTCACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTACTACTCCTTATCAAc                                                                                     >  2:501945/1‑251 (MQ=255)                                                                                                                                                                  aaaaaaGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCt                                                                                  <  1:192550/251‑1 (MQ=255)                                                                                                                                                                                   cgcgATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCg                                                                 >  2:164610/1‑251 (MQ=255)                                                                                                                                                                                           cTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTg                                                         >  2:266341/1‑251 (MQ=255)                                                                                                                                                                                             gggCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTg                                                         <  1:443252/249‑1 (MQ=255)                                                                                                                                                                                                ccTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTg                                                    >  2:580233/1‑251 (MQ=255)                                                                                                                                                                                                    cATNTGATTCATATTTCCCCGATGAAGTTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTGCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGATAGCGAGCGTCGCATGGGGCGGCCTg                                                <  2:596790/251‑1 (MQ=255)                                                                                                                                                                                                          aTTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCg                                          >  2:173385/1‑251 (MQ=255)                                                                                                                                                                                                                agttCACCGATTAATGTAATCCGAATGTCTGACACATATCTTCCTGTGTCAAGCGAATTAAAACCGTTCGCTCAGCATTACAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCTCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCTATCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCAATc                                    <  2:434135/249‑1 (MQ=255)                                                                                                                                                                                                                 tttCACCGATGAAGGTAAACCGACTGTCTGACACATATATTCCTGTGGCAAGCGAAGAAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAACACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACACAAGTTCTTTCGCCAGCAGCAGACAGTATGACTCATTATCACAGCTTTGGCTGAGATCGAGCAGCGCCTTGGGCTGCCTGATTTAGCAAATCg                                   >  2:661543/1‑251 (MQ=255)                                                                                                                                                                                                                        gATGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCAATCGGTAAttt                            >  2:82954/1‑251 (MQ=255)                                                                                                                                                                                                                         aTGAAGGTAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGAATTGGGCTGCCTGATTTCGCCACTCGGTAAttt                            >  2:545855/1‑250 (MQ=255)                                                                                                                                                                                                                                tAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGt                                                                <  2:410988/207‑1 (MQ=255)                                                                                                                                                                                                                                tAAACCGAATGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGt                                                                >  1:410988/1‑207 (MQ=255)                                                                                                                                                                                                                                        aTGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTg                                                         >  2:4427/1‑206 (MQ=255)                                                                                                                                                                                                                                        aTGTCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTg                                                         <  1:4427/206‑1 (MQ=255)                                                                                                                                                                                                                                           tCTGACACATATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCAATCGGTAATTTAATGATGTGCGTCGTCGGc         >  1:342284/1‑251 (MQ=255)                                                                                                                                                                                                                                             tGACACATATCTTACTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCACTCGGTAATTTAATGATGTGCGTCGTCGGCGt       >  2:498490/1‑251 (MQ=255)                                                                                                                                                                                                                                              gACACAGATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGGCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCGGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCAATCGGTAATTTAATGATGTGCGTCGTCGGc         <  2:257349/248‑1 (MQ=255)                                                                                                                                                                                                                                                 acaTATCTTCCGGTGGCAAGCGAAGTAAAGCAGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCTCATGTCCCGCTTGAATGATTTCTGCGTCCGGAACATCCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCAATCGGTAATTTAATGATGTGCGTCGTCGGCGTTAtt   <  2:216708/251‑1 (MQ=255)                                                                                                                                                                                                                                                  caTATCTTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCAATCGGTAATTTAATGATGTGCGTCGTCGGCGTTATTc  >  2:445293/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     CTGAAATGATGTCGTAAAATGGCGTGAGTCGATAACTGCCGCCAGCCTGAATAAATACGGAGAAGTTTTTTGCATGACCGTCCGTTGCGCCAATCAACCACTGGAAGACCTGGAATTTCATAAAATCATAGCGATCTTTCAGCGCCTCGCTGGACCCCATCAAAAAAGCCATGATCCGCGCGATGCCTGGGCCTCCATCTGATTCATATTTCACCGATGAAGGTAAACCGAATGTCTGACACATATCCTCCTGTGGCAAGCGAAGTAAAACCGTTCGCTCAGCATTCCAACGCCTGTCAAAACGTTCGACCGCTAACGCGCGCACATTTCCCGCTTTAATGATTTCTGCGTCCGGAACATTCAACCCAAGTTCTTTCGCCAGCAGCAGACAGTAATACTCATTATCAACGCTTTGGCTGAGATCGAGCGTCGCATTGGGCTGCCTGATTTCGCCAATCGGTAATTTAATGATGTGCGTCGTCGGCGTTATTC  >  NC_000913/1591171‑1591662

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