Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,194,999	A→G	H367H (CAT→CAC)	thiC ←	phosphomethylpyrimidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,194,999	0	A	G	100.0%	105.8 / NA	32	H367H (CAT→CAC)	thiC	phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (19/13);  total (19/13)

TCGGCAAACTGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCT  >  NC_000913/4194865‑4195125                                                                                                                                        |                                                                                                                               tCGGCAAACTGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGGAATTTTCCT‑Ggtgatg                                                                                                                           >  2:497450/1‑140 (MQ=255) tCGGCAAACTGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑Ggtgatg                                                                                                                           <  2:560492/140‑1 (MQ=255) tCGGCAAACTGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑Ggtgatg                                                                                                                           >  2:567379/1‑140 (MQ=255)   ggCAAACTGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATggg                                                                                                                         <  2:21427/140‑1 (MQ=255)    gCAAACTGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGga                                                                                                                        <  1:174842/140‑1 (MQ=255)    gCAAACTGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGga                                                                                                                        <  1:470296/140‑1 (MQ=255)    gCAAACTGCGCTTCATCATTGGCGTCCTGAAAAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGga                                                                                                                        <  2:267818/140‑1 (MQ=255)          tGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCa                                                                                                                  <  1:491880/140‑1 (MQ=255)          tGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACACCGAAACGTCATAAGCGGCACAGATTTCACAAAATTCGCGGAAGTGGTGATAGAGGAAACTTTCGT‑GGTGATGGGAGAGGCa                                                                                                                  >  2:299568/1‑140 (MQ=255)           gcgcTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCAc                                                                                                                 >  1:117584/1‑140 (MQ=255)                  aTCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCg                                                                                                          >  1:63816/1‑140 (MQ=255)                      ttGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCAt                                                                                                      >  2:831196/1‑140 (MQ=255)                              cTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAAc                                                                                              >  1:506416/1‑140 (MQ=255)                                     gAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCcgcg                                                                                       <  1:651694/140‑1 (MQ=255)                                     gAACCGGGGCGCAGACCGTCGCACAGCGAAAGCGAAATGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCcgcg                                                                                       <  1:303638/140‑1 (MQ=255)                                     aacccGGGGCGCAGTCCGTCGCCAAGGGACAGGAAAAAGCCAAAAGCGGCCCGGATTTCAAAAATTTCGCGGAAGTGTTGAAAGGGGAAATTTTCCT‑GGTGAGGGGTGGGGAACCATTTCGCCATAAGCGAACCGCcgcg                                                                                       <  2:484293/137‑1 (MQ=255)                                        ccGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCgaga                                                                                    >  2:571097/1‑140 (MQ=255)                                          ggggCGCAGCCCGTCGCTCAGCTACAGCGAAACGTCATACGCGGCACACATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTCCCG‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCgagaga                                                                                  <  2:217078/140‑1 (MQ=255)                                           gggCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGAc                                                                                 >  1:441952/1‑140 (MQ=255)                                           gggCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGAc                                                                                 >  2:621826/1‑140 (MQ=255)                                                cAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGGCATGATGGAACCGGCGCGAGAGACGATAc                                                                            >  1:217078/1‑140 (MQ=255)                                                cAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATAc                                                                            >  1:70599/1‑140 (MQ=255)                                                   aCCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCgg                                                                         >  1:603563/1‑140 (MQ=255)                                                             agcgacagcgaAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGttt                                                               >  1:125956/1‑140 (MQ=255)                                                                      aaaaCGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGAAAATTTTCCT‑GGTGAGGGGAGAGGCACCATTCCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCAt                                                      <  2:641091/139‑1 (MQ=255)                                                                       aaaCGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATc                                                     >  2:710319/1‑140 (MQ=255)                                                                        aaCGTCATAAGGGGCACAGATTTCAAAAATTTCGCGGAAGTGTTAATAGAGAAAATTTCCCT‑GGTGAGGGGAGAGCCCCCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCg                                                    <  2:125956/140‑1 (MQ=255)                                                                           gTCATAAGCGGCCAAAATTTCAAAAATTTCGGGGAGGTGTTGATGGAGGAATTTTTCCG‑GGTGTGGGGAGAGGCCCCTTTTCCCCAAAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGca                                                 <  2:770514/140‑1 (MQ=255)                                                                           gTCAAAAGCGGCACAGTTTCCAAAATTTTCGCGGAAGTGTTAATAGAGGAAATTTTCCG‑GGTGAGGGGAGGGGCCCCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGca                                                 <  2:70599/140‑1 (MQ=255)                                                                             aaTAAGCGGCACAGTTTTCACAATTTTCGCGGAAGTGTTGATAGAGGAATTTTTCCG‑GGTGATGGGAGAGGCCCCTTTCCGCAATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGcaca                                               <  2:168482/139‑1 (MQ=255)                                                                               tAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATa                                             >  2:466067/1‑140 (MQ=255)                                                                               tAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATa                                             >  2:663621/1‑140 (MQ=255)                                                                                   cGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGAGGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAgcgc                                         <  2:117584/140‑1 (MQ=255)                                                                                      cacaGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCgcagc                                      >  1:144718/1‑140 (MQ=255)                                                                                      cacaGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCgcagc                                      >  1:168482/1‑140 (MQ=255)                                                                                            tttCACAAATTTCGGGGAAGTGTTGATAGAGGAAATTTTCCG‑GGTGATGGGAGAGGCACCTTTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACg                                <  1:66788/140‑1 (MQ=255)                                                                                                 caAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGc                           <  1:663621/140‑1 (MQ=255)                                                                                                 aaaaTTCCGGGAAGGTGTTGAAAGGGGAAATTTCCCT‑GGTGAGGGGAGAGGCACCTTTCCCCAATAATCGAGCCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGc                           <  2:506416/139‑1 (MQ=255)                                                                                                         gcgGAAGTGTTGATAGAGGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAg                   <  1:710319/140‑1 (MQ=255)                                                                                                                   tGATAG‑GGAAATTTTCCTGGGGGATGGGAGAGGACCGATTCCGCTATAAAAGATCCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATc         <  2:603563/140‑1 (MQ=255)                                                                                                                         agGAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACAcc   >  1:445798/1‑140 (MQ=255)                                                                                                                          ggAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCt  >  1:633273/1‑140 (MQ=255)                                                                                                                          ggAAATTTTCCT‑GGTGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCt  >  2:49289/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                               TCGGCAAACTGCGCTTCATCGTTGGCGTCCTGAATAGAACCGGGGCGCAGACCGTCGCCCAGCGACAGCGAAACGTCATAAGCGGCACAGATTTCACAAATTTCGCGGAAGTGTTGATAGAGGAAATTTTCCT‑GATGATGGGAGAGGCACCATTTCGCCATAATCGAACCGCCGCGAGAGACGATACCGGTCAGGCGTTTCGCGGTCATCGGCACATAGCGCAGCAGTACGCCCGCATGGATAGTGAAGTAATCCACACCT  >  NC_000913/4194865‑4195125

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: ‑