Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,453,265	C→T	pseudogene (1737/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,453,265	0	C	T	100.0%	163.2 / NA	50	pseudogene (1737/2646 nt)	yfcU	pseudogene, FimD fimbrial export usher family,putative membrane, Not classified, putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/23);  total (27/23)

GCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAT  >  NC_000913/2453130‑2453395                                                                                                                                        |                                                                                                                                   gCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCGAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTAc                                                                                                                                 <  1:561015/139‑1 (MQ=255)         tCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCtt                                                                                                                         >  1:119858/1‑139 (MQ=255)                aTAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGcc                                                                                                                  >  2:709110/1‑139 (MQ=255)                      aGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTg                                                                                                            >  1:154285/1‑139 (MQ=255)                           cgTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTATCATa                                                                                                       >  2:785266/1‑139 (MQ=255)                               aTCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCg                                                                                                   >  2:643534/1‑139 (MQ=255)                                          tAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGcc                                                                                        <  2:119858/139‑1 (MQ=255)                                                   cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATCCTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGa                                                                               <  2:572214/139‑1 (MQ=255)                                                   cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCAt                                                                                                        >  1:92775/1‑114 (MQ=255)                                                   cGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCAt                                                                                                        <  2:92775/114‑1 (MQ=255)                                                         ggCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTATCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  1:785266/139‑1 (MQ=255)                                                              ctgTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCa                                                                                                         <  2:264556/102‑1 (MQ=255)                                                              ctgTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCa                                                                                                         >  1:264556/1‑102 (MQ=255)                                                               tgTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATAtt                                                                   >  2:829025/1‑139 (MQ=255)                                                                 tccgtccCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTaa                                                                 <  2:522740/139‑2 (MQ=255)                                                                           tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       <  1:643534/139‑1 (MQ=255)                                                                           tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       >  2:767108/1‑139 (MQ=255)                                                                           tGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       <  1:413729/139‑1 (MQ=255)                                                                             ccATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCa                                                     <  1:829025/139‑1 (MQ=255)                                                                                tAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCatat                                                  >  2:576733/1‑139 (MQ=255)                                                                                     gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGaaa                                             >  1:65327/1‑139 (MQ=255)                                                                                     gCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGaaa                                             <  2:65327/139‑1 (MQ=255)                                                                                                 aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGga                                     <  2:641907/135‑1 (MQ=255)                                                                                                 aaCGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGga                                     >  1:641907/1‑135 (MQ=255)                                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         >  1:212677/1‑97 (MQ=255)                                                                                                   cGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  2:212677/97‑1 (MQ=255)                                                                                                      tACTGTTGTCGCCCCACGGCATACTGTGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       <  1:407339/112‑1 (MQ=255)                                                                                                      tACTGTTGTCGCCCCACGGCATACTGTGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGcc                                                       >  2:407339/1‑112 (MQ=255)                                                                                                         tgttgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGAtgtt                         >  1:284479/1‑139 (MQ=255)                                                                                                          gttgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         <  2:256474/90‑1 (MQ=255)                                                                                                          gttgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCt                                                                         >  1:256474/1‑90 (MQ=255)                                                                                                            tgtCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTa                      >  1:311345/1‑139 (MQ=255)                                                                                                              tCGCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGt                    <  2:813648/139‑1 (MQ=255)                                                                                                                gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCat                                                    >  2:45996/1‑105 (MQ=255)                                                                                                                gCCCCACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCat                                                    <  1:45996/105‑1 (MQ=255)                                                                                                                  cccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCat                                                                      >  1:564998/1‑85 (MQ=255)                                                                                                                  cccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCat                                                                      <  2:564998/85‑1 (MQ=255)                                                                                                                   ccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGaaa                                             <  1:826780/109‑1 (MQ=255)                                                                                                                   ccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGaaa                                             >  2:826780/1‑109 (MQ=255)                                                                                                                   ccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGt                <  1:376072/138‑1 (MQ=255)                                                                                                                   ccACGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGt                >  2:376072/1‑138 (MQ=255)                                                                                                                     aCGGCATACTGAGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTg             >  1:142796/1‑139 (MQ=255)                                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  <  2:556477/139‑1 (MQ=255)                                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  >  1:787984/1‑139 (MQ=255)                                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  >  1:556477/1‑139 (MQ=255)                                                                                                                                aGCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAt  >  1:400699/1‑139 (MQ=255)                                                                                                                                 gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgctgc                                                        >  1:470274/1‑84 (MQ=255)                                                                                                                                 gCGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAAtgctgc                                                        <  2:470274/84‑1 (MQ=255)                                                                                                                                  cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCa      <  2:136431/133‑1 (MQ=255)                                                                                                                                  cGAAATTTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCa      >  1:136431/1‑133 (MQ=255)                                                                                                                                        |                                                                                                                                   GCCGATGTTCAACTGATAGTGAGTCGCGTCATCGACACGGCTAAAATAACCGACCTGGCTGCTGTCCGTCCCACTGCCATAGTTGCCGTTATAGCTAACGGTACTGTTGTCGCCCCACGGCATACTGAGCGAAATCTACATGCCCTTATCCGCCCGGTTGTCATACTCGTAGCGGTAGCCAGTCAGGGAAACGCTCATATTACGAATGCTGCCCATATTGAAATAGTGGGAGAGCATGATGTTGTAGTTTGTCTGCTCCTCACGAT  >  NC_000913/2453130‑2453395

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