Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,125,066	G→C	E729Q (GAG→CAG)	pycA →	pyruvate carboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,125,066	0	G	C	100.0%	126.4 / NA	34	E729Q (GAG→CAG)	pycA	pyruvate carboxylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (22/12);  total (22/12)

TCCTGAGCGATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTGAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCAATGAGTGGTTTAACAA  >  CP000730/1124927‑1125207                                                                                                                                            |                                                                                                                                              tCCTGAGCGATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTaaaa                                                                                                                                       >  2:342114/1‑148 (MQ=255)         gATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGt                                                                                                                               >  2:172918/1‑148 (MQ=255)                             gagTATTATGTCAAACTAGCTAANGAGTTAGAACGTNAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTc                                                                                                           >  1:399993/1‑148 (MQ=255)                             gagTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTc                                                                                                           >  2:274834/1‑148 (MQ=255)                             gagTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTc                                                                                                           >  1:352851/1‑148 (MQ=255)                                         aaaCTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATa                                                                                               >  1:21479/1‑148 (MQ=255)                                               gCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTg                                                                                         <  2:298190/148‑1 (MQ=255)                                               gCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTg                                                                                         <  2:160479/148‑1 (MQ=255)                                                   aaGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTaa                                                                                     >  1:253682/1‑148 (MQ=255)                                                   aaGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTaa                                                                                     >  1:253962/1‑148 (MQ=255)                                                      agTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATgg                                                                                  >  1:83440/1‑148 (MQ=255)                                                        ttAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGtt                                                                                <  1:353991/148‑1 (MQ=255)                                                        ttAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGtt                                                                                <  1:309120/148‑1 (MQ=255)                                                         tAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGttt                                                                               <  1:11118/148‑1 (MQ=255)                                                              cGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTaa                                                                          <  2:306165/148‑1 (MQ=255)                                                                                      aTTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTg                                                  <  1:182093/148‑1 (MQ=255)                                                                                              tatGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGat                                          <  1:76308/148‑1 (MQ=255)                                                                                                  gCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCa                                      >  2:374566/1‑148 (MQ=255)                                                                                                  gCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCa                                      >  2:305875/1‑148 (MQ=255)                                                                                                  gCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCa                                      >  2:214405/1‑148 (MQ=255)                                                                                                    aGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCAtt                                    >  1:385718/1‑148 (MQ=255)                                                                                                    aGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCAtt                                    >  1:241839/1‑148 (MQ=255)                                                                                                         tattaAAACCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATAc                               >  2:378641/1‑148 (MQ=255)                                                                                                               aaCCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAGTCAGCTGTAGATTTACCAATTCATCTGCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGt                         >  2:145312/1‑148 (MQ=255)                                                                                                               aaCCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCGCACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGACGCTGGTGTCGATCTCATTGAAActacggt                         >  2:163009/1‑143 (MQ=255)                                                                                                                aCCTAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGtt                        <  2:379576/148‑1 (MQ=255)                                                                                                                   tAAAGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCt                     >  1:68436/1‑148 (MQ=255)                                                                                                                     aaGCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTc                   >  1:190250/1‑148 (MQ=255)                                                                                                                       gCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCaa                 >  1:372096/1‑148 (MQ=255)                                                                                                                       gCCGCTTACGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCaa                 >  1:386872/1‑148 (MQ=255)                                                                                                                               cGAATTGATTGGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCAATGAGTGGt         <  2:253962/148‑1 (MQ=255)                                                                                                                                   ttgattgGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCAATGAGTGGTTTaa     >  1:147029/1‑148 (MQ=255)                                                                                                                                      attgGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCAATGAGTGGTTTAACaa  <  2:361166/148‑1 (MQ=255)                                                                                                                                      attgGTCAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCAATGAGTGGTTTAACaa  <  2:369394/148‑1 (MQ=255)                                                                                                                                            |                                                                                                                                              TCCTGAGCGATCAAACATTTATACTTTAGAGTATTATGTCAAACTAGCTAAAGAGTTAGAACGTGAAGGTTTCCATATTTTAGCGATTAAAGATATGGCAGGCTTATTAAAACCTAAAGCCGCTTACGAATTGATTGGTGAGTTAAAATCAGCTGTAGATTTACCAATTCATCTTCACACTCATGATACAAGTGGTAATGGTTTATTAACATACAAACAAGCAATAGATGCTGGTGTCGATATCATTGATACTGCTGTTGCTTCAATGAGTGGTTTAACAA  >  CP000730/1124927‑1125207

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