Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,395,058	C→T	G331G (GGC→GGT)	nnr →	bifunctional NAD(P)H‑hydrate repair enzyme, C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,395,058	0	C	T	100.0%	148.6 / NA	46	G331G (GGC→GGT)	nnr	bifunctional NAD(P)H‑hydrate repair enzyme, C‑terminal domain ADP‑dependent (S)‑NAD(P)H‑hydrate dehydratase and N‑terminal domain NAD(P)H‑hydrate epimerase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/21);  total (25/21)

GCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGATCACG  >  NC_000913/4394934‑4395184                                                                                                                             |                                                                                                                               gCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTgggg                                                                                                                  >  2:917264/1‑139 (MQ=255)  cAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTggggg                                                                                                                 >  2:1367608/1‑139 (MQ=255)  cAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTggggg                                                                                                                 >  2:1412318/1‑139 (MQ=255)  cAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTggggg                                                                                                                 >  2:164669/1‑139 (MQ=255)      gAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGaaaa                                                                                                             <  2:517834/139‑1 (MQ=255)       aaaaCATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGaaaaa                                                                                                            >  2:821214/1‑139 (MQ=255)           cATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCa                                                                                                        <  1:1367608/139‑1 (MQ=255)               gcgcCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGc                                                                                                    >  2:569391/1‑139 (MQ=255)                    gctgctGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCaaaaa                                                                                               >  2:1461325/1‑139 (MQ=255)                          gACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTgag                                                                                         <  1:917264/139‑1 (MQ=255)                               cACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAAttt                                                                                    >  2:680651/1‑139 (MQ=255)                                       gAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAAc                                                                            <  1:569391/139‑1 (MQ=255)                                       gAATTGATGCTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAAc                                                                            <  1:1108189/139‑1 (MQ=255)                                               ggTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGAtgttg                                                                    >  1:1547349/1‑139 (MQ=255)                                                 tGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTtgtg                                                                  <  1:1540233/139‑1 (MQ=255)                                                  gCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTgg                                                                 >  2:1209133/1‑139 (MQ=255)                                                    aTGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGc                                                                                                         <  1:184295/97‑1 (MQ=255)                                                    aTGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGc                                                                                                         >  2:184295/1‑97 (MQ=255)                                                         cTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCg                                                          <  1:1461325/139‑1 (MQ=255)                                                           gACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCa                                                                                                        <  2:542476/91‑1 (MQ=255)                                                           gACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCa                                                                                                        >  1:542476/1‑91 (MQ=255)                                                           gACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGAt                                                        >  1:891406/1‑139 (MQ=255)                                                            aCGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATg                                                       <  1:696852/139‑1 (MQ=255)                                                                  gACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCa                                                                                                        >  2:601401/1‑84 (MQ=255)                                                                  gACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCa                                                                                                        <  1:601401/84‑1 (MQ=255)                                                                  gACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGa                                                 <  1:1209133/139‑1 (MQ=255)                                                                       tctTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGa                                                         <  1:934368/126‑1 (MQ=255)                                                                       tctTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGa                                                         >  2:934368/1‑126 (MQ=255)                                                                                 aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTgg                                                                                                                    >  1:511391/1‑57 (MQ=255)                                                                                 aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTgg                                                                                                                    <  2:511391/57‑1 (MQ=255)                                                                                 aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGc                                                                                                         >  2:630400/1‑68 (MQ=255)                                                                                 aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGc                                                                                                         <  1:630400/68‑1 (MQ=255)                                                                                 aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGAt                                                        <  1:848434/117‑1 (MQ=255)                                                                                 aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGAt                                                        >  2:848434/1‑117 (MQ=255)                                                                                 aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaatca                                  >  2:1416550/1‑139 (MQ=255)                                                                                   ccTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGAtgttg                                                                    <  1:1140778/103‑1 (MQ=255)                                                                                   ccTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGAtgttg                                                                    >  2:1140778/1‑103 (MQ=255)                                                                                          tGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATa                         <  1:1252392/139‑1 (MQ=255)                                                                                            ggCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCg                                                          <  1:954252/104‑1 (MQ=255)                                                                                            ggCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCg                                                          >  2:954252/1‑104 (MQ=255)                                                                                                      gtggtgATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATc             >  1:147797/1‑139 (MQ=255)                                                                                                              tGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGCCAATCAATCCCGATa                         >  2:1402758/1‑119 (MQ=255)                                                                                                              tGGTCCCGGTCTGGGTCAGCAAAAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATa                         <  1:1402758/119‑1 (MQ=255)                                                                                                                gTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGATCAc   >  1:311551/1‑139 (MQ=255)                                                                                                                 tCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGTCAATCAATCCCGATAAGCGTCACAATCGCGTGATCAcg  >  2:297970/1‑139 (MQ=255)                                                                                                                 tCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGATCAcg  <  2:891406/139‑1 (MQ=255)                                                                                                                             |                                                                                                                               GCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGATCACG  >  NC_000913/4394934‑4395184

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