Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,430,216	A→G	D118G (GAC→GGC)	cycA →	serine/alanine/glycine/cycloserine:H(+)symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,430,216	0	A	G	98.2%	192.6 / ‑4.6	55	D118G (GAC→GGC)	cycA	serine/alanine/glycine/cycloserine:H(+)symporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (22/32);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTT  >  NC_000913/4430079‑4430348                                                                                                                                          |                                                                                                                                     ggCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCg                                                                                                                                    >  2:131808/1‑140 (MQ=255)   cAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCgtg                                                                                                                                  <  1:378651/140‑1 (MQ=255)   cAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCgtg                                                                                                                                  <  1:373033/140‑1 (MQ=255)      tGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCgtggtg                                                                                                                               <  1:829403/140‑1 (MQ=255)      tGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCgtggtg                                                                                                                               <  1:834302/140‑1 (MQ=255)      tGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCgtggtg                                                                                                                               >  2:376417/1‑140 (MQ=255)            aaTTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATc                                                                                                                         <  1:131808/140‑1 (MQ=255)            aaTTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATc                                                                                                                         <  1:376417/140‑1 (MQ=255)              ttGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCAc                                                                                                                       <  1:824072/140‑1 (MQ=255)              ttGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCAc                                                                                                                       >  2:335905/1‑140 (MQ=255)                   gctTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCtt                                                                                                                  <  1:750784/140‑1 (MQ=255)                      ttCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATg                                                                                                               <  1:487639/140‑1 (MQ=255)                      ttCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATg                                                                                                               >  2:352236/1‑140 (MQ=255)                        cGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCt                                                                                                             <  1:18439/140‑1 (MQ=255)                         gAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTc                                                                                                            <  1:374339/140‑1 (MQ=255)                              tGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTc                                                                                                       >  1:182182/1‑140 (MQ=255)                                gAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTg                                                                                                     <  2:593581/140‑1 (MQ=255)                                        aTCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCg                                                                                             >  1:853635/1‑140 (MQ=255)                                                  gACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGa                                                                                   <  1:682103/140‑1 (MQ=255)                                                     ttCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTg                                                                                >  2:368902/1‑140 (MQ=255)                                                        gCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGt                                                                             <  1:69043/140‑1 (MQ=255)                                                           tCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGc                                                                          >  1:369014/1‑140 (MQ=255)                                                             cGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCt                                                                        >  1:10005/1‑140 (MQ=255)                                                                 ttACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCt                                                                    <  1:10788/140‑1 (MQ=255)                                                                   aCTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTgg                                                                  <  2:404135/140‑1 (MQ=255)                                                                    cTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGc                                                                 <  1:673858/140‑1 (MQ=255)                                                                       gggCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGt                                                              >  2:365289/1‑140 (MQ=255)                                                                        ggCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTg                                                             <  1:352236/140‑1 (MQ=255)                                                                             tGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAg                                                         >  1:241851/1‑139 (MQ=255)                                                                             tGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAg                                                         <  2:241851/139‑1 (MQ=255)                                                                                   ggATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGtgc                                                      <  1:205100/136‑1 (MQ=255)                                                                                   ggATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGtgc                                                      >  2:205100/1‑136 (MQ=255)                                                                                    gATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAg                                                         <  2:757976/132‑1 (MQ=255)                                                                                    gATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAg                                                         >  1:757976/1‑132 (MQ=255)                                                                                      tatTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTgctgctgct                                               <  1:365289/140‑1 (MQ=255)                                                                                         ttCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGAc                                            >  1:583160/1‑140 (MQ=255)                                                                                           cACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGc                                          <  1:535463/140‑1 (MQ=255)                                                                                           cACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGc                                          <  2:182182/140‑1 (MQ=255)                                                                                             ccGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTc                                        <  2:10005/140‑1 (MQ=255)                                                                                             ccGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTc                                        <  1:335905/140‑1 (MQ=255)                                                                                              cGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCa                                       >  1:304131/1‑140 (MQ=255)                                                                                                gCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCaaa                                     >  1:573593/1‑139 (MQ=255)                                                                                                gCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAAt                                     >  2:290054/1‑140 (MQ=255)                                                                                                gCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAAt                                     <  1:290054/140‑1 (MQ=255)                                                                                                       ttACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCa                              <  1:816922/140‑1 (MQ=255)                                                                                                        tACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCAc                             <  2:369014/140‑1 (MQ=255)                                                                                                               tCTGCTGGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGaaa                      >  2:447643/1‑140 (MQ=255)                                                                                                                   ctgGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGt                  >  1:632546/1‑140 (MQ=255)                                                                                                                   ctgGGTTGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGt                  >  1:696123/1‑140 (MQ=255)                                                                                                                         tGTAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTg            <  2:386752/140‑1 (MQ=255)                                                                                                                           tAACCGGTATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTgag          >  1:105272/1‑140 (MQ=255)                                                                                                                                ggTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGa     <  1:405631/140‑1 (MQ=255)                                                                                                                                  tATGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGt   >  1:93423/1‑140 (MQ=255)                                                                                                                                   aTGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGtt  <  1:368902/140‑1 (MQ=255)                                                                                                                                   aTGGCAGGCGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGtt  <  2:304131/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                     GGCAATGGGGGAATTGCTGCTTTCGAATCTGGAATACAAATCTTTTAGTGACTTCGCTTCCGATTTACTCGGGCCGTGGGCAGGATATTTCACCGGCTGGACTTACTGGTTCTGCTGGGTTGTAACCGGTATGGCAGACGTGGTGGCGATCACGGCTTATGCTCAGTTCTGGTTCCCCGATCTCTCCGACTGGGTCGCCTCGCTGGCGGTGATAGTGCTGCTGCTGACGCTCAATCTCGCCACCGTGAAAATGTTCGGTGAGATGGAGTT  >  NC_000913/4430079‑4430348

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: ‑