Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,539,524	A→G	intergenic (‑23/‑1433)	nanC ← / → fimB	N‑acetylneuraminic acid outer membrane channel/regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,539,524	0	A	G	100.0%	122.8 / NA	34	intergenic (‑23/‑1433)	nanC/fimB	N‑acetylneuraminic acid outer membrane channel/regulator for fimA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/15);  total (19/15)

CGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAAT  >  NC_000913/4539390‑4539657                                                                                                                                       |                                                                                                                                      cGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATATGTAGGATGCAGaa                                                                                                                                  >  2:331409/1‑140 (MQ=255)    ggCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGc                                                                                                                               >  2:287058/1‑140 (MQ=255)    ggCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGc                                                                                                                               <  1:139067/140‑1 (MQ=255)    ggCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGc                                                                                                                               <  1:257536/140‑1 (MQ=255)                   atCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCa                                                                                                                <  2:490195/140‑1 (MQ=255)                        ccacGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAAc                                                                                                           <  1:331409/140‑1 (MQ=255)                            gtacgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAAtt                                                                                                       <  2:412700/140‑1 (MQ=255)                              acgtCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAg                                                                                                     >  2:480320/1‑140 (MQ=255)                                 tCCAGTGTCGCAGCCTGAGAAATTAATGGGGACTAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAAt                                                                                                        <  1:36868/134‑1 (MQ=255)                                 tCCAGTGTCGCAGCCTGAGAAATTAATGGGGACTAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAAt                                                                                                        >  2:36868/1‑134 (MQ=255)                                      tgtCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtg                                                                                             >  2:200537/1‑140 (MQ=255)                                        tCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgtt                                                                                           >  1:108074/1‑140 (MQ=255)                                          gCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTAtgttgt                                                                                         >  2:389336/1‑140 (MQ=255)                                           cAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTtgtg                                                                                        >  2:288247/1‑140 (MQ=255)                                                  agaAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATAtt                                                                                 >  1:31609/1‑140 (MQ=255)                                                          aTGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGcc                                                                         <  1:480320/140‑1 (MQ=255)                                                                                     cGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAAcc                                                 <  2:537915/137‑1 (MQ=255)                                                                                     cGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAAcc                                                 >  1:537915/1‑137 (MQ=255)                                                                                     cGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGt                                              >  1:79556/1‑140 (MQ=255)                                                                                       ccAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCgg                                               >  2:439376/1‑137 (MQ=255)                                                                                       ccAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCgg                                               >  1:504979/1‑137 (MQ=255)                                                                                       ccAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCgg                                               <  1:439376/137‑1 (MQ=255)                                                                                       ccAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCgg                                               <  2:504979/137‑1 (MQ=255)                                                                                       ccAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAg                                            >  2:161760/1‑140 (MQ=255)                                                                                       ccAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAg                                            >  2:474364/1‑140 (MQ=255)                                                                                       ccAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAg                                            >  1:427533/1‑140 (MQ=255)                                                                                         aGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCt                                          <  1:51574/140‑1 (MQ=255)                                                                                           aaaGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGt                                        >  1:108553/1‑140 (MQ=255)                                                                                                    tAGCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTatat                               >  2:108555/1‑140 (MQ=255)                                                                                                      gCCTTTTTCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATAttt                             >  1:432195/1‑140 (MQ=255)                                                                                                        cTTTTTCATAAATTTCACTCGTTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATAttttt                           <  1:395658/140‑1 (MQ=255)                                                                                                         tttttCATAAATTTCACTCATTTGTAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATAtttttt                          <  1:474364/140‑1 (MQ=255)                                                                                                                                 tAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAAt  <  2:432195/140‑1 (MQ=255)                                                                                                                                 tAGGATGCAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAAt  <  2:292049/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      CGTGGCTTCCACTACGATATCCACCACGTACGTCCAGTGTCGCAGCCTGAGAAATTAATGGGGACGAAAAGCACAGTAATAATACGCCAGAAAGTATTTTAGCCTTTTTCATAAATTTCACTCATTTGTAGGATACAGAAAGCAATACAAAGCCCGCATAAACAATTAGCATTTATGTTGTGTAATATTTTTTTGCCAGGCTTATAGTGTCTTTGGCAACCGGTAGCTGTATTTTATATTTTTTTGTATAAGGTCTCCTGTGAAAAAT  >  NC_000913/4539390‑4539657

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