Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,540,114	A→T	intergenic (‑613/‑843)	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,114	0	A	T	100.0%	59.0 / NA	22	intergenic (‑613/‑843)	nanC/fimB	N‑acetylneuraminate outer membrane channel/Type 1 fimbriae regulatory protein FimB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (12/10);  total (12/10)

GTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCG  >  NC_000913/4539870‑4540328                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         gTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCGTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCATATTTGCAATTCGGGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑c                                                                                                                                                                                                                  <  2:619914/251‑1 (MQ=255)      gTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAAc                                                                                                                                                                                                             <  1:551027/251‑1 (MQ=255)                                         aCTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTgg                                                                                                                                                                          <  1:182096/251‑1 (MQ=255)                                            aTCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTg                                                                                                                                                                        >  1:635182/1‑250 (MQ=255)                                                             taCCTGTTATAACAGATCAAAAATAACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTATTCTAAGAGGAAGATGCCATGTGAAGCCAGACGACCACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAAA‑GA‑CTAAACAAGGGGAGCTTTGCACGCTATCTCAGTGAGCTTGGTGACAATCAGCGTTTACCCg                                                                                                                                                      >  2:124372/1‑251 (MQ=255)                                                                              cAAAAATCTCGCAATCCACAAAATAAAACCAGATGTGCAATTCGTGTCACAAAATTGGTCGATTTTTTACGAAGAGGAAGATGCCATGTGAAGCCATACGAACCCTTGCGGTGGTCTTCAAAAACCAAAGATCTTAGTTTAACTATTTTTTTTATAAATAATTTGTTTGGATGT‑‑TAAACAAGGGGAGCTTTGCAAGCTAACTCATTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATa                                                                                                                                     <  1:630669/251‑1 (MQ=255)                                                                              cAAAAATCACGCAATCCATACAACAAACCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCAt                                                                                                                                      <  1:66997/250‑1 (MQ=255)                                                                               aaaaaTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTTAAGCCAGACGAACTCTTGCTGTTGTCTTCAAAAACTAAAGATCTTATTTTAACTATTTGTTTTATAAATAATTTATTaa                                                                                                                                                                                                                         <  2:384843/168‑1 (MQ=255)                                                                               aaaaaTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTaa                                                                                                                                                                                                                         >  1:384843/1‑168 (MQ=255)                                                                                                tACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAatat                                                                                                                   >  2:501490/1‑251 (MQ=255)                                                                                                  caacaaAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGt                                                                                                                 <  1:501490/251‑1 (MQ=255)                                                                                                                   cAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCtt                                                                                                >  1:170373/1‑251 (MQ=255)                                                                                                                      ttCGTGTCACAACATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGGCTTCAAAAACTAAAGATATTAGTTTAACTCTTTGTTTTATAAATAATTTGGTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCATTGTTTACCCGCCATCAGGCTGACCATAATTCTCATCATGAAATATGTTTCATGGttttgtgcttgt                                                                                              >  2:22011/1‑241 (MQ=255)                                                                                                                        cGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAAc                                                                                           >  1:315130/1‑251 (MQ=255)                                                                                                                           gtCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCCAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTgg                                                                                        >  1:640229/1‑251 (MQ=255)                                                                                                                                              tttttCTAAGAGGAAGATGCCATGTGGAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTCTACATAATTTGTTTAGA‑GT‑CTCAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAACTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATCATTCTCATCATGCAATATGTTTCCTGTGTTGTGGATTTGAACTGGGCACttactaattcaggct                                                                     >  2:477422/1‑238 (MQ=255)                                                                                                                                               ccgccgctcttccgatccccATGTGAAGCCAGACGAACACTTGCGGTGTTTTTCAAAAACTAAAGATCTTAGTTTAATTATTTTTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTCCCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTg                                                                    <  2:138046/233‑1 (MQ=255)                                                                                                                                                       gagGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCAt                                                                                                                                      >  1:608618/1‑177 (MQ=255)                                                                                                                                                            aGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTCAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑ATAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTATACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTTTGTCTTGTTACTGGTCACTTCTGAAGTCGACATGGAGAGGCttgttt                                                       >  2:267766/1‑250 (MQ=255)                                                                                                                                                                 ccATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTg                                                                    >  1:138046/1‑233 (MQ=255)                                                                                                                                                                                                           aaaaGATCTTAGTTTAATTATTTGTTTTATAAATAATTTGGTTAGA‑GT‑CAAAACAAGGGGAGCTTTTCAAGTTAACTAAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTAt        <  2:190291/250‑1 (MQ=255)                                                                                                                                                                                                            aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGTGCTTGGTGAAAATCAGAGTTTACCCGACATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGAAACTGGACACTTCTGAATTAGATCTGGAGAGGATTGTTGATGTTGGTGGTTTCAGGATGATGGTTTAATTAatttgtttacggtatt       >  2:658018/1‑243 (MQ=255)                                                                                                                                                                                                            aaaGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAG‑‑GTACTAAACAAGGGGAGCTTGGCAAGCTAACGCAGTGAGCTTGTTGAAAATCAGTGTTTCCCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTTGCTTGTAACTGGTCACTGCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGAGGATGTTTCACTTAGTTTGTTTGCCGTATc       <  2:46542/251‑1 (MQ=255)                                                                                                                                                                                                                    tAGTTTAACTATTTGTTTTATAAATTATTTTTTTAGA‑TT‑CTAAACAAGGGGAGCTTTGCAGGTTAACTAAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAAGTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTATGCCGTATCGCCCg  <  1:547063/248‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         GTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGA‑GT‑CTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCG  >  NC_000913/4539870‑4540328

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