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%	114.8 / NA	33	S612S (AGT→AGC)	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/15);  total (18/15)

GCGCAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTCC  >  NC_000913/186698‑186944                                                                                                                      |                                                                                                                                  gcgcAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAg                                                                                                                       >  2:408934/1‑130 (MQ=255) gcgcAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAg                                                                                                                       <  1:408934/130‑1 (MQ=255)          tgccagtgccTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCg                                                                                                     >  1:418167/1‑139 (MQ=255)           gccagtgccTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGt                                                                                                    >  2:719439/1‑139 (MQ=255)            ccagtgtcTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGtt                                                                                                   >  2:522662/1‑139 (MQ=255)               gtgccTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGc                                                                                                                      <  2:470455/117‑1 (MQ=255)               gtgccTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGc                                                                                                                      >  1:470455/1‑117 (MQ=255)                 gccTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtggtgg                                                                                              >  2:480648/1‑139 (MQ=255)                    tGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCtt                                                                                                                           <  1:765244/107‑1 (MQ=255)                    tGGAGTTGGTGATGGCGAACCCGTACGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCtt                                                                                                                           >  2:765244/1‑107 (MQ=255)                         ttGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATg                                                                                      <  1:773780/139‑1 (MQ=255)                           ggTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTc                                                                                    >  2:242067/1‑139 (MQ=255)                            gTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCa                                                                                   >  2:399364/1‑139 (MQ=255)                                 ggCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCGCAGCGTTTCGTTGGTGGCGCAAATGTCAGCcac                                                                              >  1:146221/1‑139 (MQ=255)                                    gAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGt                                                                           <  1:399364/139‑1 (MQ=255)                                       cccGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCGCAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAg                                                                        <  2:146221/139‑1 (MQ=255)                                        ccGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAgg                                                                       <  1:522662/139‑1 (MQ=255)                                            tCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGt                                                                                     <  2:523127/121‑1 (MQ=255)                                            tCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGt                                                                                     >  1:523127/1‑121 (MQ=255)                                                          cGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGa                                                     >  1:423406/1‑139 (MQ=255)                                                                       ccGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGttt                                        >  1:351050/1‑139 (MQ=255)                                                                        cgcgTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTg                                       >  1:839731/1‑139 (MQ=255)                                                                         gcgTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGc                                      <  2:423406/139‑1 (MQ=255)                                                                                  ctgctTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtgg                                                                                                 <  1:381093/71‑1 (MQ=255)                                                                                  ctgctTTTCGGTGGCAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtgg                                                                                                 >  2:381093/1‑71 (MQ=255)                                                                                                cAAAGTAGAGCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTtgat               <  2:839731/139‑1 (MQ=255)                                                                                                       gagCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTcc        <  1:242067/139‑1 (MQ=255)                                                                                                         gCTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCgg      <  2:150953/139‑1 (MQ=255)                                                                                                          cTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAg                                                                            <  1:762506/68‑1 (MQ=255)                                                                                                          cTCACGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAg                                                                            >  2:762506/1‑68 (MQ=255)                                                                                                             aCGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtggt                                                                                                <  2:211997/45‑1 (MQ=255)                                                                                                             aCGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTtggt                                                                                                >  1:211997/1‑45 (MQ=255)                                                                                                             aCGCAACAGGCTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTcc  >  1:778958/1‑139 (MQ=255)                                                                                                                      |                                                                                                                                  GCGCAGTAGTGCCAGTGCCTGGAGTTGGTGATGGCGAACCCGTTCGCGCATATCCGGCGTGTTTTGCATTCCGCGTCGTAGCTGCTTTTCGGTGGCAAAGTAGAGCTCACGCAACAGACTTTGCTTCCAGCTATTCCACAGCGTTTCGTTGGTGGCGCAAATGTCAGCCACAGTCAGGCATACCAGATAGCGCAGACGATTTTCCGTTTGCACTTCTTCGGCAAACTGCTTGATGACTTCCGGGTCC  >  NC_000913/186698‑186944

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