Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_hsa_pgi	4,233,746	A→G	intergenic (‑513/‑2188)	lysC ← / → yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_hsa_pgi	4,233,746	0	A	G	100.0%	185.2 / NA	55	intergenic (‑513/‑2188)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/24);  total (31/24)

CCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAACCACTTCAGCTTGACCC  >  NC_000913_3_hsa_pgi/4233600‑4233887                                                                                                                                                   |                                                                                                                                              ccACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagagga                                                                                                                                             >  2:1059708/1‑149 (MQ=255)  cACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAgaggag                                                                                                                                            >  1:595143/1‑149 (MQ=255)       cTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGAGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTg                                                                                                                                       >  2:480132/1‑149 (MQ=255)             gTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATg                                                                                                                                 >  2:485841/1‑149 (MQ=255)                 gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCg                                                                                                                             <  1:1059708/149‑1 (MQ=255)                 ctagtTCCGCTAATAAATGCTTCACTACGCTAAGTGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTATGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGCGCAGAATTGCTAAtcgtcg                                                                                                                             <  2:645559/148‑6 (MQ=255)                  tagtTCAGCTAAGAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGc                                                                                                                            <  2:364431/149‑1 (MQ=255)                     ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGctct                                                                                                                         <  2:544182/149‑1 (MQ=255)                     ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTATCAAAATCAGAGGAGTATTTCTAATGGCCGctct                                                                                                                         >  1:407955/1‑149 (MQ=255)                      tCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCCAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCtctc                                                                                                                        <  1:316939/149‑1 (MQ=255)                        aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCtctccc                                                                                                                      >  1:541780/1‑147 (MQ=255)                         gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCAcc                                                                                                                     <  1:1100629/149‑1 (MQ=255)                         gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTACTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATGGCTAATGGCCGCACACAcc                                                                                                                     >  2:658031/1‑149 (MQ=255)                           taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGTCCACCCTCACCCg                                                                                                                   >  2:543525/1‑149 (MQ=255)                             ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCggg                                                                                                                 >  2:1016221/1‑149 (MQ=255)                             ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCCCTCACCCggg                                                                                                                 >  2:1016312/1‑149 (MQ=255)                               aaaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGAc                                                                                                               >  2:386507/1‑149 (MQ=255)                                 aTGCTTCCCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGAccc                                                                                                             <  2:1045106/149‑1 (MQ=255)                                   gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGAccccc                                                                                                           >  1:144509/1‑149 (MQ=255)                                     ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAACGTCACAATTCTCAAAATCAGAGGCGTATTGCTAATGGCCGCTCTCCCCCGGCCCCCCGAg                                                                                                         >  1:1060587/1‑149 (MQ=255)                                          tGCGCAAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCa                                                                                                    <  1:543525/149‑1 (MQ=255)                                             gcTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCATAATTCTCAAAATCAGAGGAGTAGTGCTAATGGCCGCGCTCACCCGGGACCCCCAGTTCCAGaa                                                                                                 <  1:173648/149‑1 (MQ=255)                                                   ggTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATCTCTGCCTTGCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTgca                                                                                           <  1:32713/149‑1 (MQ=255)                                                          acTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGCAGTATTGCTAATGGCCGCTCTCACCCGGCACCCCCAGTTCCAGAAACTGCAGCAATgg                                                                                    >  2:546989/1‑149 (MQ=255)                                                            tCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTAACAATTCTCAAAATAAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCGATGGTa                                                                                  >  1:411005/1‑149 (MQ=255)                                                                  ttACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGa                                                                            >  2:313065/1‑149 (MQ=255)                                                                  ttACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGa                                                                            >  1:196052/1‑149 (MQ=255)                                                                   tACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAg                                                                           <  1:250087/149‑1 (MQ=255)                                                                               aCTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAg                                                               >  2:1102350/1‑149 (MQ=255)                                                                                 tAAAACCATCACGTTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCt                                                             <  2:972101/149‑1 (MQ=255)                                                                                 tAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCt                                                             >  1:464128/1‑149 (MQ=255)                                                                                 tAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCt                                                             <  1:485841/149‑1 (MQ=255)                                                                                      ccATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTgcag                                                                                          <  2:1031508/115‑1 (MQ=255)                                                                                      ccATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTgcag                                                                                          >  1:1031508/1‑115 (MQ=255)                                                                                      ccATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAAcc                                                        >  1:54392/1‑149 (MQ=255)                                                                                         tCACATTTTTTTTTTTCTTTCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTgc                                                     >  2:129052/1‑149 (MQ=255)                                                                                         tCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGAcccc                                                                                                            <  1:351699/94‑1 (MQ=255)                                                                                         tCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGAcccc                                                                                                            >  2:351699/1‑94 (MQ=255)                                                                                         tCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCt                                                                    <  1:21731/134‑1 (MQ=255)                                                                                         tCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCt                                                                    >  2:21731/1‑134 (MQ=255)                                                                                         tCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTgc                                                     >  1:440667/1‑149 (MQ=255)                                                                                                       gACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGAt                                       <  1:480132/149‑1 (MQ=255)                                                                                                        aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCGCCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATg                                      <  2:143709/149‑1 (MQ=255)                                                                                                        aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTTAACCTGCGCCGCCTCTTCGATg                                      <  2:483203/149‑1 (MQ=255)                                                                                                        aCTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATg                                      >  1:94798/1‑149 (MQ=255)                                                                                                            gcgcTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGCACCTCCGCCGCCCCTCCGATGcaac                                  >  2:379747/1‑146 (MQ=255)                                                                                                             cgcTACACTCTGTCAAACTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCcaac                                 <  2:26088/149‑1 (MQ=255)                                                                                                                tACAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAg                              <  2:196052/149‑1 (MQ=255)                                                                                                                  cAATCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGa                            >  1:716912/1‑149 (MQ=255)                                                                                                                   aaTCTTCCAAAGTCACAATTCTCAAAATCAGAGGAGTATTCCTCACGGCCACCCTCACACGGGACCCCCATTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGGGCTGAACCTGCGCCGCCTCTTCGCTGCCCAGCAGGAc                           >  1:268402/1‑149 (MQ=255)                                                                                                                            aaGTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAAc                  <  2:94798/149‑1 (MQ=255)                                                                                                                              gTCACAATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAACCa                <  1:1102350/149‑1 (MQ=255)                                                                                                                                cacaATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAACCACt              >  1:75116/1‑149 (MQ=255)                                                                                                                                cacaATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAACCACt              >  1:227545/1‑149 (MQ=255)                                                                                                                                 acaATTCTCAAAATCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAACCACtt             >  1:191220/1‑149 (MQ=255)                                                                                                                                           aaaTCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAACCACTTCAGCTTGAcc   <  2:144509/149‑1 (MQ=255)                                                                                                                                            aaTCAGAGGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAACCACTTCAGCTTGAccc  >  2:673776/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                              CCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGGCCGCTCTCACCCGGGACCCCCAGTTCCAGAAGCTGCAGCAATGGTACCGCGAGCACCGCTCCGAGCTGAACCTGCGCCGCCTCTTCGATGCCAACAAGGACCGCTTCAACCACTTCAGCTTGACCC  >  NC_000913_3_hsa_pgi/4233600‑4233887

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