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%	131.1 / NA	41	intergenic (‑105/‑100)	nudF/tolC	ADP‑ribose pyrophosphatase/transport channel
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/22);  total (19/22)

CGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAG  >  NC_000913/3177886‑3178149                                                                                                                                  |                                                                                                                                       cGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACg                                                                                                                               >  1:115035/1‑139 (MQ=255)            tCTGGCTTAAGCATTGTTAATGTCCTGGCTCTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTGTTg                                                                                                                    <  2:692692/139‑1 (MQ=255)                  ttAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAg                                                                                                              <  1:351035/139‑1 (MQ=255)                  ttAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAg                                                                                                              <  1:434612/139‑1 (MQ=255)                   tAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGACTTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGc                                                                                                             >  2:761961/1‑139 (MQ=255)                      gCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                          <  1:761961/139‑1 (MQ=255)                          tGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCt                                                                                                      >  2:232963/1‑139 (MQ=255)                                tGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCaa                                                                                                                           <  2:534211/112‑1 (MQ=255)                                tGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCaa                                                                                                                           >  1:534211/1‑112 (MQ=255)                                     tGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCg                                                                                           >  1:115526/1‑139 (MQ=255)                                      ggCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGc                                                                                          >  2:729139/1‑139 (MQ=255)                                            aaTAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCaataat                                                                                    >  2:210887/1‑139 (MQ=255)                                                   aaTTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAAcc                                                                                                                         <  2:801012/95‑1 (MQ=255)                                                   aaTTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAAcc                                                                                                                         >  1:801012/1‑95 (MQ=255)                                                       aaaTGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTg                                                                         >  2:950695/1‑139 (MQ=255)                                                          tgtgAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                      <  1:950695/139‑1 (MQ=255)                                                             gAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATcgcg                                                                   <  1:439806/139‑1 (MQ=255)                                                               aTTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGaa                                                                                               >  2:514048/1‑109 (MQ=255)                                                               aTTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGaa                                                                                               <  1:514048/109‑1 (MQ=255)                                                               aTTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATGATTTTACAGTTTGATCGCGCt                                                                 <  2:598895/139‑1 (MQ=255)                                                                  tCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTaaa                                                              <  1:861540/139‑1 (MQ=255)                                                                           tttGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTc                                                     >  2:531039/1‑139 (MQ=255)                                                                             tGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTcac                                                   <  2:115526/139‑1 (MQ=255)                                                                              gACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTcacc                                                  <  1:232963/139‑1 (MQ=255)                                                                                   ccGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAg                                             <  2:406020/139‑1 (MQ=255)                                                                                          aGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCa                                      >  2:230718/1‑139 (MQ=255)                                                                                               tCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAgg                                            >  2:760785/1‑128 (MQ=255)                                                                                               tCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAgg                                            <  1:760785/128‑1 (MQ=255)                                                                                                         aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCt                       <  1:431595/139‑1 (MQ=255)                                                                                                         aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCt                       <  1:230718/139‑1 (MQ=255)                                                                                                           ttGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                      <  2:1048873/90‑1 (MQ=255)                                                                                                           ttGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                      >  1:1048873/1‑90 (MQ=255)                                                                                                           ttGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTcc                     <  1:636092/139‑1 (MQ=255)                                                                                                              aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCa                  >  1:829327/1‑139 (MQ=255)                                                                                                              aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCa                  >  2:985395/1‑139 (MQ=255)                                                                                                                gCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTg                         <  2:482097/130‑1 (MQ=255)                                                                                                                gCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTg                         >  1:482097/1‑130 (MQ=255)                                                                                                                   caTTAACGCCCTATGACACGTAACGCCAACCTTTTg                                                                                                                    <  2:237729/36‑1 (MQ=255)                                                                                                                   caTTAACGCCCTATGACACGTAACGCCAACCTTTTg                                                                                                                    >  1:237729/1‑36 (MQ=255)                                                                                                                              tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAg  <  2:287303/139‑1 (MQ=255)                                                                                                                              tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAg  >  1:287303/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                       CGGGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAG  >  NC_000913/3177886‑3178149

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