Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	314,071	A→G	intergenic (‑266/‑286)	ykgP ← / → eaeH	pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology,factor, Not classified

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	314,071	0	A	G	100.0%	215.8 / NA	61	intergenic (‑266/‑286)	ykgP/eaeH	pseudogene, oxidoreductase family/pseudogene, attaching and effacing protein homology,factor, Not classified
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/33);  total (28/33)

ACAGTCCGATGAAGGAAAATGAATATCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTG  >  NC_000913/313943‑314205                                                                                                                                 |                                                                                                                                       aCAGTCCGATGAAGGAAAATGAATATCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACt                                                                                                                              >  1:648636/1‑139 (MQ=255)      ccGATGAAGGAAAATGAATATCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTGACTCTGAATAGTGGATTATTCACTGTCCg                                                                                                                         <  1:634485/139‑1 (MQ=255)      ccGATGAAGGAAAATGAATATCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCg                                                                                                                         <  2:648636/139‑1 (MQ=255)      ccGATGAAGGAAAATGAATATCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCg                                                                                                                         <  1:1265531/139‑1 (MQ=255)        gATGAAGGAAAATGAATATCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCa                                                                                                                       >  1:1436976/1‑139 (MQ=255)                       atatCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATaa                                                                                                        <  2:1140987/139‑1 (MQ=255)                         atCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACt                                                                                                      >  2:536127/1‑139 (MQ=255)                             aTAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAgg                                                                                                  <  2:787886/139‑1 (MQ=255)                                aaaaaTATATTGGTTTATCCTGGCATATATACCTATTTCGAGGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTAt                                                                                               >  2:944823/1‑139 (MQ=255)                                aaaaaTATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTAt                                                                                               <  2:1436976/139‑1 (MQ=255)                                aaaaaTATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAACTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTAt                                                                                               <  1:944823/139‑1 (MQ=255)                                       tatTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCg                                                                                                                         >  2:1374975/1‑106 (MQ=255)                                       tatTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCg                                                                                                                         <  1:1374975/106‑1 (MQ=255)                                            gTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGAc                                                                                   <  1:348420/139‑1 (MQ=255)                                               tATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTg                                                                                >  2:833575/1‑139 (MQ=255)                                                 tCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTAACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGca                                                                              >  2:235045/1‑139 (MQ=255)                                                      gCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAata                                                                                                           <  2:836309/105‑1 (MQ=255)                                                      gCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAata                                                                                                           >  1:836309/1‑105 (MQ=255)                                                      gCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTAt                                                                                               >  2:1299554/1‑117 (MQ=255)                                                      gCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTAt                                                                                               <  1:1299554/117‑1 (MQ=255)                                                      gCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACaaaa                                                                         >  2:491999/1‑139 (MQ=255)                                                      gCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACaaaa                                                                         >  1:985876/1‑139 (MQ=255)                                                           tataCCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTa                                                                    >  1:188919/1‑139 (MQ=255)                                                                cTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATc                                                                                              <  2:1498080/108‑1 (MQ=255)                                                                cTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATc                                                                                              >  1:1498080/1‑108 (MQ=255)                                                                   tttCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCg                                                                                                                         <  2:32455/78‑1 (MQ=255)                                                                   tttCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCg                                                                                                                         >  1:32455/1‑78 (MQ=255)                                                                   tttCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACt                                                                                                      <  1:768794/97‑1 (MQ=255)                                                                   tttCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACt                                                                                                      >  2:768794/1‑97 (MQ=255)                                                                   tttCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTcc                                                            >  2:1001644/1‑139 (MQ=255)                                                                    ttCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCa                                                                                                                       >  1:1400507/1‑79 (MQ=255)                                                                    ttCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCa                                                                                                                       <  2:1400507/79‑1 (MQ=255)                                                                    ttCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTccc                                                           >  1:54376/1‑139 (MQ=255)                                                                      cGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAAc                                                                                                       <  1:1084289/93‑1 (MQ=255)                                                                      cGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTccccc                                                         <  1:894251/139‑1 (MQ=255)                                                                      cGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATGCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAAc                                                                                                       >  2:1084289/1‑93 (MQ=255)                                                                                       ttttAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGt                                                                                                 >  1:1218105/1‑82 (MQ=255)                                                                                       ttttAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGt                                                                                                 <  2:1218105/82‑1 (MQ=255)                                                                                             ttaaAGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCcgcg                                  <  1:235045/139‑1 (MQ=255)                                                                                                aaGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGccc                               >  2:772784/1‑139 (MQ=255)                                                                                                 aGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTAt                                                                                               <  2:189400/74‑1 (MQ=255)                                                                                                 aGGCAGGTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTAt                                                                                               >  1:189400/1‑74 (MQ=255)                                                                                                       gTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCata                                                      <  2:344333/109‑1 (MQ=255)                                                                                                       gTCATTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCata                                                      >  1:344333/1‑109 (MQ=255)                                                                                                          aTTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCata                                                      <  1:1186774/106‑1 (MQ=255)                                                                                                          aTTGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCata                                                      >  2:1186774/1‑106 (MQ=255)                                                                                                            tGTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCtgatg                   <  1:536127/139‑1 (MQ=255)                                                                                                             gTTATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCtgatga                  <  2:188919/139‑1 (MQ=255)                                                                                                              ttATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTGGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATTAATAAAAATATAATCCCGCGCCCAACCACCTGATGAg                 <  2:421778/139‑1 (MQ=255)                                                                                                              ttATTCACTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAg                 <  2:54376/139‑1 (MQ=255)                                                                                                                    aCTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAAtt                                                                >  2:401126/1‑86 (MQ=255)                                                                                                                    aCTCTGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAAtt                                                                <  1:401126/86‑1 (MQ=255)                                                                                                                      tctGAATAGTGGATTATTCACTGTCCGCTGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCtgatga                  <  1:1429660/130‑1 (MQ=255)                                                                                                                      tctGAATAGTGGATTATTCACTGTCCGCTGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCtgatga                  >  2:1429660/1‑130 (MQ=255)                                                                                                                      tctGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCtat          >  1:1340140/1‑138 (MQ=255)                                                                                                                      tctGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCtat          <  2:1340140/138‑1 (MQ=255)                                                                                                                        tGAATAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAgg       >  2:89490/1‑139 (MQ=255)                                                                                                                          aaTAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCa     <  2:439971/139‑1 (MQ=255)                                                                                                                          aaTAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCa     <  1:89490/139‑1 (MQ=255)                                                                                                                          aaTAGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCa     <  2:1262553/139‑1 (MQ=255)                                                                                                                             aGTGGATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTg  <  2:609909/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                                       ACAGTCCGATGAAGGAAAATGAATATCCATAAAAAATATATTGGTTTATCCTGGCATATATACCTATTTCGACGTATTTCCAATAGTTTTAATTAAAGGCAGGTCATTGTTATTCACTCTGAATAGTGAATTATTCACTGTCCGCAGAGTAAGAAATATAACTTAGGTATCTATTTAATGACTTGCACAAAAAGCTAAATTTTCCCCCATAAATAAAAATATAATCCCGCGCCCAACCACCTGATGAGTGGCTATAGGCACTG  >  NC_000913/313943‑314205

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