Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	632,971	C→T	100%	A169A (GCG→GCA)	ybdH ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	632,971	0	C	T	98.8%	265.2 / ‑4.3	81	A169A (GCG→GCA)	ybdH	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (37/43);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

TTATTGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCCGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAGCGCC  >  NC_000913/632876‑633064                                                                                                |                                                                                              ttattGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCca                                                                                           >  1:1833870/1‑100 (MQ=255)   attGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCcagc                                                                                         >  1:155087/1‑100 (MQ=255)   attGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCcag                                                                                          <  1:1979727/99‑1 (MQ=255)    ttGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCcagca                                                                                        <  2:2321892/100‑1 (MQ=255)    ttGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCcagca                                                                                        <  1:651761/100‑1 (MQ=255)     tGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCcagcag                                                                                       >  1:1517530/1‑100 (MQ=255)     tGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCcagcag                                                                                       >  1:1915047/1‑100 (MQ=255)      gatggccAGTCGGTCGGTGAGCGGCAACGTTTGTGGTTGCGGAGCCAGCACCACCGCTGCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGa                                                                                      >  2:1992487/6‑100 (MQ=255)      gATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGa                                                                                      >  1:2367110/1‑100 (MQ=255)      gATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGa                                                                                      >  2:475897/1‑100 (MQ=255)      gATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGa                                                                                      <  1:1897760/100‑1 (MQ=255)      gATCCCCAGTCGCACGGTTAGCGGAAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATAACATTTCGCCAGCGTGTCACCGATCTCTGCcagca                                                                                        >  1:824649/1‑98 (MQ=255)        tCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGata                                                                                    >  2:325417/1‑100 (MQ=255)              gTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATAttgttg                                                                              >  2:1053459/1‑100 (MQ=255)              gTCGCACGGTTAGCGGCAACGTTTCGGGTTGCGGAGCCAGCTCCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATAttgttg                                                                              <  1:1844183/100‑1 (MQ=255)                cGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCg                                                                            <  1:888453/100‑1 (MQ=255)                cGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCg                                                                            >  2:1040642/1‑100 (MQ=255)                cGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCg                                                                            >  2:1040772/1‑100 (MQ=255)                 gCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGc                                                                             >  2:1016116/1‑98 (MQ=255)                  cACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGt                                                                          <  2:959540/100‑1 (MQ=255)                   aCGGTTAGCGGTAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCg                                                                            <  2:2386769/97‑1 (MQ=255)                   aCGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGt                                                                          >  2:786007/1‑99 (MQ=255)                   aCGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTg                                                                         >  2:786342/1‑100 (MQ=255)                     ggTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCa                                                                       <  1:326031/100‑1 (MQ=255)                      gTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCAt                                                                      <  2:279346/100‑1 (MQ=255)                           cGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAgg                                                                 >  1:659638/1‑100 (MQ=255)                            ggCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTTCGGTGCATTGAGGa                                                                >  1:681452/1‑100 (MQ=255)                               aaCGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATaa                                                             <  2:1217598/100‑1 (MQ=255)                                 cGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAAtc                                                           <  1:1007984/100‑1 (MQ=255)                                    ttCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTcc                                                        <  2:84625/100‑1 (MQ=255)                                    ttCTGGTTGCGGAGCCAGCACCACCGCTTCAGACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTcc                                                        <  2:864504/100‑1 (MQ=255)                                     tCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCg                                                       <  1:1242095/100‑1 (MQ=255)                                       tGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGt                                                     <  1:50027/100‑1 (MQ=255)                                       tGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGt                                                     <  2:49548/100‑1 (MQ=255)                                         gTTGCGGAGCCAGCACCCCCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGCGCATTCAGGATACTCGCCGGTTc                                                   >  2:633896/1‑100 (MQ=255)                                          ttGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTcc                                                  <  2:326527/100‑1 (MQ=255)                                          tcgccgggCCAGCACCACCGCTTCCTACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTcc                                                  <  2:575691/93‑1 (MQ=255)                                           tGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTcca                                                 <  1:1118401/100‑1 (MQ=255)                                            gCGGAGCCAGCACCTCCGATTCATACCATTTCGTCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTccac                                                <  1:841407/100‑1 (MQ=255)                                            gCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTccac                                                <  1:662761/100‑1 (MQ=255)                                             cgatctcAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTccacc                                               <  1:1053890/94‑1 (MQ=255)                                             cGGAGCCAGCACCACCGCTTCATACCACTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTccacc                                               <  1:644720/100‑1 (MQ=255)                                               gAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAg                                             <  2:964310/100‑1 (MQ=255)                                               gAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAg                                             <  2:1106317/100‑1 (MQ=255)                                                 gCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCa                                           >  2:580412/1‑100 (MQ=255)                                                 gCCAGCACCACCGCTTCATACCATTTCGCCAGCGCGTCACCGATCCATGCCAGCAGATATGGTTCCGGTACATCGAGGATGATCTCTTATTCCACCAGc                                            >  2:538318/1‑99 (MQ=255)                                                  ccAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCAc                                          <  1:270485/100‑1 (MQ=255)                                                   cAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAgatcg                                         >  2:1053886/1‑96 (MQ=255)                                                   cAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCAcc                                         >  2:2165895/1‑100 (MQ=255)                                                   cAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCAcc                                         <  2:2218079/100‑1 (MQ=255)                                                    aGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGTTGCATTGAGGATAATCTCCGGTTTCACCAGCa                                           >  2:708147/1‑97 (MQ=255)                                                    aGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGATTATAATCTCCGGTTCCACCAGCACCa                                        <  2:807720/100‑1 (MQ=255)                                                    aGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCa                                        <  1:2104345/100‑1 (MQ=255)                                                    aGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCa                                        <  2:2279953/100‑1 (MQ=255)                                                      caccacCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATa                                      >  1:171391/1‑100 (MQ=255)                                                        ncaccGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATaaa                                    >  2:260292/2‑100 (MQ=255)                                                         caccGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATaaaa                                   >  1:1014800/1‑100 (MQ=255)                                                         caccGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCCGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATaca                                    >  1:1837076/1‑97 (MQ=255)                                                          accGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAAt                                  <  1:885494/100‑1 (MQ=255)                                                             gCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTgg                               <  1:94685/100‑1 (MQ=255)                                                             gCTTCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCAACAGCACCATAAAATTg                                >  1:1369892/1‑99 (MQ=255)                                                               ttCATACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGcg                             <  1:88802/100‑1 (MQ=255)                                                                   tACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGcgtcgt                         <  2:132416/100‑1 (MQ=255)                                                                   tACCATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGcgtcgt                         <  1:415875/100‑1 (MQ=255)                                                                     ccATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGtcgtcg                       >  1:2346166/1‑100 (MQ=255)                                                                     ccATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGtcgtcg                       >  1:608844/1‑100 (MQ=255)                                                                     ccATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGtcgtcg                       >  1:1646492/1‑100 (MQ=255)                                                                      cATTTCGCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGa                      >  2:2235999/1‑100 (MQ=255)                                                                         ttCGCCAGCGTGTCACCTATCCCTGCCAGCATATATTGTTGCGGTGCATTTAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGaaa                    <  1:999100/99‑1 (MQ=255)                                                                           cGCCAGCGTGTCACCGATCCCTACCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAAtc                 >  2:659955/1‑100 (MQ=255)                                                                            gCCAGCGTGTCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAAtct                <  2:261567/100‑1 (MQ=255)                                                                             ccAGCGTGTCACCGATCCCTGCCAGCAGATATTGTCGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAAtct                <  2:261399/99‑1 (MQ=255)                                                                                  gtgtCACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAAt          <  2:1865687/100‑1 (MQ=255)                                                                                  gtgtCACCGATCCCTGACAGCCGATATTGTTGTTGTGCATTGAGGATAATCTCCGGTTCCACTAGCATCATAAAATTGGCGTCAGCGAAAATCTCATAAt          <  2:608921/100‑1 (MQ=255)                                                                                      cACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAg      >  2:1153184/1‑100 (MQ=255)                                                                                      cACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAg      >  1:1905167/1‑100 (MQ=255)                                                                                      cACCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAg      >  2:2233075/1‑100 (MQ=255)                                                                                       aCCGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAgc     <  2:81435/100‑1 (MQ=255)                                                                                        ccGATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAgcg    >  2:684462/1‑100 (MQ=255)                                                                                          gATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAGCGcc  <  1:1468044/100‑1 (MQ=255)                                                                                          gATCCCTGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAGCGcc  >  1:2200681/1‑100 (MQ=255)                                                                                                |                                                                                              TTATTGATCCCCAGTCGCACGGTTAGCGGCAACGTTTCTGGTTGCGGAGCCAGCACCACCGCTTCATACCATTTCGCCAGCGTGTCACCGATCCCCGCCAGCAGATATTGTTGCGGTGCATTGAGGATAATCTCCGGTTCCACCAGCACCATAAAATTGGCGTCGTCGAAAATCTCATAATGCAGCGCC  >  NC_000913/632876‑633064

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