Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,450,815	A→G	G87G (GGT→GGC)	yfcQ ←	putative fimbrial protein YfcQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,450,815	0	A	G	100.0%	179.8 / NA	51	G87G (GGT→GGC)	yfcQ	putative fimbrial protein YfcQ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/26);  total (25/26)

GCTTTTGCGTAAAATCAGTGGCGTTGATAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACTTCAAT  >  NC_000913/2450682‑2450949                                                                                                                                      |                                                                                                                                       gCTTTTGCGTAAAATCAGTGGCGTTGATAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTc                                                                                                                                  >  1:871016/1‑140 (MQ=255)       gCGTAAAATCAGTGGCGTTGATAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTACAGGTg                                                                                                                            <  2:1174172/140‑1 (MQ=255)                  gTGGCGTTGATAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCAtctc                                                                                                                 <  1:22601/140‑1 (MQ=255)                            tAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGa                                                                                                       <  2:294275/140‑1 (MQ=255)                                  gCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACgg                                                                                                 >  1:777445/1‑140 (MQ=255)                                   cGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACggg                                                                                                <  2:871016/140‑1 (MQ=255)                                        tATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGgtcgtc                                                                                           <  2:326892/140‑1 (MQ=255)                                         aTTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTCTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGtcgtcg                                                                                          >  2:151015/1‑140 (MQ=255)                                              gTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTca                                                                                     <  2:508263/140‑1 (MQ=255)                                                ttAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTcaca                                                                                   >  1:791711/1‑140 (MQ=255)                                                ttAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTcaca                                                                                   <  2:545260/140‑1 (MQ=255)                                                  aaaaCGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAg                                                                                 >  1:59616/1‑140 (MQ=255)                                                  aaaaCGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAg                                                                                 >  1:782159/1‑140 (MQ=255)                                                    aaCGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTcaca                                                                                   <  1:503568/136‑1 (MQ=255)                                                    aaCGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTcaca                                                                                   >  2:503568/1‑136 (MQ=255)                                                       gCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGt                                                                            <  2:286210/140‑1 (MQ=255)                                                          gACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCa                                                                         <  1:958108/140‑1 (MQ=255)                                                           aCCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGt                                                                            >  1:947038/1‑136 (MQ=255)                                                           aCCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGt                                                                            <  2:947038/136‑1 (MQ=255)                                                           aCCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCat                                                                        >  2:382046/1‑140 (MQ=255)                                                                 gTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCa                                                                  <  2:1040037/140‑1 (MQ=255)                                                                  tGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCAc                                                                 <  1:1126386/140‑1 (MQ=255)                                                                       aaGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCtt                                                            <  2:420481/140‑1 (MQ=255)                                                                         gCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTa                                                          >  2:764686/1‑140 (MQ=255)                                                                         gCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTa                                                          >  1:419192/1‑140 (MQ=255)                                                                           tCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACg                                                        <  2:59616/140‑1 (MQ=255)                                                                            cGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGa                                                       >  2:886794/1‑140 (MQ=255)                                                                             gATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGAc                                                      >  1:363838/1‑140 (MQ=255)                                                                              aTACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGAcc                                                     <  1:946323/140‑1 (MQ=255)                                                                                aCCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCAt                                                   <  2:419192/140‑1 (MQ=255)                                                                                         gggTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTCAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGt                                          >  2:215352/1‑140 (MQ=255)                                                                                          ggTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGtt                                         >  1:1096789/1‑140 (MQ=255)                                                                                             aCATCAGTCTCAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGAt                                      <  1:151015/140‑1 (MQ=255)                                                                                               aTCAGTTTCAATTGCTGCATCATTGAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATgt                                    <  2:413804/140‑1 (MQ=255)                                                                                                       cAATTGCTGCATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATa                            <  2:1096789/140‑1 (MQ=255)                                                                                                            gctgcATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATgag                       >  2:520923/1‑140 (MQ=255)                                                                                                             ctgcATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATgtcgt                                 >  2:1173051/1‑129 (MQ=255)                                                                                                             ctgcATCATTAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATgtcgt                                 <  1:1173051/129‑1 (MQ=255)                                                                                                                 atcatTAAAACTGGTTTGTGTGCCAGTCCAGGCGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATcgc                  >  1:521682/1‑140 (MQ=255)                                                                                                                     ttAAAACTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGaaa              >  1:750603/1‑140 (MQ=255)                                                                                                                        aaaCTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTc           >  1:875758/1‑140 (MQ=255)                                                                                                                          aCTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGa         >  1:384075/1‑140 (MQ=255)                                                                                                                          aCTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGa         >  2:180602/1‑140 (MQ=255)                                                                                                                          aCTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGa         <  2:786103/140‑1 (MQ=255)                                                                                                                          aCTGGTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGa         >  1:786103/1‑140 (MQ=255)                                                                                                                             ggTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACtt      >  2:1174924/1‑140 (MQ=255)                                                                                                                              gTTTGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACTTc     >  2:412423/1‑140 (MQ=255)                                                                                                                               tttGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACTTCa    <  1:382046/140‑1 (MQ=255)                                                                                                                               tttGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACTTCa    <  2:807137/140‑1 (MQ=255)                                                                                                                               tttGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACTTCa    <  2:1136336/140‑1 (MQ=255)                                                                                                                                ttGTGTGCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACTTCaa   <  1:447795/140‑1 (MQ=255)                                                                                                                                 tgtgtgCCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACTTCAAt  >  1:255034/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       GCTTTTGCGTAAAATCAGTGGCGTTGATAGCAAGCGGAGTATTGAGTTTAAAACGCTGACCATCGTGCTGAAGCTCGATACCAAAACCGGGTACATCAGTTTCAATTGCTGCATCATTAAAACTGGTTTGTGTACCAGTCCAGGTGAGGGTCATCTCCCAATCTGGATGACGGGTCGTCGAGTCACAAGTGAGTTCATATGGCACCTCTTTACGACCATAGTTGCCGTTGATGTCGTCGATAATGAGATCGCGAAATTCGACTTCAAT  >  NC_000913/2450682‑2450949

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: ‑