Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,358,040	T→C	H193R (CAC→CGC)	cadA ←	lysine decarboxylase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,358,040	0	T	C	95.0%	195.0 / ‑0.0	60	H193R (CAC→CGC)	cadA	lysine decarboxylase 1
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (33/24);  total (35/25)
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.72e-01

TGGTGCTGCCTGCTGGAGCAGAGTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGTGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCATGTGACCAGGAGTA  >  NC_000913/4357906‑4358174                                                                                                                                       |                                                                                                                                       tGGTGCTGCCTGCTGGAGCAGAGTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATcc                                                                                                                                   <  1:456436/140‑1 (MQ=255)         cctgcTGGAGCAGAGTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAgaga                                                                                                                           >  1:213/1‑140 (MQ=255)         cctgcTGGAGCAGAGTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGAGATATACTGTTCAGTTGCTTTGTGTGGACCGCTGCGATCgagcagctg                                                                                                                           >  2:358863/1‑137 (MQ=255)          ctgctgGAGCAGAGTACATACCTACAATTTTGTTCTCAGGGGAAGTACTGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGaa                                                                                                                          <  2:237568/140‑1 (MQ=255)            gctgGAGCAGAGTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAAcc                                                                                                                        <  1:409951/140‑1 (MQ=255)            gctgGAGCAGAGTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTATGGACCACTGCGATCCAGCAGAGAAcc                                                                                                                        <  2:683171/140‑1 (MQ=255)                   cAGAGTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCa                                                                                                                 >  1:373211/1‑140 (MQ=255)                      agTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGAt                                                                                                              <  1:481313/140‑1 (MQ=255)                          cATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTGCTTTGTGTTGACGACTGCGATCCAGGAGAGGACCCAGTTCAGATACTg                                                                                                          >  2:216361/1‑140 (MQ=255)                            tACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGaa                                                                                                        <  2:112102/140‑1 (MQ=255)                                aacaaTTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATgg                                                                                                    <  1:511134/140‑1 (MQ=255)                                aacaaTTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATgg                                                                                                    <  1:450683/140‑1 (MQ=255)                                  cagttttGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGaa                                                                                                  >  2:463905/4‑140 (MQ=255)                                  caaTTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGaa                                                                                                  >  1:573712/1‑140 (MQ=255)                                  caaTTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGaa                                                                                                  >  1:112231/1‑140 (MQ=255)                                            cGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGAtt                                                                                        <  2:373211/140‑1 (MQ=255)                                                      gTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTAtt                                                                              <  2:723855/140‑1 (MQ=255)                                                         ccGTTGGTCACCATGTAGCTGCTGTCTGCGTCAAACACGGGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCgg                                                                           <  2:229384/140‑1 (MQ=255)                                                           gTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGTGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGAc                                                                         >  1:668276/1‑140 (MQ=255)                                                             tGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCa                                                                       >  2:410791/1‑140 (MQ=255)                                                             tGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGAAACTCAAATGGAAATATGAGATTTCATGGTATTCGGACCa                                                                       >  2:269911/1‑140 (MQ=255)                                                              ggTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCaa                                                                      >  2:57770/1‑140 (MQ=255)                                                               gTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGTGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCaaa                                                                     >  2:718399/1‑140 (MQ=255)                                                                 cACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGa                                                                   >  1:96167/1‑140 (MQ=255)                                                                 cACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGa                                                                   >  2:171757/1‑140 (MQ=255)                                                                 cACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGa                                                                   <  2:112231/140‑1 (MQ=255)                                                                 cACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGa                                                                   <  1:57770/140‑1 (MQ=255)                                                                         aGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAg                                                           >  1:468546/1‑140 (MQ=255)                                                                         aGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAg                                                           >  2:527707/1‑140 (MQ=255)                                                                              cGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAg                                                      >  1:138386/1‑140 (MQ=255)                                                                              cGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTAATCGGACCAAGGAAATCATAGAACAg                                                      >  1:284459/1‑140 (MQ=255)                                                                                gTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGTAAATATCAGATTGCATGGTATTCGGGCCAACGAAATCATAGAACAGGc                                                    >  2:205017/1‑140 (MQ=255)                                                                                gTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGc                                                    >  1:52378/1‑140 (MQ=255)                                                                                gTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGc                                                    >  1:390655/1‑140 (MQ=255)                                                                                gTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGc                                                    <  1:615191/140‑1 (MQ=255)                                                                                             aCGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCt                                       >  1:486224/1‑140 (MQ=255)                                                                                                   gcgaTATACTGTTCTGCTTCTTTGTGTGGACCACTGTGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTg                                 <  1:718399/140‑1 (MQ=255)                                                                                                   gcgaTATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTg                                 <  1:410791/140‑1 (MQ=255)                                                                                                   gcgaTATACTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTg                                 <  1:600407/140‑1 (MQ=255)                                                                                                           cTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTTCCGGGCTTTTCTGGAATGCAg                         <  2:138386/140‑1 (MQ=255)                                                                                                           cTGTTCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAg                         <  1:527707/140‑1 (MQ=255)                                                                                                              ttCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAAt                             <  1:448405/133‑1 (MQ=255)                                                                                                              ttCTGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAAt                             >  2:448405/1‑133 (MQ=255)                                                                                                                 tGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGCACATGTTGCCTACCGGGCTTTTCTGGAATGCAGTACCGc                   >  2:335528/1‑140 (MQ=255)                                                                                                                 tGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAg                         >  1:468684/1‑134 (MQ=255)                                                                                                                 tGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAg                         <  2:468684/134‑1 (MQ=255)                                                                                                                 tGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGc                   >  2:426862/1‑140 (MQ=255)                                                                                                                 tGCTTCTTTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGc                   >  1:283406/1‑140 (MQ=255)                                                                                                                   cttcttTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGccc                 >  2:53840/1‑140 (MQ=255)                                                                                                                   cttcttTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGccc                 >  1:83273/1‑140 (MQ=255)                                                                                                                   cttcttTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGAGCATTTCTGAAATGCAGTACCGccc                 >  1:260583/1‑140 (MQ=255)                                                                                                                   cttcttTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAAGCTACCTGCCGGGCGTTTCTGTAATGCAGTACCGccc                 >  1:189951/1‑140 (MQ=255)                                                                                                                    ttcttTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCa                >  1:123778/1‑140 (MQ=255)                                                                                                                     tcttTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCAt               <  2:401170/140‑1 (MQ=255)                                                                                                                     tcttTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCAt               <  2:468546/140‑1 (MQ=255)                                                                                                                     tcttTGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCAt               <  2:213/140‑1 (MQ=255)                                                                                                                        ttGTGTGGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCAtgtg            >  2:444993/1‑140 (MQ=255)                                                                                                                             tgGACCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAAGGCAGTACCGCCCATGTGACCAg       >  1:323728/1‑140 (MQ=255)                                                                                                                                aCCACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCATGTGACCAGGAg    >  1:106379/1‑140 (MQ=255)                                                                                                                                  cACTGCGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCATGTGACCAGGAGTa  <  2:562842/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       TGGTGCTGCCTGCTGGAGCAGAGTACATACCAACAATTTTGTTCGCAGTGGAAGTACCGTTGGTCACCATGTAGCTGCGGTCTGCGTTAAAGACGCGAGCGATATACTGTTCTGCTTCTTTGTGTGGACCACTGTGATCCAGCAGAGAACCCAGTTCAGATACTGAAATGGAAATATCAGATTTCATGGTATTCGGACCAAAGAAATCATAGAACAGGCTACCTACCGGGCTTTTCTGGAATGCAGTACCGCCCATGTGACCAGGAGTA  >  NC_000913/4357906‑4358174

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