Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	187,607	A→G	R348R (CGT→CGC)	glnD ←	PII uridylyltransferase/uridylyl removing enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	187,607	0	A	G	100.0%	160.7 / NA	46	R348R (CGT→CGC)	glnD	PII uridylyltransferase/uridylyl removing enzyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/21);  total (25/21)

AATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCAC  >  NC_000913/187475‑187743                                                                                                                                     |                                                                                                                                         aaTGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGtt                                                                                                                                   >  2:82326/1‑140 (MQ=255)    gCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGttttt                                                                                                                                >  1:165907/1‑140 (MQ=255)      cGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTtcg                                                                                                                              <  2:833960/140‑1 (MQ=255)      cGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTtcg                                                                                                                              <  2:356150/140‑1 (MQ=255)                     tgtgCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCg                                                                                                               <  2:165907/140‑1 (MQ=255)                            cATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATgg                                                                                                        <  2:421012/140‑1 (MQ=255)                              tGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCt                                                                                                      >  2:567530/1‑140 (MQ=255)                              tGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCt                                                                                                      >  1:589651/1‑140 (MQ=255)                              tGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCt                                                                                                      >  1:255368/1‑140 (MQ=255)                                      aaCATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGaa                                                                                              >  1:948010/1‑140 (MQ=255)                                           aCGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGtt                                                                                         <  2:75070/140‑1 (MQ=255)                                           aCGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGtt                                                                                         >  1:624456/1‑140 (MQ=255)                                           aCGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGtt                                                                                         >  1:769320/1‑140 (MQ=255)                                           aCGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGtt                                                                                         >  1:97002/1‑140 (MQ=255)                                              cAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTgca                                                                                      <  2:769320/140‑1 (MQ=255)                                                aGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTgcagc                                                                                    <  2:488873/140‑1 (MQ=255)                                                 gATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTgcagca                                                                                   >  1:36176/1‑140 (MQ=255)                                                  aTGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAAc                                                                                             >  1:405619/1‑129 (MQ=255)                                                  aTGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAAc                                                                                             <  2:405619/129‑1 (MQ=255)                                                                               ttCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCgtg                                                     <  2:656778/140‑1 (MQ=255)                                                                                tcatcaCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGAcgcg                                                          >  2:973259/1‑134 (MQ=255)                                                                                tcatcaCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGAcgcg                                                          <  1:973259/134‑1 (MQ=255)                                                                                tcatcaCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCgtgt                                                    >  1:278453/1‑140 (MQ=255)                                                                                 catcaCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCa                                                                  >  1:340563/1‑125 (MQ=255)                                                                                 catcaCGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCa                                                                  <  2:340563/125‑1 (MQ=255)                                                                                           gTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGt                                         >  1:371110/1‑140 (MQ=255)                                                                                           gTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGt                                         >  2:795879/1‑140 (MQ=255)                                                                                              gATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAAt                                      <  2:842624/140‑1 (MQ=255)                                                                                              gATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAAt                                      >  1:842624/1‑140 (MQ=255)                                                                                                 tAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCt                                   <  2:36176/140‑1 (MQ=255)                                                                                                     gTACCGCGTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTtcat                               >  2:501415/1‑140 (MQ=255)                                                                                                         cgcgTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAAt                                      <  1:933201/129‑1 (MQ=255)                                                                                                         cgcgTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAAt                                      >  2:933201/1‑129 (MQ=255)                                                                                                         cgcgTAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCatcatc                           >  2:686541/1‑140 (MQ=255)                                                                                                             tAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCt                                   >  1:186353/1‑128 (MQ=255)                                                                                                             tAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCt                                   <  2:186353/128‑1 (MQ=255)                                                                                                             tAGCTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGAGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGcc                       >  1:4342/1‑139 (MQ=255)                                                                                                                cTGAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGa                    >  2:220246/1‑140 (MQ=255)                                                                                                                  gAAACTCATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGAcc                  <  2:447497/140‑1 (MQ=255)                                                                                                                        cATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCg            <  1:404238/140‑1 (MQ=255)                                                                                                                        cATCGTCGATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCg            <  2:97002/140‑1 (MQ=255)                                                                                                                          tcgtcgATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGtt          >  1:445002/1‑140 (MQ=255)                                                                                                                          tcgtcgATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGtt          <  2:445002/140‑1 (MQ=255)                                                                                                                          tcgtcgATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGtt          >  1:154542/1‑140 (MQ=255)                                                                                                                               gATTGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCtt     <  2:693250/140‑1 (MQ=255)                                                                                                                                  tGGGCGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTcac  <  2:278453/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                         AATGCCGGTGATCGCACTGTTGTGCACCATGGTGTAGAACATACGCAAGATGGCTTCCGGCTGGCGCATAAATAGTGTTTCATCACGCAGGTCGATTAGCGTACCGCGTAGCTGAAACTCATCGTCGATTGGACGTGGTTTTTCGTCGGCGGGAAGGGCGAGGATGGCTTCATCGAACAGTTGCAGCAGCATCTGGTTGAGTTCACTGACGCGGCGTGTAACGCGGAAGTAATCCTTCATCATCCGCTCGACCGGTTCGTTACCTTCAC  >  NC_000913/187475‑187743

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