Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,925,932	(T)7→6	100%	coding (493/657 nt)	yobB →	C‑N hydrolase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,925,926	0	T	.	79.7%	171.7 / 53.9	69	coding (487/657 nt)	yobB	C‑N hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (5/7);  new base . (22/33);  total (27/42)
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

ACCGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGAT  >  NC_000913/1925829‑1926024                                                                                                  |                                                                                                   acCGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTTCTGGCGTCGGCCCGTCGCTTACAAttt                                                                                                  >  2:2591703/1‑100 (MQ=255)   cGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAGCCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttt                                                                                                <  2:211803/100‑1 (MQ=255)   cGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttt                                                                                                <  2:2286766/100‑1 (MQ=255)   cGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttt                                                                                                <  2:327351/100‑1 (MQ=255)    gCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAtttttt                                                                                               >  1:1803456/1‑100 (MQ=255)     cAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttttc                                                                                              <  2:1633712/100‑2 (MQ=255)     cAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttttc                                                                                              >  2:1073162/1‑99 (MQ=255)       aGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAttttttcac                                                                                            <  2:418337/100‑4 (MQ=255)         gTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTcaca                                                                                          <  2:2224155/100‑1 (MQ=255)         gTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACAt                                                                                         >  1:1399743/1‑100 (MQ=255)          tatggatatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGCACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTCTTCACAt                                                                                         <  1:532779/100‑1 (MQ=255)            tggatatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAg                                                                                      <  1:1549430/100‑1 (MQ=255)            tggatatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAg                                                                                      >  1:1370211/1‑100 (MQ=255)              gatatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGt                                                                                     >  2:525322/1‑100 (MQ=255)              gatatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTa                                                                                    >  1:2471762/1‑100 (MQ=255)                tatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAAtttttttttattatt                                                                                      >  1:1869063/1‑89 (MQ=255)                tatggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGCGAGCACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACt                                                                                  >  1:92284/1‑100 (MQ=255)                  tggaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGTACTc                                                                                 <  2:102084/100‑1 (MQ=255)                    gaCCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCg                                                                                <  1:2594468/98‑1 (MQ=255)                       ccGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTg                                                                            >  1:400484/1‑99 (MQ=255)                         gagatGTTCACTGTTTACTACAGGCCAATTTCGCTGAGAACCCGACCTGCTTGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTg                                                                            <  1:2115644/93‑1 (MQ=255)                         gACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCg                                                                         <  1:1073158/100‑1 (MQ=255)                         gACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCg                                                                         <  1:1952088/100‑1 (MQ=255)                         gACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCg                                                                         <  2:711865/100‑1 (MQ=255)                         gACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCg                                                                         >  1:1473343/1‑100 (MQ=255)                           ccTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTg                                                                       >  1:1034249/1‑100 (MQ=255)                             tcttCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCt                                                                     <  2:884573/98‑1 (MQ=255)                              gTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCGGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGc                                                                      >  1:1385287/1‑98 (MQ=255)                                tCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCtaga                                                                  >  1:884570/1‑97 (MQ=255)                                tCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGAt                                                                  >  2:2559974/1‑100 (MQ=255)                                 cACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGGTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATg                                                                 >  1:1410571/1‑100 (MQ=255)                                    tGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCa                                                              >  2:1830584/1‑100 (MQ=255)                                      tctactacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCaaa                                                            <  1:2104898/98‑1 (MQ=255)                                         actacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCaaaag                                                          >  2:2104912/1‑97 (MQ=255)                                          ctacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGcc                                                        <  2:1679280/100‑1 (MQ=255)                                          ctacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGcc                                                        >  2:91120/1‑100 (MQ=255)                                          ctacAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGcc                                                        >  2:1127935/1‑100 (MQ=255)                                             cAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGt                                                     >  2:1395139/1‑100 (MQ=255)                                              aGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTg                                                    <  2:377164/100‑1 (MQ=255)                                               gggCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTgg                                                   <  2:1499351/100‑1 (MQ=255)                                                ggCAATGTCTTGTAGAACGCGACCTGCTGGCGTCTGCTCGTTGCTTACAA‑TTTTTTCACATCAGTATTTGATTGCCGTGTTGATGTCTAATGCCCGTGGt                                                  <  1:1151314/100‑1 (MQ=255)                                                ggCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTgg                                                   <  2:1221093/100‑1 (MQ=255)                                                ggCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGt                                                  >  2:2280945/1‑100 (MQ=255)                                                 gCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTa                                                 <  2:2380239/100‑1 (MQ=255)                                                    aTGTCTGGGGGAACCCGACCTGCTGGCGTCGGTCCGTCGCTTACAA‑ATTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTATCa                                              <  1:2456480/100‑1 (MQ=255)                                                     tGTCTGGGAGAACCCTACCTGCTGGCGGCGGCCCGTCGCTTACAATTTTTTTCACATCAGTCCTCGATTGCCGTGCTGATGGCAAACGCCCGTGGTAACa                                              <  1:1601059/100‑1 (MQ=255)                                                     tGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAg                                             <  2:2471780/100‑1 (MQ=255)                                                      gTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGTACTCGATTGCCGGGCTTATGGCAAATGCCCGTGGTAACAg                                             <  1:156020/100‑1 (MQ=255)                                                      gTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGt                                            >  1:1077405/1‑100 (MQ=255)                                                         tGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCa                                         <  2:2102481/100‑1 (MQ=255)                                                          gggAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCa                                         <  2:966008/100‑1 (MQ=255)                                                          gggAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCAc                                        <  2:147344/100‑1 (MQ=255)                                                              gaACGCGATCTGCTGGGGTCGGCCAGTCGCGTAC‑AGTGGTTACACATCAGTACTCGATTGCCGTGCTGATATCAAATGCCCGTGGTAACAGTgaa                                         <  1:778282/95‑3 (MQ=255)                                                                aCCCGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCCCATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATgg                                   >  1:265756/1‑99 (MQ=255)                                                                 cccGACCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGa                                 <  2:1801564/100‑1 (MQ=255)                                                                     aCCTGCTGGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATAGGATGa                              >  2:2201361/1‑99 (MQ=255)                                                                      ccTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTT‑TTTTTTTTTCTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAat                            >  1:1457069/1‑100 (MQ=255)                                                                          ctgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATgg                        <  2:300457/100‑1 (MQ=255)                                                                           tgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGgt                       <  2:1617011/100‑1 (MQ=255)                                                                           tgGCGTCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGgt                       <  2:340524/100‑1 (MQ=255)                                                                                tCGGCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGtct                  <  2:168298/100‑1 (MQ=255)                                                                                  ggCCCGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCa                <  2:1389007/100‑1 (MQ=255)                                                                                  ggCCCGTCGCTTAC‑AGTTTTTGCACATCCGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGGAACAGTGGTCTATGGGATGAATATGgtcg                     <  1:464845/95‑1 (MQ=255)                                                                                    cccGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCa                <  1:1937147/98‑1 (MQ=255)                                                                                     cgatctcTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCg             <  1:804990/94‑1 (MQ=255)                                                                                     cgatctcTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCg             <  1:804923/94‑1 (MQ=255)                                                                                      cGTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCg             <  1:2454454/99‑1 (MQ=255)                                                                                       gTCGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCg             <  1:2453032/98‑1 (MQ=255)                                                                                       gTCGCTTACAA‑TTTTTTCACATCAGTACTAGATGGTCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCg             <  1:1750248/98‑1 (MQ=255)                                                                                         cGCTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTcg         <  2:1693712/100‑1 (MQ=255)                                                                                          gCTTACAA‑TTTTTTCACATCAATACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTcgc        <  1:1578716/100‑1 (MQ=255)                                                                                           cTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCgagatcg       >  2:804926/1‑94 (MQ=255)                                                                                           cTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCgagatcg       >  2:804993/1‑94 (MQ=255)                                                                                           cTTACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTcgcg       <  2:1903666/100‑1 (MQ=255)                                                                                             tACAA‑TTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGcc     <  2:1883047/100‑1 (MQ=255)                                                                                               cAATTTTTTTCACATCAGTACTCGATTGCCATGATGATGGCAAATGCCTGTGGTAACAGTGCACCATGGGACGGATATGGTCGTCTCATCGTTCGCGCCg    >  2:2451956/1‑100 (MQ=255)                                                                                                 aTTTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGAt  >  2:516193/1‑100 (MQ=255)                                                                                                 aTTTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGAt  >  2:737748/1‑100 (MQ=255)                                                                                                  |                                                                                                   ACCGCAAGGTATGGATATGGACCCGACCTGTTCACTGTTTACTACAGGGCAATGTCTGGGAGAACCCGACCTGCTGGCGTCGGCCCGTCGCTTACAATTTTTTTCACATCAGTACTCGATTGCCGTGCTGATGGCAAATGCCCGTGGTAACAGTGCACTATGGGATGAATATGGTCGTCTCATCGTTCGCGCCGAT  >  NC_000913/1925829‑1926024

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