Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,449,738	+G	coding (1106/1254 nt)	ytfR →	YtfR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,449,737	1	.	G	98.9%	235.1 / ‑5.0	90	C369G (TGC→GGC)	ytfR	YtfR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (42/47);  minor base C (1/0);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCT‑GCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCGC  >  NC_000913/4449662‑4449807                                                                             |                                                                       aGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct                                                                         <  2:4482374/76‑1 (MQ=255) aGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct                                                                         >  1:1513415/1‑76 (MQ=255)  ggAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑g                                                                       <  1:2896849/76‑1 (MQ=255)  ggAAAACCGATGGCATCAGCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑g                                                                       <  2:4509308/76‑1 (MQ=255)   gAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑gg                                                                      <  1:2216716/76‑2 (MQ=255)   gAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCACCCCCCCCGCGgct‑gc                                                                      >  2:507778/1‑76 (MQ=255)     aaaCCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑ggct                                                                    >  1:449199/1‑73 (MQ=255)     aaaCCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑ggct                                                                    <  1:4645501/76‑4 (MQ=255)     aaaCCGAGGGCATCAGCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGgct‑ggct                                                                    <  2:4543105/76‑4 (MQ=255)        ccGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACg                                                                  >  2:3684898/1‑76 (MQ=255)        ccGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCACGCCCAGCGCGGCTGGCTACg                                                                  >  2:3687456/1‑76 (MQ=255)          gATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTc                                                                >  1:2281025/1‑76 (MQ=255)           aTGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTcc                                                               >  1:3316813/1‑76 (MQ=255)           aTGGCATCAGCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTcc                                                               <  2:2141357/76‑1 (MQ=255)           aTGGCAGCAGCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTcc                                                               <  2:504520/76‑1 (MQ=255)            tGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTccc                                                              >  1:2162633/1‑76 (MQ=255)             ggCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCa                                                             <  2:3062151/76‑1 (MQ=255)              gCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTcccct                                                            >  2:1888898/1‑74 (MQ=255)              gCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCAt                                                            <  2:627958/76‑1 (MQ=255)               catcatCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCCtt                                                           >  2:936514/1‑76 (MQ=255)                gtcagcGCTGCCGGCGCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCATCGCGGCTGGCTACGTCCCAttt                                                          <  2:1595676/75‑1 (MQ=255)                atcatcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCAttt                                                          >  2:657118/1‑76 (MQ=255)                 tcgtcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTc                                                         <  2:407821/73‑1 (MQ=255)                   atcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTccc                                                       <  1:4294628/76‑1 (MQ=255)                    tcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCg                                                      <  1:1161839/76‑1 (MQ=255)                    tcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCg                                                      <  1:290290/76‑1 (MQ=255)                    tcGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCg                                                      >  2:2321434/1‑76 (MQ=255)                    gcgcgGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTccag                                                      <  2:4757641/75‑3 (MQ=255)                     cGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGc                                                     <  2:1928631/76‑1 (MQ=255)                      gCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCa                                                    >  1:2366353/1‑76 (MQ=255)                      gCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCa                                                    >  1:2334964/1‑76 (MQ=255)                      cccgcccccCCGGGGCGGGAAAAAATCCTCCCCCCCCCCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCa                                                    <  1:4662739/73‑1 (MQ=255)                       cTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCaa                                                   <  1:2447418/76‑1 (MQ=255)                        tGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCaaa                                                  >  1:2496221/1‑76 (MQ=255)                        tGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCCCaaa                                                  >  2:2671499/1‑76 (MQ=255)                           cgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGaa                                               <  2:3162923/76‑1 (MQ=255)                           cgccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGaa                                               >  2:3953093/1‑76 (MQ=255)                            gccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAAc                                              <  1:1961749/76‑1 (MQ=255)                            gccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGGCCCATTTCCCGCAAAGAAc                                              >  1:3458422/1‑76 (MQ=255)                             ccTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACa                                             >  1:4564710/1‑76 (MQ=255)                              cTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACCTCCCATTTCCCCCAAAGAACac                                            >  2:870249/1‑75 (MQ=255)                                tggTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCa                                          >  2:4183561/2‑76 (MQ=255)                                cGGTGCTGGAAAATATCATCCTCGCTCTCCAGGCCCAGCTCGGCTGGCTACGTCCCATTTCCCGCACATAACCGCa                                          <  1:473007/76‑1 (MQ=255)                                 ggTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCaa                                         <  1:656946/76‑1 (MQ=255)                                   tGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAaga                                       >  1:1336579/1‑76 (MQ=255)                                   tGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAaga                                       <  1:1486896/76‑1 (MQ=255)                                   tGCGGGAAAATATCATCCTCGCTCTACAGCCCCAGCGCGGCTGGTTAGGTCCCATTTACCGCAAATAACACGAaga                                       >  1:3247411/1‑76 (MQ=255)                                    gCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAagag                                      <  2:771621/76‑1 (MQ=255)                                      gtgAAAATATCATCCTCGCTCTCAAGGCCAAGCGCTGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGAt                                    <  2:3588062/74‑1 (MQ=255)                                      gggACAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGAt                                    <  2:2170114/76‑1 (MQ=255)                                      gggAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGAt                                    <  2:3561360/76‑1 (MQ=255)                                       ggAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAACAGAtt                                   >  2:1830165/1‑76 (MQ=255)                                        gAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTg                                  <  2:4164342/76‑1 (MQ=255)                                         aaaaTATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGc                                 >  1:3699219/1‑76 (MQ=255)                                          aaaTATCTTCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGcc                                >  1:554332/1‑76 (MQ=255)                                          aaaTATCATCCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATGTCCCGCAAAGAACAGCAAGAGATTGcc                                >  2:1705851/1‑76 (MQ=255)                                            atatCATCCTCGCTCTCCAGGCCCAGCGCGGCTCGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGa                              >  2:151393/1‑76 (MQ=255)                                            atatCAGCCGCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGa                              <  2:1679590/76‑1 (MQ=255)                                              atcatcCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAAc                            <  1:2321435/76‑1 (MQ=255)                                               tcatcCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGCTTGCCGCAcc                           >  2:4408375/1‑75 (MQ=255)                                               tcatcCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACg                           <  2:2001391/76‑1 (MQ=255)                                               tcatcCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACg                           <  1:4024898/76‑1 (MQ=255)                                                 atcCTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCt                         <  1:4608093/76‑1 (MQ=255)                                                  tcCTCGCGCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCtt                        <  1:2824351/76‑1 (MQ=255)                                                   ccTCGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCttt                       <  1:4656888/76‑1 (MQ=255)                                                     ccGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTAt                     <  1:18731/75‑1 (MQ=255)                                                      cGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATc                    >  1:1208196/1‑76 (MQ=255)                                                      cGCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATc                    <  2:1076514/76‑1 (MQ=255)                                                       gCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATcc                   >  1:3563344/1‑76 (MQ=255)                                                       gCTCTCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCACGCAACGAACAGCAAGAGATTGCCGAACGCGTTAtac                   >  1:3589673/1‑74 (MQ=255)                                                        ctctCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCg                  <  1:1454694/76‑1 (MQ=255)                                                        ctctCCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCg                  >  1:1786891/1‑76 (MQ=255)                                                        ctctCCAGGCCCAGCGCGGCTGGCTACGGCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCg                  >  1:1933670/1‑76 (MQ=255)                                                         tctcCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGc                 <  1:3693445/76‑1 (MQ=255)                                                          ctcCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGcc                >  1:328211/1‑76 (MQ=255)                                                          ctcCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGcc                >  1:4123523/1‑76 (MQ=255)                                                          ctcCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAACAACAGCAAGAGATTGCCCAACCCCTTAcccccc                >  2:874040/1‑76 (MQ=255)                                                           tcCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCa               <  1:1672746/76‑1 (MQ=255)                                                           tcCAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCa               >  2:2370628/1‑76 (MQ=255)                                                            ccAGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAg              <  2:2451954/76‑1 (MQ=255)                                                             cAGGCCCAGCGCGGCTGGCTACGCCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGc             <  2:3569736/76‑1 (MQ=255)                                                              aGGCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCt            <  2:1407192/76‑1 (MQ=255)                                                               ggCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCtt           <  1:1394102/76‑1 (MQ=255)                                                               ggCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCtt           <  2:679772/76‑1 (MQ=255)                                                               ggCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCtt           <  1:1038434/76‑1 (MQ=255)                                                                gCCCAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTg          <  2:546298/76‑1 (MQ=255)                                                                 cccAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCCGCTTgg         >  2:4731997/1‑76 (MQ=255)                                                                 cccAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTgg         >  1:234186/1‑76 (MQ=255)                                                                  ccAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGc        <  1:3153633/76‑1 (MQ=255)                                                                   cAGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCa       >  1:1014753/1‑76 (MQ=255)                                                                    aGCGCGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCAt      <  1:2671501/76‑1 (MQ=255)                                                                     gcgcGGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCAtt     <  1:791817/76‑1 (MQ=255)                                                                      cgcgGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTc    >  1:4861477/1‑76 (MQ=255)                                                                      cgcgGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTc    <  1:4294877/76‑1 (MQ=255)                                                                      cgcgGCTGGCGAGGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGACCGCTTTATCCGCCAGCTTGGCATTc    <  2:3641366/76‑1 (MQ=255)                                                                       gcgGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCg   >  1:1897124/1‑76 (MQ=255)                                                                       gcgGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCg   >  1:443481/1‑76 (MQ=255)                                                                       gcgGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTCGCATTCg   >  2:2564316/1‑76 (MQ=255)                                                                        cgGCTGGCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCGc  >  1:2801041/1‑76 (MQ=255)                                                                             |                                                                       AGGAAAACCGATGGCATCATCGCTGCCGCCTCGGTGCGGGAAAATATCATCCTCGCTCTCCAGGCCCAGCGCGGCT‑GCTACGTCCCATTTCCCGCAAAGAACAGCAAGAGATTGCCGAACGCTTTATCCGCCAGCTTGGCATTCGC  >  NC_000913/4449662‑4449807

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