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%	193.7 / NA	58	K423K (AAG→AAA)	tolC	outer membrane channel TolC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/33);  total (25/33)

AGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGAACG  >  NC_000913/3179144‑3179629                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       aGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTggg                                                                                                                                                                                                                                             >  2:303572/1‑251 (MQ=255)          cAGACCGCGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGa                                                                                                                                                                                                                                     >  1:358258/1‑250 (MQ=255)               cGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTCACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCATGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGACCGAGc                                                                                                                                                                                                                               >  2:391124/1‑251 (MQ=255)                gTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCa                                                                                                                                                                                                                              <  1:217989/251‑1 (MQ=255)                    gTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGAt                                                                                                                                                                                                                          <  2:112322/251‑1 (MQ=255)                      tcctccTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATc                                                                                                                                                                                                                         <  1:375699/250‑1 (MQ=255)                       cctcctTCAACAACATTAATGCATCTATCAGTAGCGTTAACGCCTGCAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAGCGAGCAGGATct                                                                                                                                                                                                                        <  2:4260/250‑1 (MQ=255)                          cctTCAACAACATTAATGCATCTAGCAGGAGCATTGACGCCTACAAACAAGCCGTAGTTTCCGCGCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCATGATctgctg                                                                                                                                                                                                                    <  2:61640/251‑1 (MQ=255)                           ctTCAACACCGTTAATGCATCTATCAGCAGAATTACCGCCTACAAACAAGCCGTCTTTTCCGCTCAACGCTCATTAGAAGCGCTGGAAGCGGGCTAATCGGTCGGTACGCGTACCATTGTTGATGTGTTGGCTGCGACCACCACGTTGTACAACACCAAGCACGAGCTGGAGCATGAGCGTTATAACTACCTGATTAATAAGCTGAAAATTAAATAAGATCTGGGTACGTTGAACGAGCAGGATCTGCTgg                                                                                                                                                                                                                   >  2:296625/1‑251 (MQ=255)                                         aTGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATgcg                                                                                                                                                                                                     >  1:240930/1‑251 (MQ=255)                                              tCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGAAGCGATGGAAGCGGGCTACTAGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTg                                                                                                                                                                                                  >  2:293648/1‑249 (MQ=255)                                                 aTCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGAGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGAATAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGCGCAGGATCTGCTGGCACTGAACAATGAGCttagaac                                                                                                                                                                                             >  2:120878/1‑245 (MQ=255)                                                 aTCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCACAAAGCTTATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTTCCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGTCACTGAACAATGCGCTGAGCaa                                                                                                                                                                                             <  2:372867/251‑1 (MQ=255)                                                    aGTAGCATTAACGCCTACAAACAAGCCATAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTGGGATGCGAACACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAAcc                                                                                                                                                                                          <  1:240947/251‑1 (MQ=255)                                                    aGTAGCATTAACGCCTAAAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACAATTGTTGATGTGTTGGATGCGCCCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTCACTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGACCGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCaaaac                                                                                                                                                                                          >  2:158328/1‑249 (MQ=255)                                                          aTTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGACGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATCTTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAACCCGGTTTc                                                                                                                                                                                    >  1:409888/1‑251 (MQ=255)                                                               cGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCgg                                                                                                                                                                                        >  2:39867/1‑242 (MQ=255)                                                               cGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCgg                                                                                                                                                                                        <  1:39867/242‑1 (MQ=255)                                                                 ccTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAAt                                                                                                                                                                             >  2:132017/1‑251 (MQ=255)                                                                    acaaacaaGCCGTAGTTTCCGCTCAAAGCTCATTAGACTCGATGGAATCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAGCTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCg                                                                                                                                                                          <  2:53446/251‑1 (MQ=255)                                                                      aaacaaGCCGAAGAATCAGCTCCAGGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGAACCATCGTTGATGTGTTGGCTGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCCGCTCTGGGCACGTTGAACGAGCAGTATCTGCTGGCACTGAACAAGGAGCTGAGCAAACCGATTTCCACTACTCCgg                                                                                                                                                                         >  2:82648/1‑250 (MQ=255)                                                                             ccGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCACGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGa                                                                                                                                                                        >  2:240597/1‑244 (MQ=255)                                                                             ccGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGa                                                                                                                                                                        <  1:240597/244‑1 (MQ=255)                                                                                   tttCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCg                                                                                                                                                           >  2:26205/1‑251 (MQ=255)                                                                                              aGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCgg                                                                                                                                                 >  1:371689/1‑250 (MQ=255)                                                                                                    ttAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGa                                                                                                                                           >  1:203805/1‑250 (MQ=255)                                                                                                            gATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGCGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTg                                                                                                                                  <  1:88935/251‑1 (MQ=255)                                                                                                                          aCTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCTCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATgc                                                                                                                     <  2:226088/250‑1 (MQ=255)                                                                                                                           cTCGGTCGGTACGCGTACCATTGTTGATGGGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATgcgc                                                                                                                   >  2:117924/1‑251 (MQ=255)                                                                                                                                  ggtACGCGTACCATTGTTGATGTGTTGTATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCATCTCTGTGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAAGGCGCTGAGCAAACCGGTTTCCTCTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg                                                                                                            <  1:45917/251‑1 (MQ=255)                                                                                                                                   gtACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGc                                                                                                           >  2:108932/1‑251 (MQ=255)                                                                                                                                        cgTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGGTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGc                                                                                                      <  2:120839/251‑1 (MQ=255)                                                                                                                                             cATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAg                                                                                                 <  1:429863/251‑1 (MQ=255)                                                                                                                                                          ttGGATGCGACCACCCAGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGGGCAGGATCTGCTGGCACTGAACAATGCGCTGGGCAAAACGGTTTCCACTAATCCGGAAAACGTTGCACAGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCATGATAGCCCGGCAcacactggtccacaacc                                                                                    >  2:392014/1‑235 (MQ=255)                                                                                                                                                           tGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAACCAGCTCTGGGTACGGTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCCCTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTAGGCGCCTGATAGCCCGGCACCATGCGTTCAGctaaa                                                                                    <  2:91847/250‑5 (MQ=255)                                                                                                                                                            ggATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACAt                                                                                  >  2:323630/1‑251 (MQ=255)                                                                                                                                                                  gACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCcgcac                                                                            <  1:303572/251‑1 (MQ=255)                                                                                                                                                                  gACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGCTCAGCAAACATCcgcac                                                                            <  1:163018/251‑1 (MQ=255)                                                                                                                                                                     caccacGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAACCATCcgcacgc                                                                          >  2:4681/1‑250 (MQ=255)                                                                                                                                                                      accacGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCgcacgcac                                                                        >  1:213287/1‑251 (MQ=255)                                                                                                                                                                           gTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTacca                                                                   <  1:323630/251‑1 (MQ=255)                                                                                                                                                                              gTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTaccacc                                                                 <  2:358258/250‑1 (MQ=255)                                                                                                                                                                                 cAACGCTAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGAATAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTa                                                             <  2:357207/251‑1 (MQ=255)                                                                                                                                                                                 cAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCGCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACAGCCGCACGCACTGCCACCAGTa                                                             <  1:158328/251‑1 (MQ=255)                                                                                                                                                                                          gcaagAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCATGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATa                                                    <  1:314473/251‑1 (MQ=255)                                                                                                                                                                                          gcaagAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATa                                                    <  1:132017/251‑1 (MQ=255)                                                                                                                                                                                                   ggCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTcc                                           <  1:85632/251‑1 (MQ=255)                                                                                                                                                                                                   ggCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAAGCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACGAATCCGGAAAACGTTGCACCGCAAACGCCGTAACAGAATGCTATTGCTGATGGTGATGCGCCTGATAGCCCGGCACCATTCGGTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCttt                                             <  1:273780/249‑1 (MQ=255)                                                                                                                                                                                                               ttATAACTACGTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGCACGAGCAGGATCTGCTGGCAATGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAACCGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACAAGTAACGGTCATAACCCTTTCCGTAACTGATgac                               >  2:432940/1‑251 (MQ=255)                                                                                                                                                                                                                   aaCTACCTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTACCTGCTgacgacg                           >  2:141372/1‑251 (MQ=255)                                                                                                                                                                                                                        ccTGATTAATCAGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACggg                      <  1:85337/251‑1 (MQ=255)                                                                                                                                                                                                                                   agcggcatatTAAATCAGATCGGGGTACGTTGACCGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAAACGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTACGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGcgcc           <  2:358310/245‑4 (MQ=255)                                                                                                                                                                                                                                   aGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGcccc           <  1:108932/251‑1 (MQ=255)                                                                                                                                                                                                                                   aGCTGAATATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGcccc           <  1:293648/251‑1 (MQ=255)                                                                                                                                                                                                                                      cgAATAGTAAATCATCTCTGGGTACGTTGAACGAGCATGAGCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGGTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACAGCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGTCCCGTc        <  2:203805/250‑1 (MQ=255)                                                                                                                                                                                                                                        aaTATTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTg      <  1:221110/251‑1 (MQ=255)                                                                                                                                                                                                                                           atTAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGAAc   <  1:26205/251‑1 (MQ=255)                                                                                                                                                                                                                                            ttAAATCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGTAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCATCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGAAACTGATGACGACGACGGGGCTTCGGCCCCGTCTAAACg  >  1:5390/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       AGCGTCGTGCAGACCGTGCGTTCCTCCTTCAACAACATTAATGCATCTATCAGTAGCATTAACGCCTACAAACAAGCCGTAGTTTCCGCTCAAAGCTCATTAGACGCGATGGAAGCGGGCTACTCGGTCGGTACGCGTACCATTGTTGATGTGTTGGATGCGACCACCACGTTGTACAACGCCAAGCAAGAGCTGGCGAATGCGCGTTATAACTACCTGATTAATCAGCTGAATATTAAGTCAGCTCTGGGTACGTTGAACGAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGTCGTTCAGCAAACATCCGCACGCACTACCACCAGTAACGGTCATAACCCTTTCCGTAACTGATGACGACGACGGGGCTTCGGCCCCGTCTGAACG  >  NC_000913/3179144‑3179629

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