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%	120.6 / NA	38	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 (20/18);  total (20/18)

GTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGATCACG  >  NC_000913/4394924‑4395184                                                                                                                                       |                                                                                                                                gTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGc                                                                                                                             >  2:750311/1‑139 (MQ=255)      gACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAagag                                                                                                                        >  2:333716/1‑139 (MQ=255)          cGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑ggg                                                                                                                   >  1:608592/1‑139 (MQ=255)                          cgccgcTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCa                                                                                                   <  2:690373/139‑1 (MQ=255)                             cgcTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCaaaa                                                                                                <  2:995090/139‑1 (MQ=255)                               ctgctgACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAg                                                                                              >  1:322285/1‑139 (MQ=255)                                            gACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCg                                                                                 <  1:333716/139‑1 (MQ=255)                                            gACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCg                                                                                 >  2:987633/1‑139 (MQ=255)                                                 gAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAAc                                                                            <  1:263559/139‑1 (MQ=255)                                                        tGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCg                                                                          >  2:569738/1‑134 (MQ=255)                                                        tGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCg                                                                          <  1:569738/134‑1 (MQ=255)                                                           tGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTtgtg                                                                  <  1:987633/139‑1 (MQ=255)                                                            gCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTgg                                                                 >  1:760517/1‑139 (MQ=255)                                                                gAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATg                                                             <  2:760517/139‑1 (MQ=255)                                                                   cTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCg                                                          <  2:322285/139‑1 (MQ=255)                                                                    tGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGa                                                         <  2:511333/139‑1 (MQ=255)                                                                             aCTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAagag                                                                                                                        >  1:298091/1‑68 (MQ=255)                                                                             aCTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAagag                                                                                                                        <  2:298091/68‑1 (MQ=255)                                                                                         aaaGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCg                                                          >  2:190228/1‑117 (MQ=255)                                                                                         aaaGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCg                                                          <  1:190228/117‑1 (MQ=255)                                                                                         aaaGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaat                                    >  2:966316/1‑139 (MQ=255)                                                                                         aaaGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaat                                    >  1:562713/1‑139 (MQ=255)                                                                                           aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaatca                                  >  1:875534/1‑139 (MQ=255)                                                                                             ccTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaatcaat                                >  2:137405/1‑139 (MQ=255)                                                                                             ccTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaatcaat                                <  1:1046291/139‑1 (MQ=255)                                                                                             ccTGGAATGGGCCGATGTGGTGGTGAGTGGGCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaatcaat                                <  2:671228/139‑1 (MQ=255)                                                                                                    tGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATa                         <  1:137405/139‑1 (MQ=255)                                                                                                      ggCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAg                       >  2:1011977/1‑139 (MQ=255)                                                                                                       gCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAg                       >  2:546903/1‑139 (MQ=255)                                                                                                                      aTTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGa       >  2:308252/1‑139 (MQ=255)                                                                                                                      aTTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGa       >  2:761775/1‑139 (MQ=255)                                                                                                                      aTTGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGa       <  1:761775/139‑1 (MQ=255)                                                                                                                       ttGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGAt      <  1:966316/139‑1 (MQ=255)                                                                                                                        tGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCg                                                                                 <  1:859778/63‑1 (MQ=255)                                                                                                                        tGGTCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCg                                                                                 >  2:859778/1‑63 (MQ=255)                                                                                                                           tCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCa                                                     <  1:912845/88‑1 (MQ=255)                                                                                                                           tCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCa                                                     >  2:912845/1‑88 (MQ=255)                                                                                                                           tCCCGGTCTGGGTCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGATCAcg  >  1:796081/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                GTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGT‑GGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAGCGTCACAATCGCGTGATCACG  >  NC_000913/4394924‑4395184

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: ‑