Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,541,261	C→T	A102V (GCT→GTT)	fimB →	Type 1 fimbriae regulatory protein FimB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,541,261	0	C	T	100.0%	138.1 / NA	40	A102V (GCT→GTT)	fimB	Type 1 fimbriae regulatory protein FimB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/22);  total (18/22)

TCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAG  >  NC_000913/4541018‑4541500                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                tCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGc                                                                                                                                                                                                                                          >  1:412247/1‑251 (MQ=255)      aaaGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTg                                                                                                                                                                                                                                     <  2:316815/251‑1 (MQ=255)      aaaGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTg                                                                                                                                                                                                                                     <  1:453491/251‑1 (MQ=255)                        ccTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGGTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTaa                                                                                                                                                                                                                   <  2:112054/251‑1 (MQ=255)                               cagcaCGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGCGTGGGTATTTTTATCACGTAAGGGGAAt                                                                                                                                                                                                            >  2:325601/1‑251 (MQ=255)                                                       tgctttgTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTa                                                                                                                                                                                     <  2:108775/250‑1 (MQ=255)                                                                    tCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTAATGAATAAAGAAGTTAAGGCTTTAAAAAACTGGATGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTAtctc                                                                                                                                                                         >  2:455755/1‑251 (MQ=255)                                                                     cATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGAACTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGAc                                                                                                                                                                      <  1:97047/251‑1 (MQ=255)                                                                            tCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTgg                                                                                                                                                                                                                                    >  1:277375/1‑182 (MQ=255)                                                                            gccGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTgg                                                                                                                                                                                                                                    <  2:277375/181‑1 (MQ=255)                                                                             ccGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGAGATTGATCTTAAGTCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGTGCGAGTGGGTATTTTTATCACGTAAGGGGAATTCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCggtg                                                                                                                                                              <  1:244289/251‑1 (MQ=255)                                                                               gggCGAGTGAAATTTGTCGATTGAGGATTTCGGCTATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCggtggt                                                                                                                                                            >  2:410178/1‑251 (MQ=255)                                                                               gggCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCggtggt                                                                                                                                                            >  1:510111/1‑251 (MQ=255)                                                                                       gAAATTTGTCGTTTGAGGAGTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATTACCATCGATTAAAAAAAGGCTTTTCAACATCGCTCCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCGTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGcctg                                                                                                                                                    <  2:438179/251‑3 (MQ=255)                                                                                                 gATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTg                                                                                                                                          <  2:297855/251‑1 (MQ=255)                                                                                                         aTTTCGGGTAATGATCTTCAGGCCAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAAAGCACACGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGACGGGTTGTCACTGGAGATTc                                                                                                                                   >  2:249665/1‑250 (MQ=255)                                                                                                             cGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGGATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTGTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATcc                                                                                                                               <  2:530356/250‑1 (MQ=255)                                                                                                                      aTCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCTGCTATTGACTAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGCGTGGGTATTTTTCTCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGtt                                                                                                                      >  2:575982/1‑250 (MQ=255)                                                                                                                                     gtgtATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGtgtg                                                                                                        <  1:7435/251‑1 (MQ=255)                                                                                                                                           atatatCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGtttt                                                                                                   <  2:209995/250‑1 (MQ=255)                                                                                                                                               atCCATCGATTAAAAAAAGGCTTTTCGACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCttt                                                                                              >  1:626/1‑251 (MQ=255)                                                                                                                                                tCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTg                                                                                             >  2:446263/1‑251 (MQ=255)                                                                                                                                                  cATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGc                                                                                           >  2:25240/1‑251 (MQ=255)                                                                                                                                                       aTTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAata                                                                                      <  1:573108/251‑1 (MQ=255)                                                                                                                                                                    ttttCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATAc                                                                         >  1:314854/1‑251 (MQ=255)                                                                                                                                                                          acaacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACt                                                                   >  1:504098/1‑251 (MQ=255)                                                                                                                                                                            aacGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATTCGCGACTTa                                                                 <  2:173655/251‑1 (MQ=255)                                                                                                                                                                              cGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTGTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTTCCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGGGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATc                                                               <  2:48894/251‑1 (MQ=255)                                                                                                                                                                               gCACCCGCAATTGAATAAAGAAGTTCAGGCTTTAAACAACTGGTTGAGAATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCTCTGGAGATTCATCCGCACATGTTACGCCGTTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATcc                                                              >  2:231011/1‑251 (MQ=255)                                                                                                                                                                                       tATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGcc                                                                                                                                                      <  2:267796/155‑1 (MQ=255)                                                                                                                                                                                       tATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGcc                                                                                                                                                      >  1:267796/1‑155 (MQ=255)                                                                                                                                                                                              tAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGc                                               >  1:90176/1‑251 (MQ=255)                                                                                                                                                                                               aaaGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCa                                              <  1:25240/251‑1 (MQ=255)                                                                                                                                                                                                       tCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATTCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATGCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATGATCTTGGGCATCGCAat                                       <  2:40327/250‑1 (MQ=255)                                                                                                                                                                                                         aGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATAtt                                    <  1:11735/251‑1 (MQ=255)                                                                                                                                                                                                                aaaacaCTGGTTGAGTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGTATTATCTTGGACATCGCAATATTCGTCAt                              >  2:59009/1‑250 (MQ=255)                                                                                                                                                                                                                              gTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTGTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATAcc               <  2:38140/251‑1 (MQ=255)                                                                                                                                                                                                                              gTATCCGTACTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTGACCATATAATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCGCATGTTACGCCATTCGTGTGGTTTTGCTGTGGCGAATATGGGAATAGATACGCGACGTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATAcc               <  1:575982/251‑1 (MQ=255)                                                                                                                                                                                                                                       cTTCGTACCCGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTAGACTTCCGGTGGTCATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAAt      >  1:578114/1‑251 (MQ=255)                                                                                                                                                                                                                                           gtccccGCATGTTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAg  <  1:446263/248‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                TCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAA‑‑GTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTACTTCGTACCCGCATGCTGAGAGCGAGTGGGTATTTTTATCACGTAAGGGGAATCCGCTTTCTCGGCAACAGTTTTACCATATTATCTCGACTTCCGGTGGTAATGCCGGGTTGTCACTGGAGATTCATCCGCACATGTTACGCCATTCGTGTGGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAG  >  NC_000913/4541018‑4541500

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