Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,072,718	A→G	V49A (GTC→GCC)	rutB ←	ureidoacrylate amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,072,718	0	A	G	100.0%	149.3 / NA	49	V49A (GTC→GCC)	rutB	ureidoacrylate amidohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (28/21);  total (28/21)

TGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTC  >  NC_000913/1072572‑1072861                                                                                                                                                   |                                                                                                                                                tGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAAAGGCCgg                                                                                                                                             >  2:735473/1‑151 (MQ=255)   gTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGc                                                                                                                                           <  1:384172/151‑1 (MQ=255)    ttttCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTTGCAATGGCCggg                                                                                                                                            >  1:167443/1‑149 (MQ=255)            gCGTTCGATTTATGAAAATTCGGTGCGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTg                                                                                                                                  <  2:453419/151‑1 (MQ=255)              gTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTgag                                                                                                                                <  1:175924/151‑1 (MQ=255)                tCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTgaga                                                                                                                               <  2:239307/150‑1 (MQ=255)                  gATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTGATGTTGGCAATGGCCCGGCGGGGGGTTGTGACAt                                                                                                                            >  2:102929/1‑151 (MQ=255)                       aTGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGTAAAGGGCGGGGGGGGGGGTGGGACATCaaa                                                                                                                        >  2:100889/1‑150 (MQ=255)                            aaTTCGGTGAGCAGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGGAAGTTGGGGATTGCCGGGGGGGGTGTTGAGATATCAAAacc                                                                                                                     >  1:587667/1‑148 (MQ=255)                             aTTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGTATGTTTGGAATGGCCGGGCGGGGGGTTGTGGCATCAAACCCGGc                                                                                                                  >  2:23024/1‑150 (MQ=255)                             aTTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAAAGGCCGGGCGGGGGGTTGAGACATCAAACCCcgcg                                                                                                                 >  2:890421/1‑151 (MQ=255)                                cggggggCCGGGTCCGCCAGCCTCGACATACTGTTCACCCCATCCATTTTGAAACAAGAGGATCAGCAGCCCTGCCGCGCGCGCTGCGGTCACTGCGGTTTGAATGTTGACCATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCgag                                                                                                               <  2:127830/144‑1 (MQ=255)                                cGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTTGCGATGGCCCGGCGGGGGGTTTGGACATCCAAACCGGCcag                                                                                                               >  1:553556/1‑150 (MQ=255)                                          ggTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTTGCCATGGCCGGGGGGGTGGGTGAGACATCAAAACCGGCGGGATCTAAGTa                                                                                                      >  1:509362/1‑149 (MQ=255)                                          ggTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAgc                                                                                                    <  2:46477/151‑1 (MQ=255)                                                 cAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGGAATGGGCGGGGGAGGGGGTGAGAAATCAAAACCCGCGAGGTCTAAGTAGCCGCCTgg                                                                                             >  1:733770/1‑151 (MQ=255)                                                     cTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAAAGTTTGCCATGGCCGGGCGGGGGGGTTGGAACTCAAACCCGGCGAGATTTTAATAgccgccccgcg                                                                                           >  2:453809/1‑144 (MQ=255)                                                       cGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGCAAAGGGCGGGCGGGTGGTTGAGACATCAAACCCGGCGAGAACTAAGTAGGCGCCTGGCGTgg                                                                                        >  2:499120/1‑150 (MQ=255)                                                         aCATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTGATGTTGGGAATGGGCGGGCGGGGGGGTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCCGGCGGGGCa                                                                                      >  2:229723/1‑150 (MQ=255)                                                          cATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATa                                                                                    <  2:528334/151‑1 (MQ=255)                                                                 ttCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCg                                                                                <  2:456170/148‑1 (MQ=255)                                                                     tCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGtttt                                                                            <  2:482330/148‑1 (MQ=255)                                                                      cccAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCa                                                                         <  2:635987/150‑1 (MQ=255)                                                                          gCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCCATGGCCGGGCGGGGGGGTGGGACATCAAACCCGGCGAGGTCTAAGTAGCCGCCCGGCGTGGCATAAGCGGTTTTCATATcc                                                                    >  2:18257/1‑151 (MQ=255)                                                                             aTTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTTGGAATGGCCGGGGGGGGGGTTGGGAAATCAAACCCGGCGCGATCTTAATAGCCGCCCGGCGGGGGACAAGAGGTTTTCATATccacc                                                                 >  2:803349/1‑151 (MQ=255)                                                                               tttGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTGATGGTGGCAACAGCCGGGCCGGCGGTTGGGACATCAAAACCCGCGGGGTCTAAAGAGCCCCCCGGCGGGGGCTAAACGGTTTTCATATccacc                                                                 >  2:819070/1‑149 (MQ=255)                                                                               tttGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGTGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATccacca                                                                <  1:503169/150‑1 (MQ=255)                                                                                 tGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGCATGTTGGGCATGGCCGGGGGAGTGGGTGAGACATCAAAACCGGCGGGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGGTTTTCATATccacccc                                                               >  1:881594/1‑147 (MQ=255)                                                                                    aaCCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCa                                                          <  2:675027/151‑1 (MQ=255)                                                                                         gatgatCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGGTTGGACGTTGGCCATGGCCGGGCGGGGGGGGGGGGAATTAAAACCCGGGAGGTCTAAGCAGCCCCCCGGCGGGGGGTAAGCGGTTTGCATAaccacccccagcaacgc                                                      >  2:70168/1‑139 (MQ=255)                                                                                             atCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTTAATGTTGGCAAAGGCCGGGCGAGTGGGTGAGACATCAAAACCGGCGAGATCTAAATAGCCGCCCGGCGTGGCATAAGCGGTTTGCATAACCACCACGATCAGCGCACtt                                                  >  2:716435/1‑150 (MQ=255)                                                                                              tCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACtt                                                  <  2:726175/149‑1 (MQ=255)                                                                                               cAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGGTTTAATGTTGGCAATGGGCGGGCGAGTGGGTGGGAGATCAAACCCCGCCGGCTCTAAGAAACCCGCTGGCGTGGGATAAGCGTTTTGCATATCCAACCCGAACAACGCACTTtg                                                >  1:155225/1‑150 (MQ=255)                                                                                                aGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTtgct                                              <  1:454921/151‑1 (MQ=255)                                                                                                  ctccccTGCCGCTCGCGCTCCTGTCACGGCGGTTTGATTGTTGGCAATGGCCGGGCTATTGGTTGAGACATCAATCCCGGCGAGTTCTAAGTAGCCGCCTGTCGTTGCATAAGCGTTTTGCATATCCCCCGCGAGCAGCGTACTTtgctg                                             <  2:549086/147‑1 (MQ=255)                                                                                                  cATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTtgctgc                                            <  1:179232/151‑1 (MQ=255)                                                                                                   aTCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTTGCAATGGCCGGGCGGGTGGGTGAGACATCAAACCCCGCGAGATCTAAGTAGCCGCCTGGCGGGGCGTAAGCGGTTTGCATATCCACCCCGATCAGCGCACttt                                                 >  2:805470/1‑145 (MQ=255)                                                                                                   aTCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAAAGGCCGGGCGAGGGGGTGAGACATCAAAACCCGCGAGATCTAAAGAGCCGCCTGGCGTGGCATAAGCGTTTTTCATATCCACCACGATCAGCGCACTTttctg                                             >  1:396924/1‑149 (MQ=255)                                                                                                    tCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCg                                           <  2:744701/150‑1 (MQ=255)                                                                                                          ccGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGGGTGGGTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGGATAAGCGtt                                                                              >  1:697261/1‑109 (MQ=255)                                                                                                          ccGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGtt                                                                              <  2:697212/109‑1 (MQ=255)                                                                                                                    tGCGGTCACGGCGGGTTGAATGTTGGGAATGGCCGGGGGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCCGGCGGGGCATAAGCGGTTTGCATATCCACCCCGATCAGCGCACTTTGCTGCGGATCGGAGGGAATgg                            >  2:649610/1‑149 (MQ=255)                                                                                                                            cggcggTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAg                                   <  1:644505/134‑1 (MQ=255)                                                                                                                            cggcggTTTGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAg                                   >  2:644456/1‑134 (MQ=255)                                                                                                                                  tttGAATGTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGt            >  2:180811/1‑151 (MQ=255)                                                                                                                                         gTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGgtc     >  2:753141/1‑151 (MQ=255)                                                                                                                                         gTTGGCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGgtc     <  1:691133/151‑1 (MQ=255)                                                                                                                                            ggCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGgtcgtc  >  1:498182/1‑151 (MQ=255)                                                                                                                                            ggCAATGGCCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGgtcgt   <  1:836061/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                TGGTTTTCAGGGCGTTCGATTTATGAAAATTCGGTGAGCCGGGTCCGCCAGCCTCGACATACTGTTCATCCCAGCCATTTTGAAACCAGATGATCAGCATCCCTGCCGCTCGCGCTGCGGTCACGGCGGTTTGAATGTTGGCAATGACCGGGCGAGTGGTTGAGACATCAAACCCGGCGAGATCTAAGTAGCCGCCTGGCGTGGCATAAGCGTTTTGCATATCCACCACGATCAGCGCACTTTGCTGCGGATCGAAGGTAATGGCTTCCGGTCGAGCGGTTAAGGTCGTC  >  NC_000913/1072572‑1072861

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