Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,058,570	A→T	intergenic (‑160/+96)	serA ← / ← rpiA	D‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,058,570	0	A	T	100.0%	229.3 / NA	63	intergenic (‑160/+96)	serA/rpiA	D‑3‑phosphoglycerate dehydrogenase/ribose 5‑phosphate isomerase, constitutive
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (26/37);  total (26/37)

AAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAG  >  NC_000913/3058473‑3058664                                                                                                  |                                                                                               aaaaTAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGatttttt                                                                                              <  2:2450301/100‑1 (MQ=255) aaaaTAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGatttttt                                                                                              <  2:1718808/100‑1 (MQ=255)     taataTTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTaga                                                                                          <  1:127908/100‑1 (MQ=255)     taataTTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTaga                                                                                          <  2:577613/100‑1 (MQ=255)     taataTTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTaga                                                                                          <  1:577613/100‑1 (MQ=255)     taataTTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTaga                                                                                          <  2:127908/100‑1 (MQ=255)      aataTTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTagag                                                                                         <  2:2514486/100‑1 (MQ=255)        tatTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCa                                                                                       >  2:2260233/1‑100 (MQ=255)         atTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCaa                                                                                      >  2:2022547/1‑100 (MQ=255)           tGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATc                                                                                    <  2:1778623/100‑1 (MQ=255)            gCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCg                                                                                   <  2:2211852/100‑1 (MQ=255)             cGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGa                                                                                  >  2:2118181/1‑100 (MQ=255)              tccgatctGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGAc                                                                                 <  2:1237324/93‑1 (MQ=255)              ggCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGAc                                                                                 <  2:1778737/100‑1 (MQ=255)               gCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACa                                                                                <  2:210583/100‑1 (MQ=255)               gCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACa                                                                                <  1:210583/100‑1 (MQ=255)                  atatGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACagat                                                                             >  2:192882/1‑97 (MQ=255)                  atatGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACagat                                                                             >  1:192882/1‑97 (MQ=255)                   tatGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTg                                                                            <  2:2533289/100‑1 (MQ=255)                     tGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGcc                                                                          <  2:2039919/100‑1 (MQ=255)                     tGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGcc                                                                          >  2:2464032/1‑100 (MQ=255)                       aaCGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTg                                                                        <  2:33793/100‑1 (MQ=255)                       aaCGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTg                                                                        <  1:33793/100‑1 (MQ=255)                         cGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGt                                                                      >  1:155926/1‑100 (MQ=255)                         cGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGt                                                                      <  2:396809/100‑1 (MQ=255)                         cGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGt                                                                      <  1:396809/100‑1 (MQ=255)                         cGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGt                                                                      >  2:155926/1‑100 (MQ=255)                              gCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAg                                                                 >  2:2369116/1‑100 (MQ=255)                               cGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGc                                                                >  2:1335353/1‑100 (MQ=255)                                  cgcgATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTg                                                             <  2:2278529/100‑1 (MQ=255)                                      aTGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGacac                                                         <  2:1773107/100‑1 (MQ=255)                                         ttGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATg                                                      <  2:360130/100‑1 (MQ=255)                                         ttGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATg                                                      <  1:360130/100‑1 (MQ=255)                                          tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGt                                                     >  1:490465/1‑100 (MQ=255)                                          tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGt                                                     >  2:490465/1‑100 (MQ=255)                                          tGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGt                                                     >  2:2032100/1‑100 (MQ=255)                                           gAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTc                                                    <  2:842507/100‑1 (MQ=255)                                           gAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTc                                                    <  1:842507/100‑1 (MQ=255)                                                 atatCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCaaa                                              <  2:2183763/100‑1 (MQ=255)                                                  tatCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAAt                                             <  2:1108537/100‑1 (MQ=255)                                                       cGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCATATTTAAt                                        <  2:1352633/100‑1 (MQ=255)                                                        gCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATg                                       <  2:1891493/100‑1 (MQ=255)                                                           aaGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAag                                    <  2:1052776/100‑1 (MQ=255)                                                            aGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAaga                                   >  2:203943/1‑100 (MQ=255)                                                            aGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAaga                                   >  1:203943/1‑100 (MQ=255)                                                            aGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAaga                                   <  2:1870491/100‑1 (MQ=255)                                                               tAATTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGaa                                <  2:1867265/100‑1 (MQ=255)                                                                 aTTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtt                              >  2:792376/1‑100 (MQ=255)                                                                 aTTCTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtt                              >  1:792376/1‑100 (MQ=255)                                                                    cTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAttttt                           >  2:1665903/1‑100 (MQ=255)                                                                    cTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAttttt                           >  1:306954/1‑100 (MQ=255)                                                                    cTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAttttt                           <  2:2354791/100‑1 (MQ=255)                                                                    cTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAttttt                           <  2:1156616/100‑1 (MQ=255)                                                                    cTTGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAttttt                           >  2:306954/1‑100 (MQ=255)                                                                     ttGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtttttt                          >  2:1637520/1‑100 (MQ=255)                                                                     ttGGATGAATATGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAAtttttt                          <  2:1157054/100‑1 (MQ=255)                                                                             atatGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACgggg                  >  2:2545096/1‑100 (MQ=255)                                                                               atGCTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGa                >  2:1169954/1‑100 (MQ=255)                                                                                  cTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGt             >  1:661002/1‑100 (MQ=255)                                                                                  cTGGTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGt             >  2:661002/1‑100 (MQ=255)                                                                                    ggTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTc           <  2:588592/100‑1 (MQ=255)                                                                                    ggTTTAGAGGATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTc           <  1:588592/100‑1 (MQ=255)                                                                                            ggATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCa   >  2:1235170/1‑100 (MQ=255)                                                                                            ggATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCa   <  2:1025452/100‑1 (MQ=255)                                                                                             gATTTTTTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGCAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAg  <  2:1811350/100‑1 (MQ=255)                                                                                                  |                                                                                               AAAATAATATTGCGGCAATATGAACGTTTGCGTCGCGATGTTGAGGAAATATCACGCAAAGGTAATTCTTGGATGAATATGCTGGTTTAGAGGATTTATTTAGAGCAATCGACAATTGCCTGGTAAAAGCGTGACACATGTCACCAAATTTAATGAAGAGAATTTTTTTAACGGGGGAGGTTCCCCCGTCAG  >  NC_000913/3058473‑3058664

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