Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,179,383	G→A	K423K (AAG→AAA)	tolC →	outer membrane channel TolC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,179,383	0	G	A	100.0%	175.9 / NA	51	K423K (AAG→AAA)	tolC	outer membrane channel TolC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/32);  total (19/32)

CAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTG  >  NC_000913/3179153‑3179625                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   cAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTGCAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGa                                                                                                                                                                                                                                 >  2:299447/1‑250 (MQ=255)     ccGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACg                                                                                                                                                                                                                              >  2:335625/1‑249 (MQ=255)       gTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGc                                                                                                                                                                                                                           >  2:314612/1‑250 (MQ=255)              cctcctTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAAGCAGCTGAATAGTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATct                                                                                                                                                                                                                    <  2:166073/250‑1 (MQ=255)                 cctTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATctgct                                                                                                                                                                                                                 <  2:161486/250‑1 (MQ=255)                      aacaacATTAAAGCATCTATCAGTAGCATGAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACt                                                                                                                                                                                                           <  1:94588/251‑1 (MQ=255)                              aaaTGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATg                                                                                                                                                                                                   <  1:279171/250‑1 (MQ=255)                                aTGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATgcg                                                                                                                                                                                                 <  2:8155/251‑1 (MQ=255)                                   cATCAATCAGTAGCATTAACGCCAACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAAGCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGGTCTGCTGTCACTGAACAATGCGCTg                                                                                                                                                                                              <  2:150453/251‑1 (MQ=255)                                           aGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAAcc                                                                                                                                                                                      <  2:352630/251‑1 (MQ=255)                                              aGCATTAACGCCTAGAAACAAACAGTCGGTTACGATCATAGCTCATTTGACGCGAAGGAAGCGGGCTCCTCGGTAGGTACGCGCACCATTGTTGATGTGTTGGCTGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTACTCAGCTGAATATCAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGTCCTAAACAATGCGATGAGCAAACAGGt                                                                                                                                                                                   >  2:68472/1‑251 (MQ=255)                                              aGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAACCCGGt                                                                                                                                                                                   >  1:104535/1‑251 (MQ=255)                                               gCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGtt                                                                                                                                                                                  >  1:205678/1‑251 (MQ=255)                                                    aaggCCTACAAACAAGCCGTAGTTTCCGCTCAAATCTCATTAGTCGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTGGATGTGTTGGATGCGACCGCCCCGTTGTAAAACGCCAATCAAGGGCTGGCTAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCATCTCTGGGTACGTTGAGCGGGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCAc                                                                                                                                                                             <  1:25387/248‑1 (MQ=255)                                                    aaCGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCa                                                                                                                                                                              <  1:335625/250‑1 (MQ=255)                                                       gCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTa                                                                                                                                                                           >  2:235833/1‑250 (MQ=255)                                                          cacaaacaaGCCGGAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGGGTTGGATGCGAACACCACGTTGTACAACGCCAAGCAAGAGCTGTCTAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCGGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATcc                                                                                                                                                                       <  1:264876/250‑1 (MQ=255)                                                           acaaacaaGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACACCGCCACGCACGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTACTCCg                                                                                                                                                                      >  2:395464/1‑251 (MQ=255)                                                           acaaacaaGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATc                                                                                                                                                                        <  2:201555/249‑1 (MQ=255)                                                                    ccGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGtt                                                                                                                                                             <  1:318555/251‑1 (MQ=255)                                                                         gTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCAcc                                                                                                                                                        <  1:67030/251‑1 (MQ=255)                                                                                cTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACg                                                                                                                                                 >  1:100151/1‑251 (MQ=255)                                                                                         cATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAAc                                                                                                                                          <  1:1081/249‑1 (MQ=255)                                                                                           ttAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGa                                                                                                                                       >  2:287050/1‑250 (MQ=255)                                                                                                 gcgATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTAt                                                                                                                                <  1:16178/251‑1 (MQ=255)                                                                                                     tGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAACCGCCGGAACAGAATGCTATTGCt                                                                                                                            >  1:436385/1‑251 (MQ=255)                                                                                                             ggCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTa                                                                                                                    <  2:4001/251‑1 (MQ=255)                                                                                                                        cggtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGACTGCTATTGCTGATGGTTATGCGCCTGAt                                                                                                          >  2:292322/1‑250 (MQ=255)                                                                                                                         ggtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCGGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg                                                                                                        <  2:371645/251‑1 (MQ=255)                                                                                                                          gtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGGTATGCGCCTGATAGc                                                                                                       >  2:241661/1‑251 (MQ=255)                                                                                                                                gTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCa                                                                                                 <  1:151957/251‑1 (MQ=255)                                                                                                                                           gATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCa                                                                                      <  1:373627/251‑1 (MQ=255)                                                                                                                                                          accaccCCGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCATCTCTGGGTACGTTGAACGAGCAGGATCTGCGGGCACTGAGCAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGa                                                                                                                                       <  1:61331/187‑1 (MQ=255)                                                                                                                                                          accaccACGTTGTACAACGACAAGCAAGAGCTGGAGAATGCGCGTTCTAACTACCTGCTTAATCAGCTGACTATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTACACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGa                                                                                                                                       >  2:61331/1‑187 (MQ=255)                                                                                                                                                                acGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTac                                                                 <  2:310418/251‑1 (MQ=255)                                                                                                                                                                acGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTac                                                                 <  1:235833/251‑1 (MQ=255)                                                                                                                                                                       aCAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGt                                                          <  1:325868/251‑1 (MQ=255)                                                                                                                                                                        cAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTa                                                         <  1:299447/251‑1 (MQ=255)                                                                                                                                                                          aCGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAAc                                                       <  1:218292/251‑1 (MQ=255)                                                                                                                                                                                    agagCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCGCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAAcc                                             <  2:436385/251‑1 (MQ=255)                                                                                                                                                                                                  cgcgTTATAACTACCTGATTAATAAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTACTCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGACTGATAGCCCGGCACCAGTCGTTCAGCAAAAATCCGCACGCACTACCACCAGTAACGGCCATCACCCTTTCCGTaaatgt                               >  2:286612/1‑247 (MQ=255)                                                                                                                                                                                                        ataACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGCAACGGTCATAACCCTTTCCGTAACTGATgacg                          >  2:350914/1‑250 (MQ=255)                                                                                                                                                                                                          aactaactGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGGTGAACGGGCAGGATCTGCTGGCACTGAACAATGCGCTTAGCAAACCGGTTTCCACTAATCCGGGAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCATGCACTACCTCCAGTAACGGTCATAACCCTTTCCGTAACTGATgactac                        <  2:438198/250‑4 (MQ=255)                                                                                                                                                                                                              aCCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACgg                   >  2:58342/1‑251 (MQ=255)                                                                                                                                                                                                               ccTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACggg                  <  1:290773/251‑1 (MQ=255)                                                                                                                                                                                                               ccTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGAACTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCCCCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACggg                  <  2:228430/251‑1 (MQ=255)                                                                                                                                                                                                                  gATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTACACTAATCCGGCAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGACATAACCACTTCCGTAACTGATGACGACGACGGGGCt               >  2:132231/1‑251 (MQ=255)                                                                                                                                                                                                                      aaTCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCg            <  1:314612/250‑1 (MQ=255)                                                                                                                                                                                                                              gAATATTAAAGCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGt     <  2:233473/249‑1 (MQ=255)                                                                                                                                                                                                                               aaTATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTg  >  1:413551/1‑251 (MQ=255)                                                                                                                                                                                                                               aaTATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTg  <  1:287050/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   CAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTG  >  NC_000913/3179153‑3179625

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