Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,103,388	C→T	H189H (CAC→CAT)	yiiM →	2‑amino‑6‑N‑hydroxylaminopurine resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,103,388	0	C	T	100.0%	106.7 / NA	34	H189H (CAC→CAT)	yiiM	2‑amino‑6‑N‑hydroxylaminopurine resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/21);  total (13/21)

GTTGTATAGCGTGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTG  >  NC_000913/4103253‑4103515                                                                                                                                        |                                                                                                                                gtTGTATAGCGTGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGtt                                                                                                                             <  2:810544/140‑1 (MQ=255)   tGTATAGCGTGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTa                                                                                                                           >  2:688153/1‑140 (MQ=255)        aGCGTGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACtctc                                                                                                                      <  1:527562/140‑1 (MQ=255)            tGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGTAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCt                                                                                                                  >  2:29332/1‑140 (MQ=255)            tGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCt                                                                                                                  >  2:91100/1‑140 (MQ=255)            tGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCt                                                                                                                  >  2:596329/1‑140 (MQ=255)            tGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATACCGTTCGATGACGATCTATATCATCGTTTActccgctgc                                                                                                                  >  2:348522/1‑134 (MQ=255)             gATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTg                                                                                                                 <  2:20018/140‑1 (MQ=255)                              tgtCTGCGGACGCGCTGCTGGACTTGGTTTCCCGTGTGAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGc                                                                                                <  1:270148/140‑1 (MQ=255)                                  tGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATAATCCTTTACTCTCCGCCGCCGTGTTATCGAAAAGCTggc                                                                                            >  2:270148/1‑139 (MQ=255)                                         gcgcCGCTGGAGTTGTTTTCCCGTGTCAGAGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTA‑‑CTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGAcct                                                                                   <  1:348522/140‑3 (MQ=255)                                                tGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCaa                                                                              <  2:798701/140‑1 (MQ=255)                                                     ttGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAgcg                                                                         >  2:576291/1‑140 (MQ=255)                                                          ttCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGa                                                                                      >  1:521086/1‑122 (MQ=255)                                                          ttCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGa                                                                                      <  2:521086/122‑1 (MQ=255)                                                                   cAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCa                                                           <  1:91100/140‑1 (MQ=255)                                                                    aGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCaa                                                          >  2:618050/1‑140 (MQ=255)                                                                           tgACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGaa                                                   >  1:422029/1‑140 (MQ=255)                                                                                 tGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGAtttt                                             >  2:495734/1‑140 (MQ=255)                                                                                   cAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAg                                                  <  1:525612/133‑1 (MQ=255)                                                                                   cAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAg                                                  >  2:525612/1‑133 (MQ=255)                                                                                       aggCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTCCTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCccgc                                       <  1:817585/140‑1 (MQ=255)                                                                                           tGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGt                                   >  2:469893/1‑140 (MQ=255)                                                                                                       aTGGCATATGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAAt                       <  1:688153/140‑1 (MQ=255)                                                                                                            atatGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAAccc                  <  1:596329/140‑1 (MQ=255)                                                                                                               tGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCg               <  1:495734/140‑1 (MQ=255)                                                                                                               tGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCg               <  1:508348/140‑1 (MQ=255)                                                                                                               tGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCg               <  1:806923/140‑1 (MQ=255)                                                                                                               tGCCGTTTGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGAGAAAATAACCCCCg               <  1:544041/140‑1 (MQ=255)                                                                                                                    tttGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAgtg          <  1:469893/140‑1 (MQ=255)                                                                                                                     ttGATGACGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAgtgt         <  2:663089/140‑1 (MQ=255)                                                                                                                      tgatgaCGATCAATATCATCGTATACTCTTCGCTGCCCGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTa        <  1:245296/140‑1 (MQ=255)                                                                                                                          gaCGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGcc    <  1:516384/140‑1 (MQ=255)                                                                                                                            cGATCAATATCATCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTg  <  1:576291/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                GTTGTATAGCGTGATAGCGCCAGGGAAAGTGTCTGCGGACGCGCCGCTGGAGTTGGTTTCCCGTGTCAGCGATGTGACCGTGCAGGAGGCTGCCGCCATCGCATGGCATATGCCGTTTGATGACGATCAATATCACCGTTTACTCTCCGCTGCCGGGTTATCGAAAAGCTGGACGCGGACGATGCAAAAGCGCCGACTGAGCGGCAAGATCGAAGATTTTTCCCGCCGGTTGTGGGGAAAATAACCCCCGGAGTGTAGGCCTG  >  NC_000913/4103253‑4103515

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