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%	78.4 / NA	27	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 (10/17);  total (10/17)

CTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAG  >  NC_000913/4540189‑4540659                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     cTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGGTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATaa                                                                                                                                                                                                                               >  1:26868/1‑250 (MQ=255)                atcatGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCCAAGCTATATTTACCAACGAATTTAGATGAAAAAATCATTTCCTGCGTTCCCCCATATCTCCAGGATAAAAAAGAATGTaaacat                                                                                                                                                                                                               >  2:452782/1‑247 (MQ=255)                    tGAAATATGTTTCCTGGTTTGTGGCTTGTAACTTGGCACTTCTGAAGTCGATCTTGAGGGGCTTGGTGATGTTTGTTTTTTCAGGATGATGTTTTACTTAGTTTGTTTGCCGGATATCCCGGGGGATGGAGGTGGTTGGGGGAAGTTTAATGCGTATTGACCAAAGCTATATTTACCAAAGAATGTAGATTGAAAAATTAGTTCCTTCGTTCCCCCATATCTCCCTGATAAAAAAGAAtatataaatcct                                                                                                                                                                                                            >  2:145070/1‑239 (MQ=255)                     gAAATATGTTTCCTGGTTTGTGGCTTGGAACTGGTCACTTCTGAAGTCGATCTGGAGGGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGGATCGCCCGGCGCATGGATGTGATTTCGGAAGGTTTCAGGCGTAGCGACCACAGCTTTATTTACCAACGCATGAAGATGAAAAAAACCTCTCCTGCGATCCCCCATATCTCCAGGGTAAAAAGGAATATAACAAtctcc                                                                                                                                                                                                           >  1:328029/1‑249 (MQ=255)                                   ggTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGa                                                                                                                                                                                              <  1:103104/249‑1 (MQ=255)                                                                      aTCTGGAGAGGCTTGTTGCGGTTGGTGTTTTCAGGCTGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGGGATTGAGGAAGGTTTAAGTCGGAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCTGGATAAAAAGGAATGTAACAATCTCCTGGCGTAAGCTGACGAATCAGCAGGAAAAATCGCTAGGGACCTaaaaaat                                                                                                                                                         >  2:55458/1‑246 (MQ=255)                                                                        cTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATTTTTCACTTAGTTTGTTTGCCGGATCGCCCGGCGAATGGCTGTGGTTGAGGAAGGTTTAAGTCGTAGTGACCCAAGCTATATTTACCAACGAATGTAGGTGAAAAAATCATCTCCTGCGGTCCCCCATCTCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGCAGAATCAGCAGGGATAATCGCTAGGGACATAAGCAtt                                                                                                                                                         >  1:393594/1‑249 (MQ=255)                                                                          ggAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGc                                                                                                                                                      <  1:453067/250‑1 (MQ=255)                                                                                       gattttGGTGTTTTCAGGAAGATGTTTCACGTAGTTTGTTGGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCAc                                                                                                                                        <  2:428951/247‑1 (MQ=255)                                                                                                 ttttCCGGATGGTGTTTCTCTTAGTTTGTTTGCCGTATCGCCCGGCGAATGACTGTGATTGAGGAAGGTTTAAGTCGTAGGGACCAAAGATATATTTACCAAAGAATGAAGATGAAAAAATCATCTCCTGCGTTCCCCCAAATCTCCAGGATAAAAAGGAAGGTAACAATCTCATGGCGTATGCTGACGCATCAGCAGGAATAATAGCAAAGGACCTAATAATTAGCATGATAATAGCCACTAAGAAAtt                                                                                                                               >  2:441773/1‑250 (MQ=255)                                                                                                               tttCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGa                                                                                                                 >  1:59629/1‑250 (MQ=255)                                                                                                               tttCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGa                                                                                                                 <  2:59629/250‑1 (MQ=255)                                                                                                                  cACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATa                                                                                                             <  1:420574/251‑1 (MQ=255)                                                                                                                        gtttgtttGCCGTAACGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCAtt                                                                                                       <  2:310174/251‑1 (MQ=255)                                                                                                                          ttggttgCCGTATCGCCCGGCTAATGGCTGTGATTGAGTAAGGTGTAAGTCGTAGTTCCCACAACTATATTTACCAACGAATGTAGATGAAAAAAACATCTTCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCGCATGGGTTAAGCTGACGAGTCAGCAGGAATATGCGCTAGGTACCTAAGGATGAGCAGGACAATAGCCACGATTAAATTACGGTGCTCCATGGAATATCCAttt                                                                                                      <  2:301015/246‑1 (MQ=255)                                                                                                                                       cGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAg                                                                                        <  1:452782/251‑1 (MQ=255)                                                                                                                                          ccGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAAGAGCCATTTTGTGGCAAATGGAGTTg                                                                                     <  1:337242/251‑1 (MQ=255)                                                                                                                                                aaTGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACtaat                                                                               <  1:477256/251‑1 (MQ=255)                                                                                                                                                         gATTGAGGAAGGTTTAAGTCGTAGGGACCAAAGCTATATTGACCAACGAATGTAGAGGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGGCCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCat                                                                       <  2:75068/250‑1 (MQ=255)                                                                                                                                                                               aGTGACCCAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGATACACCATATCTCCAGGATAAAAAGGAATGAAACAATCTCATGGCGTAAGATGACGAATCAGAAGGAATAATCGCTAGGGACATAAGAATTCGCATGATAATAGCCACTAAGAAATTACAGCGCTCCATGCAATAGCCATTTTGTGACAAATGGAGTTAACTAATACTGTCATATGTGCGACGGATAGTTGAACg                                                 >  1:327322/1‑250 (MQ=255)                                                                                                                                                                               aGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACg                                                 <  1:55458/250‑1 (MQ=255)                                                                                                                                                                               aGTGACCAAAGCTATAATTAACAACGAATGTAGATTAAAAAATCATCTCCTGCGTTCCACCATATCTCCAGGATAAAAAGGAATGTAACAACCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGACTTAGCATGATAATAGCCTCTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGTCAAATGGAGTTGACTAATAATTTCATATGTGAGATGGCTAGTTGAACGa                                                <  2:451037/251‑1 (MQ=255)                                                                                                                                                                                            atatTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTg                                   >  2:206007/1‑251 (MQ=255)                                                                                                                                                                                              attaaCCAACGAATGTAGATGAAAAATGCATCCCCTGCGTTCCCCCATACCTCCAGTAGAAAAAGGAATGTAACAATCTCATGGCGAAAGCTGACGAATCAGCAGGAATATTCGCTAGGGCCTTAAGAATTAGCATTATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCAGTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAAGTTTGc                                  <  2:328029/246‑1 (MQ=255)                                                                                                                                                                                                   ccAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGGAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAAtt                            <  2:341697/251‑1 (MQ=255)                                                                                                                                                                                                                       aagcatcTCCTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTg        <  2:191372/248‑1 (MQ=255)                                                                                                                                                                                                                               ccTGCGTTCCCCCATATCTCCAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCGGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTGGCTGAATTGTTTATGTTGATTTTACTTGTTACAg  <  1:145070/249‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     CTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGAAGGTTTAAGTCGTAGTGACCAAAGCTATATTTACCAACGAATGTAGATGAAAAAATCATCTCCTGCGTTCCCCCATATCTCTAGGATAAAAAGGAATGTAACAATCTCATGGCGTAAGCTGACGAATCAGCAGGAATAATCGCTAGGGACCTAAGAATTAGCATGATAATAGCCACTAAGAAATTACTGCGCTCCATGAAATAGCCATTTTGTGGCAAATGGAGTTGACTAATAATGTCATATGTGAGACGGCTAGTTGAACGAATATTAAATTTTGCTGAATTTTTTATGTTGATTTTACTTGTTACAG  >  NC_000913/4540189‑4540659

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