Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,683,472	G→T	intergenic (‑25/+158)	dctA ← / ← yhjK	C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,683,472	0	G	T	100.0%	123.8 / NA	36	intergenic (‑25/+158)	dctA/yhjK	C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/19);  total (17/19)

AGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGT  >  NC_000913/3683327‑3683613                                                                                                                                                  |                                                                                                                                              aGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTc                                                                                                                                            <  6:79777/149‑1 (MQ=255)               ggTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATgg                                                                                                                              <  6:104838/149‑1 (MQ=255)                                    tCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAAt                                                                                                          <  2:120851/148‑1 (MQ=255)                                          tAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                    >  2:75973/1‑148 (MQ=255)                                           aGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                    >  1:78477/1‑147 (MQ=255)                                           aGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAcc                                                                                                    >  3:124880/1‑147 (MQ=255)                                               aTGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGac                                                                                               >  7:12918/1‑148 (MQ=255)                                                tGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAAAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGacaa                                                                                             <  5:62860/149‑1 (MQ=255)                                                     aaGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATca                                                                                        >  6:153946/1‑149 (MQ=255)                                                       ggAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcaca                                                                                      <  2:78477/149‑1 (MQ=255)                                                        gagaATACCAAAGGCTATCTCTGTCAGGACATGAAAGTAAAGGCTTTTAAACAGAGAGTTTTTAATAGGGTGTCCTTTAGTAAAACCAAATGTATTGTAAGGTTATGGGGTACCTGCTGCCTTAAAATAAAACCAAGACAACATca                                                                                        >  8:7673/1‑146 (MQ=255)                                                               ccAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGt                                                                               <  8:12918/148‑1 (MQ=255)                                                                cAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAAt                                                                                 >  6:149340/1‑145 (MQ=255)                                                                cAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAAt                                                                                 <  5:149340/145‑1 (MQ=255)                                                                cAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTAc                                                                             >  2:84747/1‑149 (MQ=255)                                                                cAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTAc                                                                             <  5:153946/149‑1 (MQ=255)                                                                cAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTAc                                                                             >  8:12955/1‑149 (MQ=255)                                                                  aTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCt                                                                           <  7:7673/149‑1 (MQ=255)                                                                     gCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGa                                                                        <  1:75973/149‑1 (MQ=255)                                                                               tCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAtt                                                                <  1:84747/147‑1 (MQ=255)                                                                                             tAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggg                                                                                                                             >  7:141426/1‑72 (MQ=255)                                                                                             tAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATggg                                                                                                                             <  8:141426/72‑1 (MQ=255)                                                                                                  gCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAAt                                                                                 <  5:83004/111‑1 (MQ=255)                                                                                                  gCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAAt                                                                                 >  6:83004/1‑111 (MQ=255)                                                                                                     tttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                         >  3:94998/1‑148 (MQ=255)                                                                                                        aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                         <  7:131563/145‑1 (MQ=255)                                                                                                        aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                         >  8:131563/1‑145 (MQ=255)                                                                                                                              tgtCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAtt                                                                >  7:93384/1‑100 (MQ=255)                                                                                                                              tgtCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAtt                                                                <  8:93384/100‑1 (MQ=255)                                                                                                                                tCCTTTAGTAAAACCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttt              <  4:124880/148‑1 (MQ=255)                                                                                                                                  cTTTAGTAAAACCACATGTCTTGAAAGGTTATGGGGTACCTGCGGCCTTAAAATAACAACCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAACGTTAAAAAGTTTGAGCAGGCTCGCACAAACCAGCACttttt            >  3:41698/1‑148 (MQ=255)                                                                                                                                          aaaCCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg   <  7:12955/149‑1 (MQ=255)                                                                                                                                          aaaCCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg   >  2:178447/1‑149 (MQ=255)                                                                                                                                           aaCCACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTTAAATAACACGCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGt  <  4:41698/149‑1 (MQ=255)                                                                                                                                             ccACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAg                    >  4:32694/1‑129 (MQ=255)                                                                                                                                             ccACATGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAg                    <  3:32694/129‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              AGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGT  >  NC_000913/3683327‑3683613

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