Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,747,584	A→C	S152R (AGC→CGC)	ydhS →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,747,584	0	A	C	97.9%	157.3 / ‑2.7	47	S152R (AGC→CGC)	ydhS	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (20/26);  total (20/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GGGCGAATATTTCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCAGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCGTGTAACGTGGGTATTATGGGAACATCC  >  NC_000913/1747436‑1747730                                                                                                                                                     |                                                                                                                                                   gggCGAATATTTCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGc                                                                                                                                                  <  1:668067/151‑1 (MQ=255)  ggCGAATATTTCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCcgcg                                                                                                                                                 >  2:29839/1‑151 (MQ=255)   gCGAATATTTCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCcgcg                                                                                                                                                 >  2:90047/1‑150 (MQ=255)     gAATATTTCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCcgcg                                                                                                                                                 >  2:900072/1‑148 (MQ=255)       atatTTCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACg                                                                                                                                            >  2:859965/1‑151 (MQ=255)         atTTCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGtt                                                                                                                                          <  2:891664/151‑1 (MQ=255)          tttCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGttt                                                                                                                                         >  1:497910/1‑151 (MQ=255)                cgATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGAttt                                                                                                                                    >  1:443549/1‑150 (MQ=255)                 gATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAg                                                                                                                                  <  2:825754/151‑1 (MQ=255)                    cAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCgg                                                                                                                               <  2:502763/151‑1 (MQ=255)                        ttttACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGAt                                                                                                                           >  1:561963/1‑151 (MQ=255)                              gACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCAc                                                                                                                     <  2:92423/151‑1 (MQ=255)                               aCTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACg                                                                                                                    >  2:735487/1‑151 (MQ=255)                                  agtagACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGt                                                                                                                 >  2:703405/1‑151 (MQ=255)                                   gtagACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTc                                                                                                                <  2:628571/151‑1 (MQ=255)                                    tagACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTc                                                                                                                >  2:2956/1‑150 (MQ=255)                                     agACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCAc                                                                                                              <  2:845460/151‑1 (MQ=255)                                      gACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACg                                                                                                             >  2:815443/1‑151 (MQ=255)                                           ggCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTgg                                                                                                        <  2:588847/151‑1 (MQ=255)                                             cACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGcc                                                                                                      >  1:916694/1‑151 (MQ=255)                                              aCGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCt                                                                                                     >  2:836327/1‑151 (MQ=255)                                              aCGACAGCAAAAATTTGCAGGGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATCTTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCAACGGGTCACGTCGTGCCt                                                                                                     >  1:308240/1‑151 (MQ=255)                                               cGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCtg                                                                                                    <  1:205189/151‑1 (MQ=255)                                                 aCAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCtgat                                                                                                  <  2:806013/151‑1 (MQ=255)                                                          aTTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATgaagaagaa                                                                                         <  2:385434/151‑1 (MQ=255)                                                                    gCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCGGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCg                                                                                 <  1:563422/149‑1 (MQ=255)                                                                       gTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAAc                                                                              <  1:660191/149‑1 (MQ=255)                                                                         ttATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTa                                                                           <  2:803495/150‑1 (MQ=255)                                                                         ttATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTa                                                                           <  1:784812/150‑1 (MQ=255)                                                                          tATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTAt                                                                          <  2:330216/150‑1 (MQ=255)                                                                           aTGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTAttt                                                                        >  1:60058/1‑151 (MQ=255)                                                                                       aGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTg                                                            >  1:481286/1‑151 (MQ=255)                                                                                       aGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTg                                                            >  1:784055/1‑151 (MQ=255)                                                                                         gTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGt                                                          >  1:810763/1‑151 (MQ=255)                                                                                         gTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGt                                                          <  1:200250/151‑1 (MQ=255)                                                                                                    gCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTGACATGTGGCCTGATGAAGAAGAAGCAACCTGAACGTATTTTCCCAGCCCGTGGTCAgtctggatg                                               <  1:773327/151‑9 (MQ=255)                                                                                                     cAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATgg                                              <  2:745941/151‑1 (MQ=255)                                                                                                             aaaTGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCaa                                        >  2:35963/1‑149 (MQ=255)                                                                                                             aaaTGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAg                                      <  1:619454/151‑1 (MQ=255)                                                                                                               aTGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGt                                    <  2:568859/151‑1 (MQ=255)                                                                                                                       gTCATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCgt                            >  2:80957/1‑151 (MQ=255)                                                                                                                         cATGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCgtgt                          <  1:924037/151‑1 (MQ=255)                                                                                                                          aTGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCgtgt                          >  2:259687/1‑150 (MQ=255)                                                                                                                          aTGCTCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCGTGTa                         <  2:642776/151‑1 (MQ=255)                                                                                                                              tCGCTACAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCGTGTAACg                      >  1:872683/1‑150 (MQ=255)                                                                                                                                    cAAATCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCGTGTAACGTGGGtat               <  1:90045/151‑1 (MQ=255)                                                                                                                                     aaaTCAGGATTTACCCAGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGAGGGAAGCAAAGGTCGATGCGTGTAACGTGGGtat               <  1:925292/150‑1 (MQ=255)                                                                                                                                       aTCAGGATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCGTGTAACGTGGGtattat            >  2:1000/1‑151 (MQ=255)                                                                                                                                           ggATTTACCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCGTGTAACGTGGGTATTATGGGa        <  1:862180/151‑1 (MQ=255)                                                                                                                                                 aCCCCGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCGTGTAACGTGGGTATTATGGGAACATcc  <  2:152308/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGGCGAATATTTCCGCGATCAATTTTTACGACTAGTAGACCAGGCACGACAGCAAAAATTTGCAGTGGCTGTTTATGAATCATGCCAGGTTACCGATCTGCAAATTACAAATGCTGGCGTCATGCTCGCTACAAATCAGGATTTACCCAGCGAGACGTTTGATTTAGCGGTGATCGCCACGGGTCACGTCTGGCCTGATGAAGAAGAAGCAACCCGAACGTATTTTCCCAGCCCGTGGTCAGGCCTGATGGAAGCAAAGGTCGATGCGTGTAACGTGGGTATTATGGGAACATCC  >  NC_000913/1747436‑1747730

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