Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	186,815	A→G	S612S (AGT→AGC)	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	186,815	0	A	G	100.0%	218.3 / NA	61	S612S (AGT→AGC)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/30);  total (31/30)

ATCCATGCGCAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGG  >  NC_000913/186692‑186940                                                                                                                            |                                                                                                                              aTCCATGCGCAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTa                                                                                                                >  2:280253/1‑139 (MQ=255)  tCCATGCGCAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAg                                                                                                                   <  1:204841/135‑1 (MQ=255)  tCCATGCGCAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAg                                                                                                                   >  2:204841/1‑135 (MQ=255)             tagtGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCt                                                                                                                             >  1:288950/1‑114 (MQ=255)             tagtGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGtt                                                                                                    >  2:1146189/1‑139 (MQ=255)                    agtgccTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtgg                                                                                             >  2:1509336/1‑139 (MQ=255)                    agtgccTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtgg                                                                                             <  1:1433565/139‑1 (MQ=255)                     gtgccTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtggt                                                                                            >  1:1522860/1‑139 (MQ=255)                          tGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGgcgc                                                                                       >  2:628293/1‑139 (MQ=255)                           ggAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCca                                                                                                           >  2:475926/1‑118 (MQ=255)                           ggAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCca                                                                                                           <  1:475926/118‑1 (MQ=255)                                 ggTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTc                                                                                >  1:1068539/1‑139 (MQ=255)                                  gTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCa                                                                               >  1:864767/1‑139 (MQ=255)                                         cGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAg                                                                        <  2:1050624/139‑1 (MQ=255)                                         cGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAg                                                                        <  1:851372/139‑1 (MQ=255)                                             cccGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    >  2:15395/1‑139 (MQ=255)                                             cccGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    <  1:1146189/139‑1 (MQ=255)                                             cccGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    <  1:209091/139‑1 (MQ=255)                                             cccGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    <  1:628293/139‑1 (MQ=255)                                                           tCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAgc                                                      >  2:965652/1‑139 (MQ=255)                                                                cGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    >  1:1606884/1‑120 (MQ=255)                                                                cGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    <  2:1606884/120‑1 (MQ=255)                                                                cGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCa                                                   >  1:25273/1‑137 (MQ=255)                                                                cGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCa                                                   <  2:25273/137‑1 (MQ=255)                                                                    ttttGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtgg                                                                                             >  1:104358/1‑91 (MQ=255)                                                                    ttttGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtgg                                                                                             <  2:104358/91‑1 (MQ=255)                                                                    ttttGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCca                                                                           >  1:369025/1‑109 (MQ=255)                                                                    ttttGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCca                                                                           <  2:369025/109‑1 (MQ=255)                                                                    ttttGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGAt                                             >  2:259587/1‑139 (MQ=255)                                                                         cATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTcc                                                                                                                     <  1:171390/62‑1 (MQ=255)                                                                         cATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTcc                                                                                                                     >  2:171390/1‑62 (MQ=255)                                                                            tCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAg                                                                                                                   >  1:659951/1‑61 (MQ=255)                                                                            tCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAg                                                                                                                   <  2:659951/61‑1 (MQ=255)                                                                            tCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGtt                                     >  1:1470413/1‑139 (MQ=255)                                                                              cgcgTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGc                                                                             <  1:706331/97‑1 (MQ=255)                                                                              cgcgTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGc                                                                             >  2:706331/1‑97 (MQ=255)                                                                              cgcgTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    >  2:157780/1‑106 (MQ=255)                                                                              cgcgTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    <  1:157780/106‑1 (MQ=255)                                                                              cgcgTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    <  1:427331/106‑1 (MQ=255)                                                                              cgcgTCGTAGCTGCTTTTCAGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    >  2:427331/1‑106 (MQ=255)                                                                                           ctTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCa                      <  1:259587/139‑1 (MQ=255)                                                                                             tttCGGTGGTAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCaaa                    <  2:1068539/139‑1 (MQ=255)                                                                                             tttCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAg                                                                                                                   <  1:536309/44‑1 (MQ=255)                                                                                             tttCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAg                                                                                                                   >  2:536309/1‑44 (MQ=255)                                                                                                      ctaaGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTtgat           <  1:980751/137‑1 (MQ=255)                                                                                                      cAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTtgat           <  2:1470413/139‑1 (MQ=255)                                                                                                      cAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTtgat           <  2:864767/139‑1 (MQ=255)                                                                                                             gagCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAg                                                                        <  1:116540/71‑1 (MQ=255)                                                                                                             gagCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAg                                                                        >  2:116540/1‑71 (MQ=255)                                                                                                               gCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    >  1:1654132/1‑73 (MQ=255)                                                                                                               gCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                    <  2:1654132/73‑1 (MQ=255)                                                                                                                cTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCgg  >  1:1191996/1‑138 (MQ=255)                                                                                                                cTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCgg  <  2:1191996/138‑1 (MQ=255)                                                                                                                   aCGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGt                                                                                 <  1:1611356/56‑1 (MQ=255)                                                                                                                   aCGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGt                                                                                 >  2:1611356/1‑56 (MQ=255)                                                                                                                     gCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGGTAGCGCAGACGATTTTCCGt                                      >  1:1197544/1‑97 (MQ=255)                                                                                                                     gCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGt                                      <  2:1197544/97‑1 (MQ=255)                                                                                                                     gCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCTCTTCTTCGGCaaa                    <  1:1120658/115‑1 (MQ=255)                                                                                                                     gCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCaaa                    >  2:1120658/1‑115 (MQ=255)                                                                                                                         cAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTcc                                        >  1:795477/1‑91 (MQ=255)                                                                                                                         cAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTcc                                        <  2:795477/91‑1 (MQ=255)                                                                                                                            |                                                                                                                              ATCCATGCGCAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGG  >  NC_000913/186692‑186940

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