Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	305,439	G→A	S460F (TCC→TTC)	ecpD ←	fimbrial adhesin EcpD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	305,439	0	G	A	100.0%	132.5 / NA	39	S460F (TCC→TTC)	ecpD	fimbrial adhesin EcpD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/19);  total (20/19)

CGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGG  >  NC_000913/305195‑305683                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      cGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGGACGCAGCGAAATGGCGTTATCCAGTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGTTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAgt                                                                                                                                                                                                                                                  <  1:92178/250‑1 (MQ=255)       aaCATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCg                                                                                                                                                                                                                                           <  1:396298/251‑1 (MQ=255)        aCATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCgg                                                                                                                                                                                                                                          >  1:378761/1‑251 (MQ=255)             aaTCTCTTCTGAAGCGCTGACTTCGCCAAACCAGCAGTTATCATCTACCGTACGCAGCGCAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGCCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCTCACCCGGCCTCGTAGGTTTTTGTAGCTCCGTAGCGGGTAATAACGAAACGTGTAGCGGGGAACg                                                                                                                                                                                                                                     >  2:49768/1‑251 (MQ=255)             aaTCTCTCCTGAAGCGCTTACTTCGCCAAACCAGCCGTTCTCATCTTCCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTTGTCTTATCCATCTGCCCCACTTCGCCAGAGCTATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCTCCAGCTATCATCGCACCCGGCATCGTAGGTTTTGGTAGCTCCGTTGAGGGTAATAAAGAAAAGTGTCGCCGGTAACg                                                                                                                                                                                                                                     <  1:247930/251‑1 (MQ=255)              aTCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATCAAGAAAAGTGTCGCCGGGAACgg                                                                                                                                                                                                                                    >  1:220865/1‑251 (MQ=255)                    ccTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑Ac                                                                                                                                                                                                                                <  1:166759/248‑1 (MQ=255)                                        aaaCCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGc                                                                                                                                                                                                         >  1:263198/1‑251 (MQ=255)                                                           aCCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAgg                                                                                                                                                                                          <  2:354993/247‑1 (MQ=255)                                                              gTACGCAGCGAAATGGCGTTATCCATTGGAATCGCAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATc                                                                                                                                                                                                              >  1:378560/1‑224 (MQ=255)                                                              gTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGTGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCGTACCGTCATc                                                                                                                                                                                                              <  2:378560/224‑1 (MQ=255)                                                                                 tATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTAGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGt                                                                                                                                                                >  2:205200/1‑251 (MQ=255)                                                                                       ttGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTg                                                                                                                                                                       <  1:10440/238‑1 (MQ=255)                                                                                       ttGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTg                                                                                                                                                                       >  2:10440/1‑238 (MQ=255)                                                                                              cGAAAATTTGACTGTGGTCTTCTCCACCGACCCCACTTCGCAAGCGATATCAGTCAACGGTGTGGTCAACCCCAGCGAACCGGTCATATCCCGCCAGTCATCATCGCACCCGGCATCGTATGTTTTTGTAGCCCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑CCTTTCGCCTTACCTTCTTCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTg                                                                                                                                                    <  2:106205/250‑1 (MQ=255)                                                                                                       gACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACt                                                                                                                                           <  1:205200/250‑1 (MQ=255)                                                                                                              gTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCAACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGc                                                                                                                                   >  1:152446/1‑251 (MQ=255)                                                                                                                                 cTTCGCCAGAGATCTCAGTCCACGGTGTGGTCAACCACCGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCtggtg                                                                                                                >  2:269664/1‑251 (MQ=255)                                                                                                                                               tCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGcc                                                                                                  >  1:267238/1‑251 (MQ=255)                                                                                                                                                     cACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATc                                                                                            >  1:237709/1‑251 (MQ=255)                                                                                                                                                           gtgGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ATTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCTCCTTGATCATCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGc                                                                                      <  1:269664/251‑1 (MQ=255)                                                                                                                                                             ggTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCg                                                                                    >  1:201973/1‑251 (MQ=255)                                                                                                                                                                cAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACt                                                                                 >  2:313782/1‑251 (MQ=255)                                                                                                                                                                         cgcATCGGTCATATCCCGCCAGCTATCATCGCCCCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGAAATAAAGAAAAGTGTCGCCGGGAAAGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGc                                                                        >  1:18530/1‑251 (MQ=255)                                                                                                                                                                                 tCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACatat                                                                >  1:72924/1‑251 (MQ=255)                                                                                                                                                                                       cccGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATAtccctcc                                                          <  1:243338/251‑1 (MQ=255)                                                                                                                                                                                       cccGCCAGCTATCAGCGCCCCCGGCATCGTAGGTTTTTGTAGCTCCGTTTCGGGTAATAAATAAAAGGGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCCTCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATAtccctc                                                           <  1:354151/250‑1 (MQ=255)                                                                                                                                                                                            cAGCTATCATCGAACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAACGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGCGCTGAAGACACAATAGGAGCGCCCGCCAATCTCCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATTTAGCCGAACTCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCcgcg                                                      >  2:341005/1‑250 (MQ=255)                                                                                                                                                                                              gctttctccgcccccGGCTTCGTAGGTTTTTGTAGCCCTGTGGCGGGTAAAAAAGAAAAGTGTCGCCGGTAATGGT‑CCTTGCGCCGTACCTGCTCCGGAGCTAAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTCCCGCTGGTGGTTACGCTATAGCCGAAATCGAGTGCGGACTCGCCGCTGCCAACATACCCCTCCCTCGACg                                                   <  2:306072/239‑1 (MQ=255)                                                                                                                                                                                                 atcatcGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑CCTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACgg                                                  <  2:267238/249‑1 (MQ=255)                                                                                                                                                                                                         aCCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATCAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGt                                         >  1:38850/1‑250 (MQ=255)                                                                                                                                                                                                         aCCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGt                                         >  2:126860/1‑250 (MQ=255)                                                                                                                                                                                                             ggCATCGTAGTTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATAtt                                    <  1:319211/251‑1 (MQ=255)                                                                                                                                                                                                              gCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTc                                   >  1:180663/1‑251 (MQ=255)                                                                                                                                                                                                                aTCGTAGGTTTTTGTATCCCCGTTGCGGGTAATAAAGAAAAGGGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAAGAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCAt                                 <  1:323426/251‑1 (MQ=255)                                                                                                                                                                                                                      gtttttttTAGCTCCGTTGCTGTTAATAAAGAAAAGTTTCGCCGGGAACGGT‑CCTTTCCCCTTCCCGTCATCGGAGCTGAAGCCACAATGGGAGCCCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGgaga                           <  2:263198/249‑1 (MQ=255)                                                                                                                                                                                                                               aGCTCCGTTGCGGGTAATAAAGAAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTg                   >  1:348561/1‑250 (MQ=255)                                                                                                                                                                                                                                  cccGTTGCGGTAAATAAAGAAAGGTGTCGCCGGGAACGGT‑CCTTCCGCCTTACCGTCATCGGAGCTGAAGACCCAAAAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGc               <  2:150481/250‑1 (MQ=255)                                                                                                                                                                                                                                               ataaAGAAAGGTTTCGCCGG‑AACGGTTACTTTCGCCTTCCCGTCACCGTGGCGGAAGACACAAAAGAAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACgggg  <  2:328957/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CGCCTGAACATGAATCTCTCCTGAAGCGCTGACTTCGCCAAACCAGCCGTTATCATCTACCGTACGCAGCGAAATGGCGTTATCCATTGGAATCGAAAATTTGACTGTGGTCTTATCCATCTGCCCCACTTCGCCAGAGATATCAGTCCACGGTGTGGTCAACCACAGCGCATCGGTCATATCCCGCCAGCTATCATCGCACCCGGCATCGTAGGTTTTTGTAGCTCCGTTGCGGGTAATAAAGGAAAGTGTCGCCGGGAACGGT‑ACTTTCGCCTTACCGTCATCGGAGCTGAAGACACAATAGGAGCGCCCGCCAATCACCTGCGCGGGTCCGGTCACCTTGATCAGCACTTCGTCGGCAGCTGTTTTACCGCTGGTGGTTACGATATAGCCGAAATCGAGTGCCGACTCGCCGCTGCCAACATATCCCTCCCGCGACGGTGTCTGAGTATATTCATCGGAGATAATACTGATGCTGAAATCACGGGG  >  NC_000913/305195‑305683

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