Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,102,930	C→T	L37L (CTG→TTG)	yiiM →	6‑N‑hydroxylaminopurine resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,102,930	0	C	T	100.0%	125.1 / NA	32	L37L (CTG→TTG)	yiiM	6‑N‑hydroxylaminopurine resistance protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/15);  total (17/15)

AATGGGGAAATTTCTCGTGGAGAGGGAACAGATGCGATATCCGGTTGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGCTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGAGTTGTTTGTTGCGCCTGCGTTTGGTG  >  NC_000913/4102791‑4103077                                                                                                                                            |                                                                                                                                                    aaTGGGGAAATTTCTCGTGGAGAGGGAACAGATGCGATATCCGGTTGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGa                                                                                                                                           >  2:103022/1‑150 (MQ=255)       gAAATTTCTCGTGGAGAGGGAACAGATGCGATATCCGGTTGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGAc                                                                                                                                    <  2:70949/151‑1 (MQ=255)                cGTGGAGAGGGAACAGATGCGATATCCGGTTGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGc                                                                                                                            <  2:388006/150‑1 (MQ=255)                   ggAGAGGGAACAGATGCGATATCCGGTTGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCgg                                                                                                                          >  1:10883/1‑149 (MQ=255)                                 tGCGATATCCGGTTGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCAcgg                                                                                                          >  2:32689/1‑151 (MQ=255)                                          cGGTTGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGcc                                                                                                    >  2:458701/1‑148 (MQ=255)                                              tGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCcagacag                                                                                              <  2:569457/150‑1 (MQ=255)                                               gATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCcagacaga                                                                                             <  1:228569/150‑1 (MQ=255)                                                                    cAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTg                                                                       <  1:53278/151‑1 (MQ=255)                                                                          tATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCa                                                                   <  1:244384/149‑1 (MQ=255)                                                                             cccGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTAt                                                               <  1:298526/150‑1 (MQ=255)                                                                             cccGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTAt                                                               >  2:459676/1‑150 (MQ=255)                                                                              ccGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTAtct                                                             <  2:492557/151‑1 (MQ=255)                                                                              ccGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTAtct                                                             >  2:350472/1‑151 (MQ=255)                                                                                 aaGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTAtctc                                                            >  1:6665/1‑149 (MQ=255)                                                                                   ggCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACt                                                        >  2:380739/1‑151 (MQ=255)                                                                                    gCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTg                                                       >  1:619841/1‑151 (MQ=255)                                                                                       aaCCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTggg                                                     >  2:100167/1‑150 (MQ=255)                                                                                               gcAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAAt                                            >  1:157332/1‑151 (MQ=255)                                                                                               gcAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAAt                                            <  1:578610/151‑1 (MQ=255)                                                                                                cAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAAtt                                           >  2:552165/1‑151 (MQ=255)                                                                                                 aaTTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAAttt                                          <  2:418622/151‑1 (MQ=255)                                                                                                    tGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAAttttc                                        >  1:635613/1‑148 (MQ=255)                                                                                                        aaaaTCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGaa                                    >  1:646599/1‑150 (MQ=255)                                                                                                        aaaaTCCAGGTTGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCAACTGGGCGCGGGAATTTCCGGa                                     <  1:476298/149‑1 (MQ=255)                                                                                                                   tGATGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGAgt                          >  1:2633/1‑149 (MQ=255)                                                                                                                      tGGTGAGTTGATGCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGGGttgttt                     >  1:291746/1‑151 (MQ=255)                                                                                                                                  gCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGAGTTGTTTGTTGCGCCTGc          <  1:457771/150‑1 (MQ=255)                                                                                                                                  gCTGACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGAGTTGTTTGTTGCGCCTGc          <  1:409258/150‑1 (MQ=255)                                                                                                                                     gACGGAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGAGTTGTTTGTTGCGCCTGCGtt       <  2:618552/150‑1 (MQ=255)                                                                                                                                        ggAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGAGTTGTTTGTTGCGCCTGCGTTTgg    >  2:240225/1‑150 (MQ=255)                                                                                                                                         gAGTTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGAGTTGTTTGTTGCGCCTGCGTTTGGTg  <  2:75005/151‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                    AATGGGGAAATTTCTCGTGGAGAGGGAACAGATGCGATATCCGGTTGATGTATACACAGGCAAGATTCAGGCTTATCCCGAAGGCAAACCCAGCGCAATTGCTAAAATCCAGGTTGATGGTGAGTTGATGCTGACGGAGCTGGGGCTGGAAGGTGACGAGCAGGCGGAGAAAAAGGTTCACGGCGGGCCAGACAGAGCGCTGTGTCATTATCCTCGTGAGCATTATCTCTACTGGGCGCGGGAATTTCCGGAACAGGCGGAGTTGTTTGTTGCGCCTGCGTTTGGTG  >  NC_000913/4102791‑4103077

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