Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	805,315	C→T	9.0%	L605L (CTG→TTG)	ybhJ →	aconitase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	805,315	0	C	T	9.0%	107.5 / 2.0	47	L605L (CTG→TTG)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base T (2/2);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

GATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAT  >  NC_000913/805167‑805460                                                                                                                                                     |                                                                                                                                                  taTTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGc                                                                                                                                                   <  1:173828/148‑1 (MQ=255)  aTTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGCCTGGCGGGGAATGTCAGCGAGCt                                                                                                                                                  >  1:540868/1‑149 (MQ=255)  aTTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGAAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGAAAAGCGGCTGCTGAGCTGGAAAGTCAGCGTCTGGGGGGGAATGTCAGCGAGca                                                                                                                                                  >  1:196968/1‑148 (MQ=255)       ttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACag                                                                                                                                             >  1:691717/1‑149 (MQ=255)       ttccGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACag                                                                                                                                             >  1:469473/1‑149 (MQ=255)       ttccGGTGAAACCTCTTCTTATCGTTAAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACag                                                                                                                                             >  1:655789/1‑149 (MQ=255)       ttccGGTGAAACCTCTTCTTAGCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGAGCCCGGTTATGTTAGCAGAAGTAAAGAGACTGCTGGCATGGAAGATCAGCATCTGGCGGGGAATGTCAGCGAGGTGGCag                                                                                                                                             >  1:47731/1‑149 (MQ=255)              gAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGtt                                                                                                                                      <  1:99582/149‑1 (MQ=255)                   ctctTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcgc                                                                                                                                 >  1:1359/1‑149 (MQ=255)                        cttCTCGTTCAAATCCGATTGGTATGGCGGAGTTTATCCAGAAACGACGCGATCCCGGTAATGCAAGATTTATTTAAGCAACAGCGCTTGAGTAAAAGCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                            <  1:771995/149‑1 (MQ=255)                          tATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGGGGTGTTTGCGCGCATTaa                                                                                                                          <  1:574449/149‑1 (MQ=255)                            tCGTACACACCCGCTTGGTCTGTCGCAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGc                                                                                                                        <  2:611384/149‑1 (MQ=255)                                   aaTCCGATTGGTCTGGCGGAGTTTACGCTGTCTCGCAGCGATCCCGGTTCTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGACAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGGATTAAGCAGGttac                                                                                                                 >  1:229973/1‑147 (MQ=255)                                       cGATTGGTCTGGCGGAGTTTACCCTGTCCCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGt                                                                                                             <  2:555684/149‑1 (MQ=255)                                         aTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGTTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCa                                                                                                           >  2:129036/1‑149 (MQ=255)                                         aTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGTTGACAGAGGTGTTTGCGCGCATTAAGCAGATCGCTGGTCa                                                                                                           >  2:544912/1‑149 (MQ=255)                                                tGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGCGCCGCCGGCGGTGTGTGTGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                                                                    >  1:604105/1‑149 (MQ=255)                                                tGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCat                                                                                                    >  1:711474/1‑149 (MQ=255)                                                  gCGGAGTTTACCCTGTCTCGCCGCGCTCCCGGTTATGTTATCAGAAGCAAAGCGACTGCTGAACTGGAAAATCAGCGTATGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCatct                                                                                                  >  1:200793/1‑147 (MQ=255)                                                  gCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCatat                                                                                                  >  1:522021/1‑149 (MQ=255)                                                      aGTTTACCCTGTCCCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                              <  1:360869/149‑1 (MQ=255)                                                      aGTTTACCCTGCCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                              <  1:223163/149‑1 (MQ=255)                                                              cTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCa                                                                                      >  2:308265/1‑149 (MQ=255)                                                                     gccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGGTTGCGCGCAGTACGAAGATTGCTGGTCAGGCGCATATTGATCCGCTGCAAACTGac                                                                               >  2:128306/1‑148 (MQ=255)                                                                      ccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaaa                                                                              <  1:639049/149‑1 (MQ=255)                                                                      ccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaaa                                                                              <  1:779459/149‑1 (MQ=255)                                                                      ccgcGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaaa                                                                              <  1:420664/149‑1 (MQ=255)                                                                                  ttATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATgg                                                                  <  2:724768/149‑1 (MQ=255)                                                                                       ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGGCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTGTAt                                                             >  1:187832/1‑149 (MQ=255)                                                                                          gCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGTTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCACACTGAAATTGGCAGCATGGTCTATgtg                                                          <  1:544912/149‑3 (MQ=255)                                                                                          gCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGTTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCg                                                          <  1:129036/149‑1 (MQ=255)                                                                                             gAAGTAAAGCTACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTg                                                       <  1:217558/149‑1 (MQ=255)                                                                                             gAAGTAAAGCGCCTGCTGAGCGGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAAGTGGCAGCATGGTCTATGCGGTg                                                       <  2:90483/149‑1 (MQ=255)                                                                                              gaGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCCGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGa                                                      <  2:194830/148‑1 (MQ=255)                                                                                              aaGTAAAGCGACTGCTGGGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGa                                                      <  1:558159/149‑1 (MQ=255)                                                                                                 tAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAc                                                   <  2:273342/149‑1 (MQ=255)                                                                                                   aaGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                 <  1:678276/149‑1 (MQ=255)                                                                                                       gACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGGATTAAGCAGATTGCTGGTCACGAGCATATTGAGCCGCTGAAAACTGAAATTGGCAGCATGGTCTATGCGGTGAACCCAgggc                                             >  2:337440/1‑147 (MQ=255)                                                                                                             tgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                       >  1:304496/1‑149 (MQ=255)                                                                                                                 cTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGCAGGTTCTgc                                   >  2:653131/1‑149 (MQ=255)                                                                                                                        aTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGaa                            >  2:658746/1‑149 (MQ=255)                                                                                                                           aGCGTCGGGGGGGGGAGGGCAGCGGTCGTACCGAGGTGGTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCAGATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTTTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAg                         <  1:698153/149‑1 (MQ=255)                                                                                                                            gCGTCTGGGGGAGCAGGTCCGCGAGCTGGCAGTGGTGGTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAgg                        <  2:611302/149‑1 (MQ=255)                                                                                                                            gCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGCGGTTTGCGCGCATTAATCAGATTGCTGGTTAGGAGATTATTTGTCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAgg                        <  2:44/149‑1 (MQ=255)                                                                                                                                 tGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGgcggcg                   <  2:304496/149‑1 (MQ=255)                                                                                                                                   gCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGCACAGgcggcgcg                 >  1:196203/1‑147 (MQ=255)                                                                                                                                     ggggAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCt               <  2:1359/149‑1 (MQ=255)                                                                                                                                                  gtgCTGATAGAGGTGTTGGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt  <  2:160846/147‑1 (MQ=255)                                                                                                                                                  gAGCGGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGCTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGTTGCCAGCGTGTGAt  <  2:604105/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAT  >  NC_000913/805167‑805460

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