Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,059,293	T→C	G11G (GGA→GGG)	rpiA ←	ribose 5‑phosphate isomerase, constitutive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,059,293	0	T	C	100.0%	230.8 / NA	65	G11G (GGA→GGG)	rpiA	ribose 5‑phosphate isomerase, constitutive
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (35/30);  total (35/30)

GGCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCATCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATAAAATA  >  NC_000913/3059179‑3059425                                                                                                                   |                                                                                                                                     ggCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                              >  2:858276/1‑139 (MQ=255) ggCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGAGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                              >  2:364066/1‑139 (MQ=255)                     tCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGcc                                                                                          >  2:1232373/1‑139 (MQ=255)                      cATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCt                                                                                         <  1:858276/139‑1 (MQ=255)                            aCCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTa                                                                                   >  2:1360339/1‑139 (MQ=255)                                aGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCttttt                                                                                                                         <  2:1321180/97‑1 (MQ=255)                                aGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCttttt                                                                                                                         >  1:1321180/1‑97 (MQ=255)                                aGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGaa                                                                               >  2:1199822/1‑139 (MQ=255)                                aGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGaa                                                                               >  2:421426/1‑139 (MQ=255)                                  cgcgTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAAt                                                                             >  1:775914/1‑139 (MQ=255)                                       cAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcac                                                                        <  1:1199822/139‑1 (MQ=255)                                         aTAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      <  1:421426/139‑1 (MQ=255)                                         aTAAAGTGTGCGGCGGTGGAACCTGTACATACACCAACAATGGTGCCCGGCTGAACATAGTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATGCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTGacac                                                                      <  1:733767/139‑1 (MQ=255)                                               tgtgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATta                                                                >  1:1324372/1‑139 (MQ=255)                                                 tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                        >  2:154381/1‑113 (MQ=255)                                                 tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                        <  1:154381/113‑1 (MQ=255)                                                 tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcac                                                                        >  1:1269457/1‑129 (MQ=255)                                                 tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcac                                                                        <  2:1269457/129‑1 (MQ=255)                                                 tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      >  1:809683/1‑131 (MQ=255)                                                 tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      <  2:1316870/131‑1 (MQ=255)                                                 tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      <  2:809683/131‑1 (MQ=255)                                                 tgCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      >  1:1316870/1‑131 (MQ=255)                                                  gcggcgGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcaca                                                                       >  1:290932/1‑129 (MQ=255)                                                  gcggcgGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcaca                                                                       <  2:290932/129‑1 (MQ=255)                                                           gAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCt                                                                                                                             >  1:1145332/1‑66 (MQ=255)                                                           gAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCt                                                                                                                             <  2:1145332/66‑1 (MQ=255)                                                            aaCCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAg                                                   <  1:15585/139‑1 (MQ=255)                                                                tGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGAt                                               >  1:41795/1‑139 (MQ=255)                                                                        cacCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      >  1:841689/1‑108 (MQ=255)                                                                        cacCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      <  2:841689/108‑1 (MQ=255)                                                                         acCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcaca                                                                       <  1:1305453/106‑1 (MQ=255)                                                                         acCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcaca                                                                       >  2:1305453/1‑106 (MQ=255)                                                                            aacaaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      <  2:1334094/104‑1 (MQ=255)                                                                            aacaaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      >  1:1334094/1‑104 (MQ=255)                                                                            aacaaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGa                                                <  1:1254761/126‑1 (MQ=255)                                                                            aacaaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGa                                                >  2:1254761/1‑126 (MQ=255)                                                                              caaTGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCATACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTgg                                 <  2:1324372/139‑1 (MQ=255)                                                                                 tGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                        >  1:1218628/1‑81 (MQ=255)                                                                                 tGGTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCtg                                                                                        <  2:1218628/81‑1 (MQ=255)                                                                                  ggTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcac                                                                        >  1:1140555/1‑96 (MQ=255)                                                                                  ggTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcac                                                                        <  2:1140555/96‑1 (MQ=255)                                                                                  ggTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAg                             <  2:41795/139‑1 (MQ=255)                                                                                   gTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGt                                                                                                <  2:894483/71‑1 (MQ=255)                                                                                   gTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGt                                                                                                >  1:894483/1‑71 (MQ=255)                                                                                   gTGCCGGGCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGc                            >  1:539263/1‑139 (MQ=255)                                                                                         ggCTGAACATACTGAAGTGCCGCCCACCCTACTGCTTTTTTGAATTCATCCGGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTg                      >  2:946781/1‑139 (MQ=255)                                                                                              aaCATACTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTaa                 >  1:467957/1‑139 (MQ=255)                                                                                                   aCTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGAccc                                       >  2:245883/1‑112 (MQ=255)                                                                                                   aCTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGAccc                                       <  1:245883/112‑1 (MQ=255)                                                                                                   aCTGAAGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCAt            <  1:1360339/139‑1 (MQ=255)                                                                                                        aGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                              >  2:1339665/1‑36 (MQ=255)                                                                                                        aGTGCCGCCCACCCTACTGCTTTTTTCAATTCATcc                                                                                                              <  1:1339665/36‑1 (MQ=255)                                                                                                        aGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcac                                                                        <  2:1377394/74‑1 (MQ=255)                                                                                                        aGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTcac                                                                        >  1:1377394/1‑74 (MQ=255)                                                                                                        aGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      >  1:759260/1‑76 (MQ=255)                                                                                                        aGTGCCGCCCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCacac                                                                      <  2:759260/76‑1 (MQ=255)                                                                                                           gccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATGGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCAt            >  2:1086643/1‑131 (MQ=255)                                                                                                           gccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCAt            <  1:1086643/131‑1 (MQ=255)                                                                                                            ccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCa                                                                   >  2:1246720/1‑75 (MQ=255)                                                                                                            ccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCa                                                                   <  1:1246720/75‑1 (MQ=255)                                                                                                            ccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCcgc                         >  2:555155/1‑117 (MQ=255)                                                                                                            ccgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCcgc                         <  1:555155/117‑1 (MQ=255)                                                                                                             cgccCACCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATAAAata  >  2:448222/1‑139 (MQ=255)                                                                                                               cccACCCTACTGCTTTTTTCAATTCATCCTGCGTCATg                                                                                                     <  2:663384/38‑1 (MQ=38)                                                                                                               cccACCCTACTGCTTTTTTCAATTCATCCTGCGTCATg                                                                                                     >  1:663384/1‑38 (MQ=255)                                                                                                                   |                                                                                                                                     GGCCCCTTCAATCTGGCCTTTCATTGTACCGAGCGCGTCAATAAAGTGTGCGGCGGTGGAACCTGTACCTACACCAACAATGGTGCCGGGCTGAACATACTGAAGTGCCGCCCATCCTACTGCTTTTTTCAATTCATCCTGCGTCATGATCGTTTCGCCTGTGGTATGAAATTTCACACGCATTATATACAAAAAAAGCGATTCAGACCCCGTTGGCAAGCCGCGTGGTTAACTCATCCATAAAATA  >  NC_000913/3059179‑3059425

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