Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,001,463	G→A	A750T (GCG→ACG)	elfC →	putative fimbrial usher protein ElfC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,001,463	0	G	A	95.8%	70.0 / NA	24	A750T (GCG→ACG)	elfC	putative fimbrial usher protein ElfC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (8/15);  minor base C (1/0);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATAGCGACGCTGGACGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGA  >  NC_000913/1001217‑1001702                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                ggCTACAACTACAGCGATAACGGCAGTCATCAACAACTGAACTACGCTCTTTCAGACAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGAcgtc                                                                                                                                                                                                                                             <  1:328226/251‑1 (MQ=255)     aCAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGACAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGtt                                                                                                                                                                                                                                          <  1:48918/250‑1 (MQ=255)                         aGTCAACAACAACTGAACTACGCTCTTTCAGACAGTTTAGTTGCCCATTCACAGGGCATTAACCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCAACTTCTTATCGGGAAAATCGAATCGCACATGATGCGCCGTCGTTAAAACGTAACGTGGATCTTg                                                                                                                                                                                                                      >  2:15253/1‑250 (MQ=255)                          gTCAACAACAACTGAACTACGCTCTTTCAGACAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGaa                                                                                                                                                                                                                    <  1:129511/251‑1 (MQ=255)                                aacaacTGAACTACGCTCTTTCAGACAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAACCTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAGATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGAGCTTGAAAATGc                                                                                                                                                                                                               >  1:116334/1‑250 (MQ=255)                                                        acAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCaa                                                                                                                                                                                       >  1:413486/2‑250 (MQ=255)                                                                             aGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAg                                                                                                                                                                 <  1:352569/251‑1 (MQ=255)                                                                                                      gggggAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAgg                                                                                                                                         <  2:5218/250‑1 (MQ=255)                                                                                                        gggAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACCGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGCTGCGACGTCGTTAAAACGTAACGTGGATCTCGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGt                                                                                                                                       >  2:138164/1‑250 (MQ=255)                                                                                                         ggAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGGAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTAt                                                                                                                                     <  2:136612/251‑1 (MQ=255)                                                                                                           aaaCTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTAtt                                                                                                                                    >  1:420310/1‑250 (MQ=255)                                                                                                             aCTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGAGCAAGGGTATTAAt                                                                                                                                 >  1:449863/1‑251 (MQ=255)                                                                                                                  cGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGaaa                                                                                                                             >  2:275099/1‑250 (MQ=255)                                                                                                                  cGTNCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCTCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGATTTCAATGCCCATGCGGGTGCAAGGGTATTAATGaaaa                                                                                                                            <  2:405743/251‑1 (MQ=255)                                                                                                                           tGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGGGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAg                                                                                                                    <  1:55066/250‑1 (MQ=255)                                                                                                                                              aaTACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGTCTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGttt                                                                                                >  2:386085/1‑251 (MQ=255)                                                                                                                                                                                  gactgGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCTCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATAGCGACGCTGGACGGCGTTCAGGCTa                                                             <  1:389962/250‑1 (MQ=255)                                                                                                                                                                                           ggANATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATAGCGACGCTGGACGGCGTACAGGCTAATAGCGGCAt                                                   <  2:420310/251‑1 (MQ=255)                                                                                                                                                                                               ataCCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATAGCGACGCTGGACGGCGTACAGGCTAATAGCGGCATAAtt                                               <  2:29371/251‑1 (MQ=255)                                                                                                                                                                                                       gttCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATAGCGACGCTGGACGGCGTACAGGCTAATAGCGGCATAATtgatgatg                                        <  2:332075/250‑1 (MQ=255)                                                                                                                                                                                                                gCCNCTTCTTATCGGGAAAATCGAATCGCACTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATATCGACGCTGGACGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGgctcgct                              <  2:409128/251‑1 (MQ=255)                                                                                                                                                                                                                            cGGGAAAATCGAATCTCACTTGATTCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTTGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGGAAGGGTATTAATGAAAGCATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCTATAGCGACGCTGGACGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCgg                   <  2:174399/250‑1 (MQ=255)                                                                                                                                                                                                                                  aatcgaatcgCACTTGATGCGACGTAGTTAAAACGTAACGTGGATCTTGAAAAGGCAGTAGTCAACGTGGTTCCCACCGAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATAGCGACGCTGGACGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGGGGGtttgcc            <  2:449863/251‑4 (MQ=255)                                                                                                                                                                                                                                             aCTTGATGCGACGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATAGCGACGCTGGACGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTAACGGGGAAGGGa  >  2:341456/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                GGCTACAACTACAGCGATAACGGCAGTCAACAACAACTGAACTACGCTCTTTCAGGCAGTTTAGTTGCCCATTCACAGGGCATTACCCTGGGGCAATCGCTGGGGGAAACTAACGTTCTGATTGCAGCACCAGGCGCAGAGAATACTCGTGTGGCGAACAGCACCGGGCTGAAAACTGACTGGCGCGGATATACCGTTGTTCCTTATGCCACTTCTTATCGGGAAAATCGAATCGCACTTGATGCGGCGTCGTTAAAACGTAACGTGGATCTTGAAAATGCAGTAGTCAACGTGGTTCCCACCAAAGGGGCGTTGGTTCTGGCGGAGTTCAATGCCCATGCGGGTGCAAGGGTATTAATGAAAACATCAAAGCAGGGTATACCGCTGCGTTTTGGCGCGATAGCGACGCTGGACGGCGTACAGGCTAATAGCGGCATAATTGATGATGATGGCTCGCTCTATATGGCGGGTTTACCGGCGAAGGGA  >  NC_000913/1001217‑1001702

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