Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,540,431	T→C	intergenic (‑930/‑526)	nanC ← / → fimB	N‑acetylneuraminate outer membrane channel/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,540,431	0	T	C	100.0%	141.9 / NA	42	intergenic (‑930/‑526)	nanC/fimB	N‑acetylneuraminate outer membrane channel/Type 1 fimbriae regulatory protein FimB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/22);  total (20/22)

TAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCA  >  NC_000913/4540196‑4540668                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              tAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAAt                                                                                                                                                                                                                                >  1:88865/1‑251 (MQ=255)  aaTTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCAATTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCCGGATGATGTTTCACTTAGTTTGTTTGCCGTATTGCCCGGCGAATGGCTGTGATTGTGGAAGGTTTTAGTCGTAGTGGCCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATTATCTCCTGCGTTCCCCCATATCTCCAGGATTAAAAGGAATg                                                                                                                                                                                                                               >  2:569692/1‑251 (MQ=255)        catcatGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTaa                                                                                                                                                                                                                            <  1:393012/248‑1 (MQ=255)         atcatGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGGTGGTGTTTCTCTTTGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGCACAAAGCTATATTTACCAAAGCATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCCTATCTCCAGGATAAAAAGGAATGTAACAAt                                                                                                                                                                                                                        >  2:421786/1‑251 (MQ=255)               aaaTATGTTTCCTGGTTTGTGGCTTGTAACTGGTCTCTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATg                                                                                                                                                                                                                  <  2:465358/251‑1 (MQ=255)                     gTTTCCTGGTTTGTGTCTTGTAACTGGGCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCTCTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTAGATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGGGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAAGCTCATGGCGTa                                                                                                                                                                                                             <  1:243844/250‑1 (MQ=255)                             gTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGCTGTTGGTGTTTTCCGGATGATGTTTCACTTAGTTTGTTTGCCGTCTCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCCAAGCTATATTTAACAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACg                                                                                                                                                                                                     >  1:593535/1‑250 (MQ=255)                               ttGTGGCTTGTAACTGGTCCCTTCTGAAGTCGATCTGGAGAGGCTTGTTGGTGTTGGTGTTTTCAGGATGATGGTTCACTTTGTTTGTTTGCCGGATCGCCCGGCGAAATGCTGGGATTGTGGAAGGTTTAAGTCGTAGGGACCCAAGCTATATTTACCAACGAATTTAGATGAAAAAATCATCTCCTGCGGTCCCCCCTATCCCCCGGATAAAAAAGAATGTAACAATCTCCTGGCCTAACCTGACCAAt                                                                                                                                                                                                  >  2:347595/1‑251 (MQ=255)                                 gtgGCTTGTAACTTGTCACTTCTGAAGTCTATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATc                                                                                                                                                                                                 <  1:426147/250‑1 (MQ=255)                                                   cTTCTGAAGTCGATATGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCCACGAATGTAGGTGAAAAAATCATCTCCGGCGTTCCCCCATATCTCCAGGATAAAAAGGAATATAACAATATCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAgg                                                                                                                                                                              >  2:580232/1‑251 (MQ=255)                                                                   ggAGAGGCTTGTTGATGTTGGTGTGTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCGAGGGACCTAAGaa                                                                                                                                                                    <  2:210613/245‑1 (MQ=255)                                                                    gagaGGCTTGTTGATGTTGGTGGTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCAt                                                                                                                                                             <  1:217044/251‑1 (MQ=255)                                                                                      ggTGTTTTCATGATGATGTTTCACTTAGTATGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGaa                                                                                                                                           <  1:249258/251‑1 (MQ=255)                                                                                                   tgatgTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGCCGAATCAGCAGGAATAATCGCTCGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTgcgcagt                                                                                                                              >  1:572618/1‑248 (MQ=255)                                                                                                   tgatgTTTCACTTAGTTTGTTTGCCGTAGCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCAAGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACCAAGAAATTACTgcgc                                                                                                                                 <  2:572618/248‑1 (MQ=255)                                                                                                           cACTTAGGTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAAGTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAAt                                                                                                                       <  2:540181/250‑1 (MQ=255)                                                                                                                 gtttgtttGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCAtt                                                                                                                >  1:334958/1‑251 (MQ=255)                                                                                                                   ttgtttgCCGTCTCGCACGGCGAATGACTGTGATTGAGGAAGGTTTAAGTCGGAGTGACCAAAGATATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCAtttt                                                                                                              >  2:541128/1‑251 (MQ=255)                                                                                                                      tttgCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTtgtg                                                                                                           >  1:63130/1‑251 (MQ=255)                                                                                                                       ttgCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTgg                                                                                                          >  1:384171/1‑251 (MQ=255)                                                                                                                       ttaccggaTCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAAGCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGTGCTCCATGAAATAGCCATTTtgtg                                                                                                           <  2:428980/243‑1 (MQ=255)                                                                                                                         gCCGGATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGGTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCa                                                                                                        <  1:475387/251‑1 (MQ=255)                                                                                                                           cGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCaaa                                                                                                      <  1:131640/251‑1 (MQ=255)                                                                                                                                          aTGGCTGTGATTGAGGAATGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTACCATGATAATAGCCACTAAGAAATTACTTCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACtaata                                                                                       <  1:575140/251‑1 (MQ=255)                                                                                                                                                      gAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATAtgtg                                                                           >  1:143148/1‑251 (MQ=255)                                                                                                                                                            gggttAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCACCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACgg                                                                     <  1:126283/248‑1 (MQ=255)                                                                                                                                                                           gACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAat                                                       <  1:358312/250‑1 (MQ=255)                                                                                                                                                                             ccAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATAtt                                                    >  2:201727/1‑251 (MQ=255)                                                                                                                                                                               aaaGCTATATTTACCAACGAATGTAGATAAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTaa                                                  <  2:64919/251‑1 (MQ=255)                                                                                                                                                                                   cTATATTTACCAACGAATGTAGCTGACAAAATCATCTCCTGCGTTCCCCCATATATCCAGGATAACAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGAAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAAttt                                              >  2:563846/1‑251 (MQ=255)                                                                                                                                                                                     atatTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTg                                            <  1:347595/251‑1 (MQ=255)                                                                                                                                                                                              aaCGAATGTAGATGAAAAAATCATCTCCGGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGGGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAAtttt                                   <  1:266165/251‑1 (MQ=255)                                                                                                                                                                                               aCGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAAttttt                                  >  2:254015/1‑251 (MQ=255)                                                                                                                                                                                               aCGAATGTAGATGAAAAAATCATCGCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAAtttt                                             <  1:270721/240‑1 (MQ=255)                                                                                                                                                                                               aCGAATGAAGATGACAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGACTGTAACAATCTCATGGCGTCAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGCATATTAAAtttt                                             >  2:270721/1‑240 (MQ=255)                                                                                                                                                                                                cGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGaa                                                        <  1:475736/228‑1 (MQ=255)                                                                                                                                                                                                cGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGaa                                                        >  2:475736/1‑228 (MQ=255)                                                                                                                                                                                                     gTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTaa                                                                                                                                                                       >  1:503450/1‑112 (MQ=255)                                                                                                                                                                                                     gTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTaa                                                                                                                                                                       <  2:503450/112‑1 (MQ=255)                                                                                                                                                                                                      tAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTg                           >  1:121950/1‑251 (MQ=255)                                                                                                                                                                                                                    atcTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTAc             >  2:144944/1‑251 (MQ=255)                                                                                                                                                                                                                                cccccATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATca  <  1:396330/250‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                              TAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAGAACATATCA  >  NC_000913/4540196‑4540668

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