Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913_3_hsa_pgi	4,235,433	Δ1 bp	100%	intergenic (‑2200/‑501)	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,235,433	0	A	.	100.0%	166.6 / NA	38	intergenic (‑2200/‑501)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (23/15);  total (23/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCAC  >  NC_000913_3_hsa_pgi/4235291‑4235573                                                                                                                                               |                                                                                                                                             cAGTGGGGCGTGGAGGTGGTAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCGCAAGTGAACTCTCACGACGCTTCTAACAATGGGCTCATCAACTTCATCAAGTCGGAGCGCGAGTCCGGAGGCCAAT‑Aatcatc                                                                                                                                       >  2:23713/1‑149 (MQ=255) cAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑Aatcatc                                                                                                                                       >  1:77653/1‑149 (MQ=255) cAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAACGCGAGGCCAGAGTCCAAT‑Aatcatc                                                                                                                                       >  1:76090/1‑149 (MQ=255) cAGTGGGGAGTAGAGCTGGCAAAGCAGCCGGCTAAGAAAACACATCCGGAGATTCATGGCTATGCTGAAGTGACCTCTCACCACGCTTCTACCAATCGGCTCATCAACTTCATCAAGCAGCACCGCGGGCCCAGAGTCAAATaaccatc                                                                                                                                        >  2:307108/1‑144 (MQ=255)   gTGGGGCGTGGAGCTGGGAAAGCAGCTGGCTAAGACAATAGAGCCTGAGCTCGCTGGCAGTGCCCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCAtcgt                                                                                                                                     <  2:177874/149‑1 (MQ=255)    tGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCTAGGCCAGAGTCCAAT‑AATCAtcgt                                                                                                                                     >  2:40097/1‑148 (MQ=255)    tGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCAtcgt                                                                                                                                     <  1:40097/148‑1 (MQ=255)         aGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTAATGGCAGTGCTCAAGTGACCGCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCTAGAGTCCCAT‑AATCATCGTCGatat                                                                                                                               >  1:600677/1‑149 (MQ=255)            ggAGCTTGGAAGGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGAGGGCAGTGCTCAAGTGACCTCGCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTa                                                                                                                            <  2:584948/149‑1 (MQ=255)             gAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTAAAGCCGCTTCTAGCAATGGTCTCATCAACTTCATCCAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATATAg                                                                                                                           >  2:405420/1‑149 (MQ=255)                      aaGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTCCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATAGTAGATATGTAGGCCGGCTaa                                                                                                                  >  2:460792/1‑149 (MQ=255)                         cagcTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGc                                                                                                               <  1:2737/149‑1 (MQ=255)                               gCTAAGAAAATAGAGCCTGAGCTTGATGGTAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCAGCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                         <  2:92794/149‑1 (MQ=255)                               gCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                         >  1:36198/1‑149 (MQ=255)                                 tAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTATACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGCGGCCAGAGTCCAAT‑AATCATCGCCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                       >  2:568025/1‑149 (MQ=255)                                                 gggCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAACTGGCCCATCAACTCCATCAAGCAGCATCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGa                                                                                       <  1:357474/147‑1 (MQ=38)                                                    cTTGATGGCAGTGATCAGGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACTCCGCATCCGGCAACCGATGc                                                                                    <  1:647894/149‑1 (MQ=255)                                                              gTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAc                                                                          <  2:406923/149‑1 (MQ=255)                                                                 cTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgcg                                                                       >  1:534948/1‑149 (MQ=255)                                                                  tCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                                                                      >  1:45035/1‑149 (MQ=255)                                                                                 cgacgCTTCTACCAATGGGCTCATCAACTTCAGCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAgg                                                       <  2:468174/149‑1 (MQ=255)                                                                                      cTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTaa                                                  >  1:394338/1‑148 (MQ=255)                                                                                          tACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                              >  1:285125/1‑149 (MQ=255)                                                                                             cAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCAAAT‑AATAATCGTAGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGACGCATGATGAGACGCGGTCGCGTCTGTTCAGGACTAAAGTTcca                                           >  2:94741/1‑147 (MQ=25)                                                                                             cAAAGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                           <  2:119547/149‑1 (MQ=255)                                                                                              aaTGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCGTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCCGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCACGCCTACAGGTCGAt                                          >  2:205512/1‑149 (MQ=255)                                                                                                tGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATGCGGCAACCGAGGCTTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                        <  1:157065/149‑1 (MQ=255)                                                                                                   gCTCATCAACTTCCTCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGa                                     <  1:299199/149‑1 (MQ=255)                                                                                                      catcgACTTCATCAAGCAGGAGCCCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGCCCTTCCCCCCGCATCCACCAACCTTTGTCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatat                                  <  1:562862/149‑1 (MQ=25)                                                                                                       atcaACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGataag                                 >  2:389876/1‑147 (MQ=255)                                                                                                           aCTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACa                             <  2:534948/149‑1 (MQ=255)                                                                                                            cTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGGCGGTATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCACCGATGCCTGATGCGACGCGGTCGCGTCTTATCGGGCCTACCGGTCGATGCCAATCTGTACAt                            >  2:205336/1‑149 (MQ=25)                                                                                                             ttCATCAAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATc                           >  1:352211/1‑149 (MQ=255)                                                                                                                atcaAGCAGCAGCGCGAGGTCAAAGTCCAAT‑AATCATCTTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTa                        <  2:448189/149‑1 (MQ=255)                                                                                                                atcaAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTCTCACGCCTACAGGTCGATGCCGCTATGTAAATCGTa                        >  2:633197/1‑149 (MQ=37)                                                                                                                atcaAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACGTCGTa                        >  1:377604/1‑149 (MQ=255)                                                                                                                atcaAGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGACAGGTACATCGTa                        >  1:103325/1‑149 (MQ=255)                                                                                                                   aaGCAGCAGCGCGAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATACGTACATCGTATTc                     >  1:500470/1‑149 (MQ=255)                                                                                                                               gAGGCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATAGTTACAACGTATTCGGCAATtactac         >  2:476886/1‑149 (MQ=25)                                                                                                                                  gCCAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATa      >  1:343429/1‑149 (MQ=255)                                                                                                                                   ccAGAGTCCAAT‑AATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg     >  1:380086/1‑149 (MQ=255)                                                                                                                                      gagTCCAAT‑AATCAACGTCGATATCTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCAc  >  1:576220/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                             CAGTGGGGAGTGGAGCTGGGAAAGCAGCTGGCTAAGAAAATAGAGCCTGAGCTTGATGGCAGTGCTCAAGTGACCTCTCACGACGCTTCTACCAATGGGCTCATCAACTTCATCAAGCAGCAGCGCGAGGCCAGAGTCCAATAAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCAC  >  NC_000913_3_hsa_pgi/4235291‑4235573

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