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%	143.0 / NA	45	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 (22/23);  total (22/23)

AGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAG  >  NC_000913/4394923‑4395163                                                                                                                                        |                                                                                                          aGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAg                                                                                                        >  1:125126/1‑139 (MQ=255)        aCCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCaaaaat                                                                                                 >  1:644138/1‑135 (MQ=255)          ccGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTgg                                                                                               <  1:336509/139‑1 (MQ=255)                  aaaaCATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGaaaaa                                                                                       >  2:1357675/1‑139 (MQ=255)                      cATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCa                                                                                   <  2:125126/139‑1 (MQ=255)                          gcgcCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGc                                                                               >  2:359002/1‑139 (MQ=255)                           cgccgcTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCa                                                                              >  2:754400/1‑139 (MQ=255)                                ctgctgACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAg                                                                         <  1:182205/139‑1 (MQ=255)                                             gACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCg                                                            <  2:1302085/139‑1 (MQ=255)                                             gACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCg                                                            >  2:1136017/1‑139 (MQ=255)                                                cGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCaaa                                                                            <  2:709148/120‑1 (MQ=255)                                                cGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCaaa                                                                            >  1:709148/1‑120 (MQ=255)                                                  gAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAAc                                                       <  1:359002/139‑1 (MQ=255)                                                    aTTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCg                                                     <  2:66674/139‑1 (MQ=255)                                                               aTGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGa                                          <  1:1357675/139‑1 (MQ=255)                                                                 gAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATg                                        <  2:1554056/139‑1 (MQ=255)                                                                    cTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCg                                     <  2:694821/139‑1 (MQ=255)                                                                      gACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGaa                                                                  >  2:1542503/1‑108 (MQ=255)                                                                      gACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGaa                                                                  <  1:1542503/108‑1 (MQ=255)                                                                       aCGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑‑GGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCaaa                                                        <  1:525590/115‑1 (MQ=255)                                                                       aCGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGT‑‑GGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCaaa                                                        >  2:525590/1‑115 (MQ=255)                                                                             gACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGa                            >  1:523087/1‑139 (MQ=255)                                                                               ctctctTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAAc                                                       >  1:668251/1‑110 (MQ=255)                                                                               ctctctTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAAc                                                       <  2:668251/110‑1 (MQ=255)                                                                                            aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGa                                    >  1:1269336/1‑116 (MQ=255)                                                                                            aGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGa                                    <  2:1269336/116‑1 (MQ=255)                                                                                                       ggCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAAcc                                                      <  2:1268776/87‑1 (MQ=255)                                                                                                       ggCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAAcc                                                      >  1:1268776/1‑87 (MQ=255)                                                                                                       ggCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAg  >  1:22683/1‑139 (MQ=255)                                                                                                        gCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAg  >  1:471114/1‑138 (MQ=255)                                                                                                        gCCGATGTGGTGGTGATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAg  <  2:471114/138‑1 (MQ=255)                                                                                                                      gATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGca                 <  1:1004280/109‑1 (MQ=255)                                                                                                                      gATTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGca                 >  2:1004280/1‑109 (MQ=255)                                                                                                                       aTTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCACTCAATCCCGATAAg  >  2:353901/1‑123 (MQ=255)                                                                                                                       aTTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATa    <  2:96071/121‑1 (MQ=255)                                                                                                                       aTTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATa    >  1:96071/1‑121 (MQ=255)                                                                                                                       aTTGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAg  <  1:353901/123‑1 (MQ=255)                                                                                                                         tGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCg                                                     <  2:1046176/70‑1 (MQ=255)                                                                                                                         tGGTCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCg                                                     >  1:1046176/1‑70 (MQ=255)                                                                                                                            tCCCGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATa    >  1:1090016/1‑116 (MQ=255)                                                                                                                            tCCCGGACTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAATCCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATa    <  2:1090016/116‑1 (MQ=255)                                                                                                                              ccGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaa                <  1:701708/102‑1 (MQ=255)                                                                                                                              ccGGTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaa                >  2:701708/1‑102 (MQ=255)                                                                                                                                 gTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaat               <  2:1091492/100‑1 (MQ=255)                                                                                                                                 gTCTGGGTCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGcaat               >  1:1091492/1‑100 (MQ=255)                                                                                                                                        |                                                                                                          AGTACTGACCCGCAGTGAAAACATTGCGCCGCTGCTGACTGCACGACCGGAATTGATGGTGCATGAACTGACGATGGACTCTCTTACCGAAAGCCTGGAATGGGCCGATGTGGTGGTGATTGGTCCCGGTCTGGGCCAGCAAGAGTGGGGGAAAAAAGCACTGCAAAAAGTTGAGAATTTTCGCAAACCGATGTTGTGGGATGCCGATGCATTGAACCTGCTGGCAATCAATCCCGATAAG  >  NC_000913/4394923‑4395163

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