Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,454,272	A→G	pseudogene (730/2646 nt)	yfcU ←	pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,454,272	0	A	G	100.0%	157.5 / NA	43	pseudogene (730/2646 nt)	yfcU	pseudogene, FimD fimbrial export usher family;putative membrane; Not classified; putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/19);  total (24/19)

AGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGC  >  NC_000913/2454125‑2454415                                                                                                                                                    |                                                                                                                                                aGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCg                                                                                                                                              >  2:305940/1‑151 (MQ=255)    cgcgCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGact                                                                                                                                           <  2:112972/151‑1 (MQ=255)       gcAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGactcc                                                                                                                                         >  1:251934/1‑150 (MQ=255)        cAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGactccac                                                                                                                                       <  1:98726/151‑1 (MQ=255)          gATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGactccactc                                                                                                                                     >  1:588375/1‑151 (MQ=255)           aTTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGactccactcc                                                                                                                                    >  2:384347/1‑151 (MQ=255)             tGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTccc                                                                                                                                   >  2:382450/1‑150 (MQ=255)                           tGATCGTCAGTACTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTgtgtgt                                                                                                                    >  1:611393/1‑151 (MQ=255)                                             cTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCAt                                                                                                  >  1:232535/1‑151 (MQ=255)                                                  accaACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtca                                                                                             <  1:345968/151‑1 (MQ=255)                                                       cATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGtcatcatc                                                                                        >  1:537090/1‑151 (MQ=255)                                                             tAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACt                                                                                  >  2:111818/1‑151 (MQ=255)                                                                 aCCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTgcgc                                                                              <  1:274754/151‑1 (MQ=255)                                                                 aCCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTgcgc                                                                              >  1:622532/1‑151 (MQ=255)                                                                  ccGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGcgcg                                                                             >  1:136718/1‑151 (MQ=255)                                                                    gTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTa                                                                           <  2:510206/151‑1 (MQ=255)                                                                      cGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTAt                                                                          <  1:603334/150‑1 (MQ=255)                                                                      cGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATg                                                                         >  1:550510/1‑151 (MQ=255)                                                                           atatCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGAt                                                                    >  2:97559/1‑151 (MQ=255)                                                                             atCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAg                                                                  <  1:343015/151‑1 (MQ=255)                                                                               cGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGtt                                                                >  2:659744/1‑151 (MQ=255)                                                                               cGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGt                                                                 <  1:60787/150‑1 (MQ=255)                                                                                 gAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTgg                                                              >  1:289887/1‑151 (MQ=255)                                                                                    ttGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCt                                                           >  1:165849/1‑151 (MQ=255)                                                                                    ttGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCt                                                           >  1:457023/1‑151 (MQ=255)                                                                                      gAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGc                                                         >  2:128168/1‑151 (MQ=255)                                                                                      gAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGc                                                         >  2:337312/1‑151 (MQ=255)                                                                                       aGGTAATCTTCGCCCAGCGCCATTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTg                                                          <  2:118562/149‑1 (MQ=255)                                                                                           aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGt                                                     <  2:587063/150‑1 (MQ=255)                                                                                           aaTCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTc                                                    >  1:591598/1‑151 (MQ=255)                                                                                                 cGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgcgc                                               <  2:528501/150‑1 (MQ=255)                                                                                                 cGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAgcgc                                               <  1:281293/150‑1 (MQ=255)                                                                                                  gCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGAGCGTATGTTGATAGTTGGTCTGCCAGTCAgcgcgc                                             >  1:628443/1‑151 (MQ=255)                                                                                                      aGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATAc                                         <  1:668284/151‑1 (MQ=255)                                                                                                          ccAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCa                                     <  2:504115/151‑1 (MQ=255)                                                                                                            aGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATg                                   >  1:245119/1‑151 (MQ=255)                                                                                                              tttGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGc                                 >  2:388325/1‑151 (MQ=255)                                                                                                               ttGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGc                                 <  2:85005/150‑1 (MQ=255)                                                                                                                     tttAATGACGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAg                            <  1:97561/149‑1 (MQ=255)                                                                                                                             cGGTAACGCCCGCCAGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCg                  <  2:107862/151‑1 (MQ=255)                                                                                                                                           aGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGcc    <  1:621387/151‑1 (MQ=255)                                                                                                                                           aGGCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGcc    >  1:542281/1‑151 (MQ=255)                                                                                                                                             gCATAATGGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGc  <  1:342933/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                AGCCGCGCAGATTGGGAGGCAACATTTGATCGTCAGTGCTGACACTGCCACCAACATAATTAAAACCGTCGAAAATATCGGAATTGAGGTAATCTTCGCCCAGCGCCAGTTTGGCTTTTAATGACGGTAACGCCCGCCAGGCATAATAGCGACTCCACTCCCATTTTTTTTGTGTGTCATCGCCGCCAAATTCATCGTCATCATCATTACTGCGCGTATGTTGATAGTTGGTCTGCCAGTCAGCGCGCATACGCCATGGCCCCAGGTTAACCCCGACCGTCCCGTTGCCGC  >  NC_000913/2454125‑2454415

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