Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,292,312	C→T	P324P (CCC→CCT)	yraK →	putative fimbrial adhesin YraK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,292,312	0	C	T	100.0%	150.3 / NA	43	P324P (CCC→CCT)	yraK	putative fimbrial adhesin YraK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/16);  total (27/16)

CCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCCGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTT  >  NC_000913/3292072‑3292558                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       ccGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCtt                                                                                                                                                                                                                                              <  1:462099/251‑1 (MQ=255)       tgccatgccaGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGc                                                                                                                                                                                                                                        >  2:399593/1‑251 (MQ=255)        gccatgccaGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGc                                                                                                                                                                                                                                        >  1:176993/1‑250 (MQ=255)          cctgccaGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGCTGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTTCGGGCAtt                                                                                                                                                                                                                                     >  2:121972/3‑250 (MQ=255)               caGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCAcc                                                                                                                                                                                                                                >  1:362016/1‑251 (MQ=255)                   cGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAAGGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCTGGCATGGCACCGTCa                                                                                                                                                                                                                            <  1:437660/251‑1 (MQ=255)                                gTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACACTGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTg                                                                                                                                                                                                               >  2:223022/1‑251 (MQ=255)                                             tCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTg                                                                                                                                                                                                  >  2:153607/1‑251 (MQ=255)                                              cACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTg                                                                                                                                                                                                  <  2:423684/250‑1 (MQ=255)                                                  tATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCgg                                                                                                                                                                                             >  1:3457/1‑251 (MQ=255)                                                       cTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGgtggtg                                                                                                                                                                                        >  1:390427/1‑251 (MQ=255)                                                                 aCTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAAc                                                                                                                                                                              >  1:148620/1‑251 (MQ=255)                                                                      aaaGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGa                                                                                                                                                                          >  1:50846/1‑250 (MQ=255)                                                                             aaaaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCCAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTc                                                                                                                                                                   >  2:124072/1‑250 (MQ=255)                                                                               aaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAgg                                                                                                                                                                <  2:374500/251‑1 (MQ=255)                                                                                  cTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGc                                                                                                                                                              >  1:420570/1‑250 (MQ=255)                                                                                     gTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAg                                                                                                                                                          <  2:50846/251‑1 (MQ=255)                                                                                         cAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGGGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGccac                                                                                                                                                      <  2:404409/251‑1 (MQ=255)                                                                                         cAAATGCGCATCGACGGTACAACTGGGGGAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGGGCACGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGATCGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGGGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGccac                                                                                                                                                      >  1:279359/1‑251 (MQ=255)                                                                                            aTGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCAt                                                                                                                                                   >  1:375830/1‑251 (MQ=255)                                                                                                 cATCGACGGTACAACTGGGGTAGTAGACCGGTACAACCTGGTCGCCAGGCGAAGAAGGTCAGACAATGTGCCCGATGTCGGTATGCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACTGTAGCCTTCCTGTTGGTCCTTCCGGGCATGGCGCCGTCAATATGCGCGCCTGGCCAGTTAATCTGGGCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCg                                                                                                                                              <  1:441268/251‑1 (MQ=255)                                                                                                     gACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCAt                                                                                                                                          >  1:10575/1‑251 (MQ=255)                                                                                                      aCGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCGGGCGAAGAAGTTCAGAAAAGGTGCCCGATGTCGGTATTCGTAATGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATc                                                                                                                                         <  1:111556/251‑1 (MQ=255)                                                                                                       cGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCg                                                                                                                                        <  2:375830/251‑1 (MQ=255)                                                                                                               cTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTa                                                                                                                                 >  2:457134/1‑250 (MQ=255)                                                                                                                            ccAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGTCGTTGc                                                                                                                   <  2:363451/251‑1 (MQ=255)                                                                                                                             cAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCg                                                                                                                  >  1:50253/1‑251 (MQ=255)                                                                                                                              aGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGt                                                                                                                 <  1:400716/251‑1 (MQ=255)                                                                                                                                     ccTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGAtt                                                                                                          <  1:91346/251‑1 (MQ=255)                                                                                                                                                    aagaagTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGCAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTCAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCgg                                                                                           >  2:79076/1‑251 (MQ=255)                                                                                                                                                        agTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATTTTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTATCGCAACGCTTCGGCACt                                                                                       <  2:353911/251‑1 (MQ=255)                                                                                                                                                                                         aTTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACa                                                      <  1:315692/251‑1 (MQ=255)                                                                                                                                                                                            gAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCCTTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAgcg                                                   >  2:382250/1‑251 (MQ=255)                                                                                                                                                                                               aaTCTCGGCGGCGGCGTTGCAAATATTCCTATTCAGAACGGTATCCATCCTGTTGCTCATTCCGGGCATGGCACCCTCAATATGCGCGCCTGGCCAGATAATCTGGTCGGTGGAGAGCTGGCAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGCTTAACGCAACGCTTCGGCACTTTTACCCATGCTGCTACATCGCATACTTATACAgcgcct                                                >  2:73795/1‑249 (MQ=255)                                                                                                                                                                                               aaTCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCAt                                                >  2:296090/1‑251 (MQ=255)                                                                                                                                                                                                           ggAGTTGAAAATATACCTGTTCAGAAATGTATCCTTCCTGTTGATCCTTAAGAACATCGCACCGTCAATATGCGCCCTTGGCAAGTTAATCAGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAAAAGACGTTGCGTTTTGATTACCGCAACGCTTCTGCACTGCTACCCCTGCTGCTCCAGCGCATAATTATACAGCGCCTTTTTCTTCACGc                                    >  2:304849/1‑251 (MQ=255)                                                                                                                                                                                                               ttGCAAATATTCCTTTTCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTg                                >  1:237333/1‑251 (MQ=255)                                                                                                                                                                                                                    aaTATTCCTTTTCAGCACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGACttt                           >  1:387993/1‑251 (MQ=255)                                                                                                                                                                                                                      tatTCCTTTTCAGAACGGTATCCTTCCTGTTGCTCCTTCCGGGCATGGCAACGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCt                         >  2:325529/1‑251 (MQ=255)                                                                                                                                                                                                                              ttCAGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCa                  <  1:248936/250‑1 (MQ=255)                                                                                                                                                                                                                                 aGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCCCCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAgcg                                                                  <  2:95504/199‑1 (MQ=255)                                                                                                                                                                                                                                 aGAACGGTATCCTTCCTGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAgcg                                                                  >  1:95504/1‑199 (MQ=255)                                                                                                                                                                                                                                             ttCCTGTTGATACTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCtttttt  >  2:157961/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       CCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGAGTTGCAAATATTCCTTTTCAGAACGGTATCCTTCCCGTTGATCCTTCCGGGCATGGCACCGTCAATATGCGCGCCTGGCCAGTTAATCTGGTCGGTGGTGAGCTGGAAACAGGAAAATTTCAGGGCACAGCCACCATTACCGTCATCGTGCGGTAAAAGAAGACGTTGCGTTTTGATTAACGCAACGCTTCGGCACTGTTACCCCTGCTGCTCCAGCGCATACTTATACAGCGCATTTTTCTTCACGCCGTGAATTTCTGCGGCCAGCGCCGCCGCTTTTTT  >  NC_000913/3292072‑3292558

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: ‑