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

AGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGC  >  NC_000913/3177891‑3178144                                                                                                                             |                                                                                                                                  aGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAAc                                                                                                                     >  2:357846/1‑139 (MQ=255)    ttGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCtt                                                                                                                  >  1:721482/1‑139 (MQ=255)             ttAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAg                                                                                                         >  1:1393938/1‑139 (MQ=255)                               cTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATcc                                                                                       >  1:983248/1‑139 (MQ=255)                               cTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATcc                                                                                       >  1:1323182/1‑139 (MQ=255)                                 ggCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCATATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGc                                                                                     >  1:390015/1‑139 (MQ=255)                                 ggCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGc                                                                                     >  2:231829/1‑139 (MQ=255)                                       aaTAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCaataat                                                                               >  1:695762/1‑139 (MQ=255)                                        aTAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAAtt                                                                              >  2:1184771/1‑139 (MQ=255)                                         tAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAAttt                                                                             <  1:1184771/139‑1 (MQ=255)                                          aGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAAtttt                                                                            <  1:231829/139‑1 (MQ=255)                                           gTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTa                                                                           <  2:983248/139‑1 (MQ=255)                                             gAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACa                                                                         >  1:1177489/1‑139 (MQ=255)                                                ttAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGtt                                                                      >  2:660503/1‑139 (MQ=255)                                                        gAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCATATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATcgcg                                                              <  2:390015/139‑1 (MQ=255)                                                             tCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTaaa                                                         <  1:660503/139‑1 (MQ=255)                                                              cAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAAt                                                        <  2:1393938/139‑1 (MQ=255)                                                               aGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     >  1:100002/1‑93 (MQ=255)                                                               aGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     <  2:100002/93‑1 (MQ=255)                                                                   aCGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     <  1:927834/89‑1 (MQ=255)                                                                   aCGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                     >  2:927834/1‑89 (MQ=255)                                                                     gTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCtt                                                 <  2:1323182/139‑1 (MQ=255)                                                                      tttGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAAttt                                                                             <  2:648125/110‑1 (MQ=255)                                                                      tttGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAAttt                                                                             >  1:648125/1‑110 (MQ=255)                                                                              ccGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGa                                                                   <  2:1118776/112‑1 (MQ=255)                                                                              ccGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGa                                                                   >  1:1118776/1‑112 (MQ=255)                                                                                 tttGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTcac                                              <  1:102670/130‑1 (MQ=255)                                                                                 tttGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTcac                                              >  2:102670/1‑130 (MQ=255)                                                                                         gTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACcaca                                          >  1:569071/1‑126 (MQ=255)                                                                                         gTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACcaca                                          <  2:569071/126‑1 (MQ=255)                                                                                          tCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATga                            <  2:1177489/139‑1 (MQ=255)                                                                                                    aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATa                                                       >  1:735828/1‑102 (MQ=255)                                                                                                    aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATa                                                       <  2:735828/102‑1 (MQ=255)                                                                                                    aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCt                  >  2:23898/1‑139 (MQ=255)                                                                                                         agtcacaTTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                 <  2:429961/84‑1 (MQ=255)                                                                                                         aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATGCGCAATAATTTTACAGTTTGATcgcg                                                              >  2:930513/1‑90 (MQ=255)                                                                                                         aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                 >  1:429961/1‑87 (MQ=255)                                                                                                         aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCa             >  1:1355737/1‑139 (MQ=255)                                                                                                         aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCa             >  1:630795/1‑139 (MQ=255)                                                                                                         aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCa             >  2:791222/1‑139 (MQ=255)                                                                                                         aGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTCTGATcgcg                                                              <  1:930513/90‑1 (MQ=255)                                                                                                                ttAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTAt      <  1:357846/139‑1 (MQ=255)                                                                                                                    cGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGc  <  2:630795/139‑1 (MQ=255)                                                                                                                    cGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGc  <  2:1355737/139‑1 (MQ=255)                                                                                                                    cGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGc  <  2:721482/139‑1 (MQ=255)                                                                                                                         tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATc                                                                                        >  2:532469/1‑48 (MQ=255)                                                                                                                         tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATc                                                                                        <  1:532469/48‑1 (MQ=255)                                                                                                                         tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGAt                                                                  <  2:157016/70‑1 (MQ=255)                                                                                                                         tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGAt                                                                  >  1:157016/1‑70 (MQ=255)                                                                                                                         tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACcaca                                          >  1:645183/1‑94 (MQ=255)                                                                                                                         tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACcaca                                          <  2:645183/94‑1 (MQ=255)                                                                                                                             |                                                                                                                                  AGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGC  >  NC_000913/3177891‑3178144

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: ‑