Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,910,662	A→T	D288E (GAT→GAA)	pstS ←	phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,910,662	0	A	T	97.1%	233.8 / ‑1.4	69	D288E (GAT→GAA)	pstS	phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base T (24/43);  total (26/43)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910570‑3910758                                                                                             |                                                                                                 gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           <  2:1289852/100‑1 (MQ=255) gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           <  2:183536/100‑1 (MQ=255) gCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgt                                                                                           <  1:1425641/100‑1 (MQ=255)  cGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTtgtg                                                                                          <  2:1553631/100‑1 (MQ=255)    tAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGa                                                                                        <  2:968191/100‑1 (MQ=255)      aTCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATc                                                                                      <  1:839081/100‑1 (MQ=255)       tCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCa                                                                                     >  1:285672/1‑100 (MQ=255)        ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    >  2:686904/1‑100 (MQ=255)        ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    <  2:565930/100‑1 (MQ=255)        ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    <  2:461769/100‑1 (MQ=255)        ccAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAg                                                                                    <  2:31463/100‑1 (MQ=255)         cAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGa                                                                                   >  2:1197041/1‑100 (MQ=255)           ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 <  1:1197046/100‑1 (MQ=255)           ggTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAAt                                                                                 >  1:404762/1‑100 (MQ=255)            gTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATg                                                                                <  2:506587/100‑1 (MQ=255)              cGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGaa                                                                              >  1:1246200/1‑100 (MQ=255)              cGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGaa                                                                              >  2:1270793/1‑100 (MQ=255)                ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  2:1331595/100‑1 (MQ=255)                ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  1:456414/100‑1 (MQ=255)                ttCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  2:844187/100‑1 (MQ=255)                ctCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACg                                                                            <  1:804625/99‑1 (MQ=255)                 tCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACga                                                                           >  2:804620/1‑99 (MQ=255)                 tCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGt                                                                           >  2:1077486/1‑100 (MQ=255)                  cGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTg                                                                          <  2:1187648/100‑1 (MQ=255)                    ccTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGt                                                                        <  1:1022215/100‑1 (MQ=255)                       gTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTaga                                                                     <  1:138974/100‑1 (MQ=255)                          tAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGt                                                                  <  2:337554/100‑1 (MQ=255)                             ccccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAAt                                                               >  2:1209921/1‑100 (MQ=255)                              cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              >  1:1546505/1‑100 (MQ=255)                              cccGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATa                                                              <  2:212507/100‑1 (MQ=255)                                cGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAgg                                                            <  1:12712/100‑1 (MQ=255)                                  gTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGcc                                                          <  1:1386925/100‑1 (MQ=255)                                    tttGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCAt                                                        <  1:447888/100‑1 (MQ=255)                                       gTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCa                                                     <  2:1325063/100‑1 (MQ=255)                                              cAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTACTAGGCCATGCATCTTCGc                                              >  2:667885/1‑100 (MQ=255)                                              cAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGc                                              <  1:238664/100‑1 (MQ=255)                                               aGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGcc                                             <  2:789413/100‑1 (MQ=255)                                                   gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  1:621461/1‑100 (MQ=255)                                                   gaagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTc                                         >  1:1389606/1‑100 (MQ=255)                                                    aagtaTTTCAGCACTTCTGTGCGTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCttt                                          <  2:1339082/98‑1 (MQ=255)                                                    aagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCt                                        >  2:1087297/1‑100 (MQ=255)                                                    aagaaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCt                                        >  2:816008/1‑100 (MQ=255)                                                       aaTTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggt                                     >  2:1103378/1‑100 (MQ=255)                                                          ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttgg                                  <  2:464122/100‑1 (MQ=255)                                                          ttCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAACAGGCCATGCATCTTCGCCTTTCtggttgg                                  <  1:21355/100‑1 (MQ=255)                                                            cAGCACTTCTGTGCTTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAGCGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                <  2:1456173/100‑1 (MQ=255)                                                            cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                >  2:1265999/1‑100 (MQ=255)                                                            cAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTc                                >  2:252211/1‑100 (MQ=255)                                                             aGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCa                               >  1:517657/1‑100 (MQ=255)                                                              gCACTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCACCTTCGCCTTTCTGGTTGGTCAg                              <  2:918779/100‑1 (MQ=255)                                                                 cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGTATCTTCGCCTTTCTGGTTGGTCAGATc                           <  2:726157/100‑1 (MQ=255)                                                                 cTTCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATc                           <  1:564109/100‑1 (MQ=255)                                                                  ttCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATcc                          <  1:43970/100‑1 (MQ=255)                                                                   tCTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCt                         >  1:171004/1‑100 (MQ=255)                                                                    cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        >  2:277967/1‑100 (MQ=255)                                                                    cTGTGCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTg                        >  2:161974/1‑100 (MQ=255)                                                                      gtgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAg                      <  1:1266004/100‑1 (MQ=255)                                                                       tgCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGc                     <  1:686906/100‑1 (MQ=255)                                                                        gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg                    <  2:464448/100‑1 (MQ=255)                                                                        gCCTTGTTCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg                    <  1:63058/100‑1 (MQ=255)                                                                               tCTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGttt             <  1:1270798/100‑1 (MQ=255)                                                                                cTGGTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGtttt            <  2:1202940/100‑1 (MQ=255)                                                                                 tagTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTg           >  1:1069641/3‑100 (MQ=255)                                                                                  ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          <  1:1571824/100‑1 (MQ=255)                                                                                  ggTTTCTTCTGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGc          <  1:816013/100‑1 (MQ=255)                                                                                         tctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTc   <  2:1246195/100‑1 (MQ=255)                                                                                         tctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGCTCAGATCCTGAGCGAAGGTTTTGCTCCAGt    >  1:1330916/1‑99 (MQ=255)                                                                                          ctGTTCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCt  >  2:425511/1‑100 (MQ=255)                                                                                          ctGTTCTCTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCt  >  1:766133/1‑100 (MQ=255)                                                                                             |                                                                                                 GCGTAATCCAGGTCGTTCGCCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCT  >  NC_000913/3910570‑3910758

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