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%	161.7 / 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)

GGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAG  >  NC_000913/3177888‑3178149                                                                                                                                |                                                                                                                                       ggCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGcc                                                                                                                             <  1:1070505/139‑1 (MQ=255)   cAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCaa                                                                                                                           >  1:23875/1‑139 (MQ=255)     ggTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAAcc                                                                                                                         >  1:882997/1‑139 (MQ=255)        tGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCttt                                                                                                                      <  2:882997/139‑1 (MQ=255)         gTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCtttt                                                                                                                     >  1:1189592/1‑139 (MQ=255)                ttAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAg                                                                                                              <  2:1189592/139‑1 (MQ=255)                  aaGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACTTTTTGCGGTAGCg                                                                                                            <  2:1030023/139‑1 (MQ=255)                  aaGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCg                                                                                                            <  2:333173/139‑1 (MQ=255)                   aGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCgg                                                                                                           >  1:564052/1‑139 (MQ=255)                    gCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                          <  1:586652/139‑1 (MQ=255)                    gCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                          <  2:23875/139‑1 (MQ=255)                        tGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGcc                                                                                                                             >  1:201844/1‑116 (MQ=255)                        tGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGcc                                                                                                                             <  2:201844/116‑1 (MQ=255)                        tGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCa                                                                                                                            <  2:1141662/117‑1 (MQ=255)                        tGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCa                                                                                                                            >  1:1141662/1‑117 (MQ=255)                                  cTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATcc                                                                                            >  2:475989/1‑139 (MQ=255)                                    ggCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGc                                                                                          >  2:836683/1‑139 (MQ=255)                                    ggCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGc                                                                                          >  1:587436/1‑139 (MQ=255)                                    ggCACATATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTTATGTGTTATATTGGGGCCCATAAAGGCCCTATGACAGTTATCTCCATCTTTTTGCGGTAGCGGCTTCTGCTAGAATCCGc                                                                                          >  2:1062898/1‑139 (MQ=255)                                         tAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCaataa                                                                                     >  1:1495052/1‑139 (MQ=255)                                          aaTAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCaataat                                                                                    >  1:20157/1‑139 (MQ=255)                                            tAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAAttt                                                                                  <  1:475989/139‑1 (MQ=255)                                              gTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTa                                                                                <  2:1495052/139‑1 (MQ=255)                                               tGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTAc                                                                               >  2:556973/1‑139 (MQ=255)                                                 aaTTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAg                                                                             >  1:738306/1‑139 (MQ=255)                                                      aaTGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGa                                                                        <  2:1067841/139‑1 (MQ=255)                                                        tgtgAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                      <  1:1462314/139‑1 (MQ=255)                                                        tgtgAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                      <  2:20157/139‑1 (MQ=255)                                                           gAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATcgcg                                                                   >  2:957934/1‑139 (MQ=255)                                                                  aGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATa                                                            >  1:65158/1‑139 (MQ=255)                                                                       cGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCt                                                       <  1:957934/139‑1 (MQ=255)                                                                         tttGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTc                                                     <  1:836683/139‑1 (MQ=255)                                                                             aCTGCCGTTTGAGCAGTAATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTcacca                                                 <  1:404745/139‑1 (MQ=255)                                                                                 ccGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAg                                             <  2:738306/139‑1 (MQ=255)                                                                                 ccGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACAACAAg                                             >  1:1137948/1‑139 (MQ=255)                                                                                           aGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTg                                                                                                                    <  2:813065/58‑1 (MQ=255)                                                                                           aGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTg                                                                                                                    >  1:813065/1‑58 (MQ=255)                                                                                              cATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCaataa                                                                                     >  1:42309/1‑86 (MQ=255)                                                                                              cATGTGTTAAATTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCaataa                                                                                     <  2:42309/86‑1 (MQ=255)                                                                                                       aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                          <  1:598483/56‑1 (MQ=255)                                                                                                       aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGc                                                                                                          >  2:598483/1‑56 (MQ=255)                                                                                                       aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAAttt                                                                                  <  1:347406/80‑1 (MQ=255)                                                                                                       aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAAttt                                                                                  >  2:347406/1‑80 (MQ=255)                                                                                                       aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                      <  2:287297/92‑1 (MQ=255)                                                                                                       aaTTGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATc                                                                      >  1:287297/1‑92 (MQ=255)                                                                                                         ttGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCt                                                                 <  2:808162/95‑1 (MQ=255)                                                                                                         ttGAGGCACATTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCt                                                                 >  1:808162/1‑95 (MQ=255)                                                                                                                 caTTAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCtt             >  1:334330/1‑139 (MQ=255)                                                                                                                   ttAACGCCCTATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTAt           >  2:793558/1‑139 (MQ=255)                                                                                                                            tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAg  >  2:1435022/1‑139 (MQ=255)                                                                                                                            tATGACACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAg  >  2:342131/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                       GGCAGGTTGTCTGGCTTAAGCATTGTTAATGTCCTGGCACTAATAGTGAATTAAATGTGAATTTCAGCGACGTTTGACTGCCGTTTGAGCAGTCATGTGTTAAATTGAGGCACATTAACGCCCTATGGCACGTAACGCCAACCTTTTGCGGTAGCGGCTTCTGCTAGAATCCGCAATAATTTTACAGTTTGATCGCGCTAAATACTGCTTCACCACAAGGAATGCAAATGAAGAAATTGCTCCCCATTCTTATCGGCCTGAG  >  NC_000913/3177888‑3178149

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