Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,131,004 +G coding (1170/2433 nt) metL → Bifunctional aspartokinase/homoserine dehydrogenase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,131,0001.G100.0% 131.6 / NA 37Q389R (CAG→CGG) metLBifunctional aspartokinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/21);  total (16/21)

CGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCTGG  >  NC_000913/4130863‑4131131
                                                                                                                                          |                                                                                                                                   
cgcgcAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGTCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAg                                                                                                                                    <  1:744068/139‑1 (MQ=255)
cgcgcAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAg                                                                                                                                    >  1:1196466/1‑139 (MQ=255)
cgcgcAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAg                                                                                                                                    >  2:762500/1‑139 (MQ=255)
cgcgcAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAg                                                                                                                                    <  2:715048/139‑1 (MQ=255)
cgcgcAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAg                                                                                                                                    >  2:1354808/1‑139 (MQ=255)
cgcgcAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAg                                                                                                                                    >  2:1041733/1‑139 (MQ=255)
                 tGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCtgttgg                                                                                                                    >  1:927294/1‑139 (MQ=255)
                               cATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCa                                                                                                      <  1:1354808/139‑1 (MQ=255)
                                         gCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCAccc                                                                                            >  1:800972/1‑139 (MQ=255)
                                                gctgcAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAAccc                                                                                     <  2:1275568/139‑1 (MQ=255)
                                                  tgcAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGc                                                                                   <  1:1343412/139‑1 (MQ=255)
                                                          tGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGc                                                                           >  2:1243421/1‑139 (MQ=255)
                                                            cTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCa                                                                         <  1:1041733/139‑1 (MQ=255)
                                                            cTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCa                                                                         <  1:1243421/139‑1 (MQ=255)
                                                                    cAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGa                                                                                                                 <  1:1086380/91‑1 (MQ=255)
                                                                    cAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGa                                                                                                                 >  2:1086380/1‑91 (MQ=255)
                                                                    cAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCt                                                                 >  1:683545/1‑139 (MQ=255)
                                                                      gAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTgg                                                               <  2:1196466/139‑1 (MQ=255)
                                                                            gCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTg                                                                                                        >  1:1317503/1‑92 (MQ=255)
                                                                            gCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTg                                                                                                        <  2:1317503/92‑1 (MQ=255)
                                                                            gCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGgcagc                                                           <  1:916364/137‑1 (MQ=255)
                                                                            gCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGgcagc                                                           >  2:916364/1‑137 (MQ=255)
                                                                            gCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCaa                                                         >  2:570024/1‑139 (MQ=255)
                                                                            gCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCaa                                                         >  1:593931/1‑139 (MQ=255)
                                                                            gCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCaa                                                         >  1:737043/1‑139 (MQ=255)
                                                                                      cgcTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGt                                                                                          <  2:272453/96‑1 (MQ=255)
                                                                                      cgcTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGt                                                                                          >  1:272453/1‑96 (MQ=255)
                                                                                                ccTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGa                                     <  1:570024/139‑1 (MQ=255)
                                                                                                ccTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGa                                     <  2:737043/139‑1 (MQ=255)
                                                                                                          gCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTgg                           <  1:762500/139‑1 (MQ=255)
                                                                                                              gATTACCTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGt                       <  2:800972/139‑1 (MQ=255)
                                                                                                                    cTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACt                                                       >  1:515068/1‑101 (MQ=255)
                                                                                                                    cTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACt                                                       <  2:515068/101‑1 (MQ=255)
                                                                                                                    cTGGCGAACTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATg                 <  2:927294/139‑1 (MQ=255)
                                                                                                                            cTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGt                                        <  2:342869/108‑1 (MQ=255)
                                                                                                                            cTGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGt                                        >  1:342869/1‑108 (MQ=255)
                                                                                                                             tGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCa                                                                               >  2:683585/1‑68 (MQ=255)
                                                                                                                             tGCGCCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCa                                                                               <  1:683585/68‑1 (MQ=255)
                                                                                                                              gcgcCTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAg       >  1:379271/1‑139 (MQ=255)
                                                                                                                                 ccTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCt    <  2:683545/139‑1 (MQ=255)
                                                                                                                                 ccTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCt    >  2:516957/1‑139 (MQ=255)
                                                                                                                                 ccTGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCt    <  2:379271/139‑1 (MQ=255)
                                                                                                                                   tGCGTCAGGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCtgg  <  2:593931/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                   
CGCGCAGGTACGCCCGCTGGCGGTTGGCGTACATAACGATCGCCAGTTGCTGCAATTTTGCTACACCTCAGAAGTGGCCGACAGTGCGCTGAAAATCCTCGACGAAGCGGGATTACCTGGCGAACTGCGCCTGCGTCAGGGGCTGGCGCTGGTGGCGATGGTCGGTGCAGGCGTCACCCGTAACCCGCTGCATTGCCACCGCTTCTGGCAGCAACTGAAAGGCCAGCCGGTCGAATTTACCTGGCAGTCCGATGACGGCATCAGCCTGG  >  NC_000913/4130863‑4131131

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: