Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,361,323	T→A	Q31L (CAG→CTG)	rhmD ←	L‑rhamnonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,361,323	0	T	A	97.3%	255.8 / ‑1.3	76	Q31L (CAG→CTG)	rhmD	L‑rhamnonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (38/35);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

CGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGAT  >  NC_000913/2361231‑2361417                                                                                             |                                                                                               cGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTg                                                                                         >  1:522058/1‑100 (MQ=255)  gTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGa                                                                                        <  2:708337/100‑1 (MQ=255)  gTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGa                                                                                        <  2:1112169/100‑1 (MQ=255)  gTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGa                                                                                        <  2:1173197/100‑1 (MQ=255)  gTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGa                                                                                        <  2:1283820/100‑1 (MQ=255)   ttAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGAt                                                                                       <  2:846392/100‑1 (MQ=255)     aaTGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATaa                                                                                     <  1:1708683/100‑1 (MQ=255)      aTGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAAt                                                                                    <  2:1543826/100‑1 (MQ=255)         ccAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCAc                                                                                 <  1:284218/100‑1 (MQ=255)          cAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCAcc                                                                                <  1:648309/100‑1 (MQ=255)          cAAATGACTGGCGTGACTGCTCGTAAACGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCAcc                                                                                >  2:1183301/1‑100 (MQ=255)              tGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCAccgcc                                                                             <  1:972953/99‑1 (MQ=255)                aCTGGCGTGACTGCTCGTAATCGCGGTATTAACTCATCGGTGTGGCAATATAATCGTCCATCCAGTGATTCGCCCCCAGGTCGTGATAATCTCCgccgcc                                                                          <  1:799147/100‑1 (MQ=255)                      gTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGcc                                                                    >  2:683163/1‑100 (MQ=255)                        gACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCtt                                                                  <  2:1706811/100‑1 (MQ=255)                          cTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCtttt                                                                <  1:233015/100‑1 (MQ=255)                           tGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCttttt                                                               >  1:1106079/1‑100 (MQ=255)                             cTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                             >  2:30888/1‑100 (MQ=255)                             cTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                             >  2:427074/1‑100 (MQ=255)                             cTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                             <  1:1225770/100‑1 (MQ=255)                             cTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                             <  1:992617/100‑1 (MQ=255)                             cTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCACCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTct                                                             >  2:867615/1‑100 (MQ=255)                              tCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctg                                                            <  2:1269481/100‑1 (MQ=255)                              tCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctg                                                            <  2:1088019/100‑1 (MQ=255)                               cGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                           >  2:143418/1‑100 (MQ=255)                               cGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgc                                                           >  2:1175669/1‑100 (MQ=255)                                gTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgct                                                          <  2:1251846/100‑1 (MQ=255)                                gTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgct                                                          >  2:1221639/1‑100 (MQ=255)                                gTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCAACGCCGCCAGCGCCTTTTTctgct                                                          >  1:1759541/1‑100 (MQ=255)                                 tAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCGGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTctgctg                                                         >  1:1363304/1‑100 (MQ=255)                                    tCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCg                                                      <  2:839991/100‑1 (MQ=255)                                      gcgGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCa                                                    <  1:539128/100‑1 (MQ=255)                                         gTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCg                                                 <  1:269463/100‑1 (MQ=255)                                           aTTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGcc                                               >  1:1093605/1‑100 (MQ=255)                                             ttACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAg                                             <  1:187203/100‑1 (MQ=255)                                                cTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATTACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTaa                                          >  1:893943/1‑100 (MQ=255)                                                  cATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAc                                        >  1:868698/1‑100 (MQ=255)                                                  cATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAc                                        >  2:953547/1‑100 (MQ=255)                                                      ggggTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAAcc                                       >  1:333102/1‑97 (MQ=255)                                                        ggTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGgcg                                  <  2:288842/100‑1 (MQ=255)                                                            gCAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGaa                              >  2:1412362/1‑100 (MQ=255)                                                             cAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAc                             >  2:779263/1‑100 (MQ=255)                                                             cAATATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCCTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAc                             >  2:1130223/1‑100 (MQ=255)                                                              aaTATGATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAAcc                            <  1:907932/100‑1 (MQ=255)                                                                tatCATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCCCGAACCt                           >  2:492144/1‑99 (MQ=255)                                                                   gATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGttt                       <  2:1738165/100‑1 (MQ=255)                                                                   gATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGttt                       >  2:902105/1‑100 (MQ=255)                                                                   gATCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGttt                       <  1:683162/100‑1 (MQ=255)                                                                    aTCGTCTATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTa                      >  1:1877501/1‑100 (MQ=255)                                                                         cTATCCAGTGATTCGCCCCCAGGTCGTGATAAGCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAttt                  <  1:649980/99‑1 (MQ=255)                                                                          tATCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAAttttt                <  1:137294/100‑1 (MQ=255)                                                                           aTCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg               >  1:300412/1‑100 (MQ=255)                                                                           aTCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg               >  1:1901229/1‑100 (MQ=255)                                                                           aTCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg               >  1:506778/1‑100 (MQ=255)                                                                           aTCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTg               <  2:1084420/100‑1 (MQ=255)                                                                            tCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg              >  1:1241287/1‑100 (MQ=255)                                                                            tCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg              <  1:293672/100‑1 (MQ=255)                                                                            tCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg              >  2:1425801/1‑100 (MQ=255)                                                                            tCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg              <  2:899794/100‑1 (MQ=255)                                                                            tCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg              <  1:584816/100‑1 (MQ=255)                                                                            tCCAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTgg              >  1:1203642/1‑100 (MQ=255)                                                                              ccgTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTa            >  1:1827081/3‑100 (MQ=255)                                                                              cAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTa            >  2:988034/1‑100 (MQ=255)                                                                              cAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTa            >  2:779658/1‑100 (MQ=255)                                                                              cAGTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTa            >  1:1824784/1‑100 (MQ=255)                                                                                gTGATTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTTTTTTTTTTTTCCGCCAGTAAACCCGGCCCGAACCTGTTTAATTTTTGGt             >  1:1325127/1‑97 (MQ=255)                                                                                   aTTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTc       >  2:1545384/1‑100 (MQ=255)                                                                                   aTTCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTc       >  2:1545192/1‑100 (MQ=255)                                                                                     tCGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCGCCGCCAGTAGACCAGGCGCGAACCAGTTTAATTTTTGGTAGGGTCat     >  1:1310317/1‑100 (MQ=255)                                                                                      cGCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatg    >  2:296895/1‑100 (MQ=255)                                                                                       gCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatga   <  2:1260019/100‑1 (MQ=255)                                                                                       gCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatga   <  2:481864/100‑1 (MQ=255)                                                                                       gCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatga   <  2:481747/100‑1 (MQ=255)                                                                                       gCCCCCAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatga   >  2:1349287/1‑100 (MQ=255)                                                                                        cccccAGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCatgat  <  2:177925/100‑1 (MQ=255)                                                                                          cccAGGTCGCGATAATCACCGCCGCCAGCGCCTTTTTTTTCTGTTTTATTTTCTTTAAACCAGGCGCGAAACTGTCTAATTTTTGGTAGGGCCatgat  >  1:82713/1‑98 (MQ=255)                                                                                             |                                                                                               CGTTAATGCCAAATGACTGGCGTGACTGCTCGTAATCGCGGTATTTACTCATCGGGGTGGCAATATGATCGTCTATCCAGTGATTCGCCCCCTGGTCGTGATAATCACCGCCGCCAGCGCCTTTTTCTGCTGTCGCACCGCCAGTAAACCAGGCGCGAACCTGTTTAATTTTTGGTAGGGTCATGAT  >  NC_000913/2361231‑2361417

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