Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,411,637	Δ1 bp	coding (453/651 nt)	yfbT ←	sugar phosphatas

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,411,634	0	C	.	90.0%	273.2 / 26.4	70	coding (456/651 nt)	yfbT	sugar phosphatas
Reads supporting (aligned to +/- strand):  ref base C (4/3);  new base . (31/32);  total (35/35)
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

GCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGA  >  NC_000913/2411536‑2411724                                                                                                   |                                                                                           gCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGcc                                                                                           >  1:442122/1‑100 (MQ=255)  ccGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGccc                                                                                          >  2:1433373/1‑100 (MQ=255)         gTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCt                                                                                  >  2:592499/1‑100 (MQ=255)          ttAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCTCGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTg                                                                                 <  1:1364738/100‑1 (MQ=255)          ttAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTg                                                                                 >  1:1155762/1‑100 (MQ=255)          ttAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTg                                                                                 <  2:904312/100‑1 (MQ=255)            aaCCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcg                                                                               <  1:16905/100‑1 (MQ=255)             aCCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcgc                                                                              >  2:1807503/1‑100 (MQ=255)                gCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGcc                                                                           <  2:452248/100‑1 (MQ=255)                       catgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGa                                                                    >  1:1377344/1‑100 (MQ=255)                        atgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGAt                                                                   <  1:458516/100‑1 (MQ=255)                         tgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATa                                                                  >  2:902880/1‑100 (MQ=255)                           acaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACg                                                                >  1:1872309/1‑100 (MQ=255)                              aCCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCAt                                                             >  1:1525226/1‑100 (MQ=255)                              aCCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCAt                                                             >  2:55657/1‑100 (MQ=255)                               cccGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATc                                                            >  2:907078/1‑100 (MQ=255)                                 cgccgcCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCan                                                          <  2:717198/100‑2 (MQ=255)                                      ccaggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCt                                                     >  1:937102/1‑100 (MQ=255)                                        aggccagAAAGCGCTCCATCGGGAGCATCTTCGACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTTCGCGCCTAACAGATACGCATCAGGTTCTgg                                                   <  2:1529808/100‑1 (MQ=255)                                         ggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                  <  2:1268710/100‑1 (MQ=255)                                         ggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                  >  2:330347/1‑100 (MQ=255)                                           ccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                >  1:1848973/1‑100 (MQ=255)                                           ccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                >  1:1524228/1‑100 (MQ=255)                                            cagAAAGCACGCCAGCGGGAGCATCTTCCACCATCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGtttt                                               >  1:221662/1‑100 (MQ=255)                                            cagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                >  1:253915/1‑100 (MQ=255)                                            cagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGtttt                                               >  2:17461/1‑100 (MQ=255)                                            cagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGtttt                                               >  1:1718020/1‑100 (MQ=255)                                            cagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                  <  1:1836532/97‑1 (MQ=255)                                             agAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTc                                              <  1:1605893/100‑1 (MQ=255)                                              naaaGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTcc                                             >  1:248937/2‑100 (MQ=255)                                                    aCGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCtt                                       <  2:973868/100‑1 (MQ=255)                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTc                                      >  2:1464038/1‑100 (MQ=255)                                                        cAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACt                                   <  1:330345/100‑1 (MQ=255)                                                           cGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGc                                <  2:759439/100‑1 (MQ=255)                                                                 cATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCCCTCGCTCAGCg                          <  2:161261/100‑1 (MQ=255)                                                                  aTCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCgg                         <  1:893761/100‑1 (MQ=255)                                                                  aTCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCgg                         >  2:1353811/1‑100 (MQ=255)                                                                   tCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGt                        <  2:501730/100‑1 (MQ=255)                                                                    cTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGtt                       >  1:1456265/1‑100 (MQ=255)                                                                    cTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGtt                       >  1:1721494/1‑100 (MQ=255)                                                                     ttCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGtt                       <  1:106742/100‑1 (MQ=255)                                                                     ttCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTa                      >  2:856226/1‑100 (MQ=255)                                                                      tCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTAc                     >  2:1549952/1‑100 (MQ=255)                                                                       ccaccaccaCACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACa                    <  1:964409/100‑1 (MQ=255)                                                                       ccaccaccaCACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGACTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACa                    <  1:1240642/100‑1 (MQ=255)                                                                       ccaccaccaCACACTCCTGCGGCGCAA‑CCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACa                    >  1:87561/1‑100 (MQ=255)                                                                        caccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACaa                   >  2:1079552/1‑100 (MQ=255)                                                                        caccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACaa                   >  2:735351/1‑100 (MQ=255)                                                                         accaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCACTTCACTCGCTCAGCGGTTACaaa                  <  1:124436/100‑1 (MQ=255)                                                                          ccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                 <  2:104543/100‑1 (MQ=255)                                                                          ccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                 >  2:1445672/1‑100 (MQ=255)                                                                          ccaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                 <  1:62979/100‑1 (MQ=255)                                                                           caccacACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAac                 >  2:386014/1‑100 (MQ=255)                                                                           caccacACACTCCTGCGGCGCAAGCCCAA‑CAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAaca                <  1:1348455/100‑1 (MQ=255)                                                                           caccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAaca                <  2:849704/100‑1 (MQ=255)                                                                             ccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACAcc              <  2:330902/100‑1 (MQ=255)                                                                                cacacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACctct           >  2:736696/1‑100 (MQ=255)                                                                                cacacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACctct           <  1:219148/100‑1 (MQ=255)                                                                                  cgcTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTg          <  2:945764/97‑1 (MQ=255)                                                                                  cacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTgg         <  2:1718031/100‑1 (MQ=255)                                                                                   acTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTgg         >  1:1590128/1‑100 (MQ=255)                                                                                   acTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGt        >  1:994560/1‑100 (MQ=255)                                                                                      ccTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCg     >  1:1590950/1‑100 (MQ=255)                                                                                      ccTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCg     <  2:1524236/100‑1 (MQ=255)                                                                                       cTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCg     <  2:253917/100‑1 (MQ=255)                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACGTCTGGTGCgg    >  1:1273419/1‑100 (MQ=255)                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg    >  1:297136/1‑100 (MQ=255)                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg    <  2:1298883/100‑1 (MQ=255)                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg    <  2:442124/100‑1 (MQ=255)                                                                                        tGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCggg   <  1:1152115/100‑1 (MQ=255)                                                                                         gcggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGa  <  1:1515455/100‑1 (MQ=255)                                                                                         gcggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGa  <  2:1469483/100‑1 (MQ=255)                                                                                                   |                                                                                           GCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGA  >  NC_000913/2411536‑2411724

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