Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,178,015	G→A	intergenic (‑105/‑100)	nudF ← / → tolC	ADP‑ribose pyrophosphatase/transport channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,178,015	0	G	A	100.0%	160.0 / NA	50	intergenic (‑105/‑100)	nudF/tolC	ADP‑ribose pyrophosphatase/transport channel
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/27);  total (23/27)

GCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGC  >  NC_000913/3177889‑3178144                                                                                                                               |                                                                                                                                  gCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCa                                                                                                                       >  2:288453/1‑139 (MQ=255)   aGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAAc                                                                                                                     >  1:935503/1‑139 (MQ=255)               ttAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAg                                                                                                         <  1:970649/139‑1 (MQ=255)                  aGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCgg                                                                                                      >  1:881328/1‑139 (MQ=255)                   gCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     <  2:881328/139‑1 (MQ=255)                   gCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     <  2:804389/139‑1 (MQ=255)                   gCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     <  2:1315053/139‑1 (MQ=255)                       tGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCt                                                                                                 >  2:1027746/1‑139 (MQ=255)                           aaTGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGt                                                                                                                              >  2:682277/1‑106 (MQ=255)                           aaTGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGt                                                                                                                              <  1:682277/106‑1 (MQ=255)                             tGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGa                                                                                           <  2:935503/139‑1 (MQ=255)                                 cTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATcc                                                                                       >  2:364323/1‑139 (MQ=255)                                          aTAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGcc                                                                                                                        <  2:1284565/97‑1 (MQ=255)                                          aTAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGcc                                                                                                                        >  1:1284565/1‑97 (MQ=255)                                              tGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     >  1:502816/1‑112 (MQ=255)                                              tGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     <  2:502816/112‑1 (MQ=255)                                              tGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGa                                                                                           >  1:387754/1‑122 (MQ=255)                                              tGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGa                                                                                           <  2:387754/122‑1 (MQ=255)                                                aaTTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAg                                                                        >  2:106551/1‑139 (MQ=255)                                                       tgtgAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                 <  1:106551/139‑1 (MQ=255)                                                          gAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATcgcg                                                              <  1:288453/139‑1 (MQ=255)                                                           aaTTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCgcgc                                                             >  1:645348/1‑139 (MQ=255)                                                            aTTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCt                                                            >  2:1033551/1‑139 (MQ=255)                                                                   cGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGc                                                                                                        <  2:8086/88‑1 (MQ=255)                                                                   cGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGc                                                                                                        >  1:8086/1‑88 (MQ=255)                                                                     aCGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGc                                                   <  2:658842/139‑1 (MQ=255)                                                                     aCGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGc                                                   >  1:173920/1‑139 (MQ=255)                                                                      cGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCt                                                  <  1:364323/139‑1 (MQ=255)                                                                       gTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCtt                                                 <  1:756652/139‑1 (MQ=255)                                                                       gTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCtt                                                 <  2:173920/139‑1 (MQ=255)                                                                              tGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACcaca                                          <  1:1033551/139‑1 (MQ=255)                                                                                ccGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCaa                                                                                   >  2:542879/1‑96 (MQ=255)                                                                                ccGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCaa                                                                                   <  1:542879/96‑1 (MQ=255)                                                                                 cGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACaa                                         >  2:1281287/1‑137 (MQ=255)                                                                                 cGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACaa                                         <  1:1281287/137‑1 (MQ=255)                                                                                  gTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAg                                                                                                         >  2:1125005/1‑72 (MQ=255)                                                                                  gTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAg                                                                                                         <  1:1125005/72‑1 (MQ=255)                                                                                                      aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGa                                                                   <  1:491868/90‑1 (MQ=255)                                                                                                      aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGa                                                                   >  2:491868/1‑90 (MQ=255)                                                                                                      aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAAt                                                        >  2:723772/1‑101 (MQ=255)                                                                                                      aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAAt                                                        <  1:723772/101‑1 (MQ=255)                                                                                                           aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                 <  2:515687/87‑1 (MQ=255)                                                                                                           aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                 <  2:1028368/87‑1 (MQ=255)                                                                                                           aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                 >  1:515687/1‑87 (MQ=255)                                                                                                           aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                 >  1:1028368/1‑87 (MQ=255)                                                                                                                caTTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCtt        >  1:785516/1‑139 (MQ=255)                                                                                                                      cGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCt                  >  2:315741/1‑123 (MQ=255)                                                                                                                      cGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCt                  <  1:315741/123‑1 (MQ=255)                                                                                                                      cGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGc  <  2:645348/139‑1 (MQ=255)                                                                                                                      cGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGc  <  1:1027746/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                  GCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGC  >  NC_000913/3177889‑3178144

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