Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,026,030	G→A	E89K (GAA→AAA)	guaD →	guanine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,026,030	0	G	A	100.0%	221.1 / NA	62	E89K (GAA→AAA)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/34);  total (28/34)

AACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAA  >  NC_000913/3025893‑3026163                                                                                                                                          |                                                                                                                                      aaCAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTaaa                                                                                                                                     >  1:193248/1‑140 (MQ=255)     gggAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtgg                                                                                                                                 >  1:702846/1‑140 (MQ=255)     gggAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtgg                                                                                                                                 >  1:1136005/1‑140 (MQ=255)      ggAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggt                                                                                                                                >  1:460167/1‑140 (MQ=255)      ggAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggt                                                                                                                                >  1:179531/1‑140 (MQ=255)       gAAATGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggtg                                                                                                                               <  2:1126821/140‑1 (MQ=255)       gAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggtg                                                                                                                               <  1:1000216/140‑1 (MQ=255)       gAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggtg                                                                                                                               >  2:1189028/1‑140 (MQ=255)       gAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggtg                                                                                                                               <  2:262696/140‑1 (MQ=255)       gAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggtg                                                                                                                               <  2:589545/140‑1 (MQ=255)       gAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggtg                                                                                                                               <  1:1189028/140‑1 (MQ=255)       gAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAAtggtg                                                                                                                               <  1:1074532/140‑1 (MQ=255)              ggAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCt                                                                                                                        <  2:460167/140‑1 (MQ=255)              ggAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCt                                                                                                                        >  2:405587/1‑140 (MQ=255)              ggAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCt                                                                                                                        <  1:977591/140‑1 (MQ=255)               gAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTa                                                                                                                       <  2:658936/140‑1 (MQ=255)                  tGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATgg                                                                                                                    <  1:49116/140‑1 (MQ=255)                   ggTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGt                                                                                                                   <  1:988335/140‑1 (MQ=255)                       tGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGc                                                                                                               >  2:1102352/1‑140 (MQ=255)                       tGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGc                                                                                                               <  1:405587/140‑1 (MQ=255)                       tGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGc                                                                                                               <  1:1102352/140‑1 (MQ=255)                       tGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGc                                                                                                               <  2:193248/140‑1 (MQ=255)                       tGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGc                                                                                                               <  1:1047451/140‑1 (MQ=255)                              tGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCt                                                                                                        >  1:928309/1‑140 (MQ=255)                              tGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCt                                                                                                        >  2:224704/1‑140 (MQ=255)                              tGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCt                                                                                                        >  2:344445/1‑140 (MQ=255)                                 gAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGa                                                                                                     <  1:3910/140‑1 (MQ=255)                                  aaaaCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAg                                                                                                    <  1:54321/140‑1 (MQ=255)                                      cGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGt                                                                                                <  2:1136005/140‑1 (MQ=255)                                           aGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAAt                                                                                           >  1:1010126/1‑140 (MQ=255)                                           aGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAAt                                                                                           >  1:508532/1‑140 (MQ=255)                                                cAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAaca                                                                                      <  2:928309/140‑1 (MQ=255)                                                 aaaTTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacac                                                                                     <  1:250797/140‑1 (MQ=255)                                                   aTTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacacac                                                                                   <  2:702846/140‑1 (MQ=255)                                                   aTTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacacac                                                                                   <  2:508532/140‑1 (MQ=255)                                                   aTTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacacac                                                                                   >  1:23271/1‑140 (MQ=255)                                                               aTTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAaca                                                                                      >  1:847144/1‑125 (MQ=255)                                                               aTTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTACCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAaca                                                                                      <  2:847144/125‑1 (MQ=255)                                                                 tCGCGTGCGCGACTATCGCGGCAAACTGATAGTGCCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAac                                                                                       >  2:780471/1‑122 (MQ=255)                                                                 tCGCGTGCGCGACTATCGCGGCAAACTGATAGTGCCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAac                                                                                       <  1:780471/122‑1 (MQ=255)                                                                    cgTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgtc                                                                  >  1:525269/1‑140 (MQ=255)                                                                         gcgACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTAt                                                             <  1:344445/140‑1 (MQ=255)                                                                                          cTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTa                                                              <  2:726984/122‑1 (MQ=255)                                                                                          cTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTa                                                              >  1:726984/1‑122 (MQ=255)                                                                                                     gggCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGt                                 >  1:719148/1‑140 (MQ=255)                                                                                                       gCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCg                               <  2:179531/140‑1 (MQ=255)                                                                                                       gCTTTGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCg                               >  2:438242/1‑140 (MQ=255)                                                                                                           tGTCGATACACATATCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGt                           >  1:8176/1‑140 (MQ=255)                                                                                                                      atatCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAag                <  2:8176/140‑1 (MQ=255)                                                                                                                      atatCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAag                <  1:224704/140‑1 (MQ=255)                                                                                                                       tatCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAagc               >  2:174814/1‑140 (MQ=255)                                                                                                                       tatCCATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAagc               <  2:1207751/140‑1 (MQ=255)                                                                                                                           cATTATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCt           >  2:1183018/1‑140 (MQ=255)                                                                                                                              tATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGtt                          <  2:821011/122‑1 (MQ=255)                                                                                                                              tATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGtt                          >  1:821011/1‑122 (MQ=255)                                                                                                                              tATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCtttt        >  1:12174/1‑140 (MQ=255)                                                                                                                              tATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCtttt        >  2:33910/1‑140 (MQ=255)                                                                                                                              tATCCGCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCtttt        >  1:384476/1‑140 (MQ=255)                                                                                                                                   gCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAag                <  1:1111649/127‑1 (MQ=255)                                                                                                                                   gCAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAag                >  2:1111649/1‑127 (MQ=255)                                                                                                                                    cAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTaa  <  2:719148/140‑1 (MQ=255)                                                                                                                                    cAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTaa  <  2:525269/140‑1 (MQ=255)                                                                                                                                    cAAAGTAAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTaa  >  2:147457/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                      AACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAA  >  NC_000913/3025893‑3026163

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