Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	616,676	C→T	6.9%	D840D (GAC→GAT)	entF →	enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	616,676	0	C	T	6.9%	179.7 / 2.7	58	D840D (GAC→GAT)	entF	enterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (31/23);  new base T (2/2);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

CATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGCCGTTACCCAC  >  NC_000913/616530‑616820                                                                                                                                                   |                                                                                                                                                 cATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGt                                                                                                                                                >  5:131245/1‑149 (MQ=255)  aTCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGCCGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGtt                                                                                                                                               >  6:124560/1‑149 (MQ=255)  aTCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGtt                                                                                                                                               <  1:72487/149‑1 (MQ=255)      ggTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGccc                                                                                                                                           >  7:91077/1‑149 (MQ=255)                    tgGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAAc                                                                                                                             >  3:40975/1‑149 (MQ=255)                              tctcTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTg                                                                                                                    >  5:30277/1‑148 (MQ=255)                                  tATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTa                                                                                                               >  1:217987/1‑149 (MQ=255)                                   aTCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGATGTTGCCCGCTGGCTGGAt                                                                                                                                >  8:209316/1‑131 (MQ=255)                                   aTCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGATGTTGCCCGCTGGCTGGAt                                                                                                                                <  7:209316/131‑1 (MQ=255)                                       cACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCg                                                                                                          >  8:159742/1‑149 (MQ=255)                                        aCTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTAcc                                                                                                             <  6:131245/145‑1 (MQ=255)                                              aTTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAg                                                                                                   >  8:196782/1‑149 (MQ=255)                                                  aaCTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGCAGTACCTCGGGCGCAGtgat                                                                                               >  4:200438/1‑149 (MQ=255)                                                     tGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTgatgat                                                                                            >  7:304336/1‑149 (MQ=255)                                                       gcgcAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATc                                                                                           <  7:159742/148‑1 (MQ=255)                                                              gCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGtgatga                                                                                             <  6:162869/139‑1 (MQ=255)                                                              gCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGtgatga                                                                                             >  5:162869/1‑139 (MQ=255)                                                                tATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAAt                                                                                 >  8:153827/1‑149 (MQ=255)                                                                  tctcGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTc                                                                               >  5:236717/1‑149 (MQ=255)                                                                   ctcGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTcg                                                                              >  6:138817/1‑149 (MQ=255)                                                                      ggACGCCCCTATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCgg                                                                           <  8:304336/149‑1 (MQ=255)                                                                      ggACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTcgc                                                                             <  4:40975/147‑1 (MQ=255)                                                                           ccccGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCTGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGc                                                                      <  6:236717/149‑1 (MQ=255)                                                                            cccGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACgg                                                                                                                           <  5:194528/95‑1 (MQ=255)                                                                            cccGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACgg                                                                                                                           >  6:194528/1‑95 (MQ=255)                                                                            cccGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCgg                                                                           <  5:287802/143‑1 (MQ=255)                                                                            cccGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCgg                                                                           >  6:287802/1‑143 (MQ=255)                                                                                     aCCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTAt                                                                  <  6:30277/143‑1 (MQ=255)                                                                                     aCCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATc                                                                 <  5:124560/144‑1 (MQ=255)                                                                                       cGCCAGCCGCTTTATTGCCGATCCTTTTGGCCCAGGTGAACGGATGTACCGGACCGGAGACGTGGCCCGCTGGCTGGATAACGGCGCGGGGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTgg                                                          <  3:200438/149‑1 (MQ=255)                                                                                       cGCCAGCCGCTTTAGTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTgg                                                          <  8:91077/149‑1 (MQ=255)                                                                                        gCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTggg                                                         >  2:246692/1‑149 (MQ=255)                                                                                          cAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGc                                                        >  3:146907/1‑148 (MQ=255)                                                                                           aGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCg                                                       >  5:261999/1‑148 (MQ=255)                                                                                               gCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAatc                                                  >  8:200699/1‑149 (MQ=255)                                                                                                 tttATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAatcga                                                >  6:221789/1‑149 (MQ=255)                                                                                                     ttGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATcgc                                            <  7:196782/149‑1 (MQ=255)                                                                                                        ccGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCa                                                                                                    >  5:10983/1‑90 (MQ=255)                                                                                                        ccGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCa                                                                                                    <  6:10983/90‑1 (MQ=255)                                                                                                          gATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGAt                                       >  5:17594/1‑149 (MQ=255)                                                                                                           aTCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATg                                      >  2:186606/1‑149 (MQ=255)                                                                                                            tCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGgc                                                                                                                          <  6:57521/64‑1 (MQ=255)                                                                                                            tCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGgc                                                                                                                          >  5:57521/1‑64 (MQ=255)                                                                                                            tCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGc                                     <  5:138817/149‑1 (MQ=255)                                                                                                               ttttGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAgg                                  >  5:234160/1‑149 (MQ=255)                                                                                                               ttttGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCACGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCCCGTGATGCAgg                                  >  1:381095/1‑149 (MQ=255)                                                                                                                     cccAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGc                            >  7:10158/1‑149 (MQ=255)                                                                                                                     cccAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGc                            >  7:36162/1‑149 (MQ=255)                                                                                                                       cAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCg                          >  6:148068/1‑149 (MQ=255)                                                                                                                        aGGTGAACGGATGTACCGTACCGGAGATGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTTAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAatcgatc                                              <  8:345866/128‑1 (MQ=255)                                                                                                                        aGGTGAACGGATGTACCGTACCGGAGATGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTGAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAatcgatc                                              >  7:345866/1‑128 (MQ=255)                                                                                                                            gAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGt                     <  4:146907/149‑1 (MQ=255)                                                                                                                              aCGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTc                    >  4:209424/1‑148 (MQ=255)                                                                                                                                 gATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCTGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAAc                <  6:261999/149‑1 (MQ=255)                                                                                                                                   tGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACaa              <  8:10158/149‑1 (MQ=255)                                                                                                                                       ccgtaccgGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGcc           <  1:186606/148‑1 (MQ=255)                                                                                                                                       ccgtaccgGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGcc           <  1:246692/148‑1 (MQ=255)                                                                                                                                               gagaCGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGCCGTTACCCAc  <  7:200699/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CATCCGGTGCCGCCGGGTGTGGCGGGTGATCTCTATCTCACTGGCATTCAACTGGCGCAGGGCTATCTCGGACGCCCCGATCTGACCGCCAGCCGCTTTATTGCCGATCCTTTTGCCCCAGGTGAACGGATGTACCGTACCGGAGACGTTGCCCGCTGGCTGGATAACGGCGCGGTGGAGTACCTCGGGCGCAGTGATGATCAGCTAAAAATTCGCGGGCAGCGTATCGAACTGGGCGAAATCGATCGCGTGATGCAGGCGCTGCCGGATGTCGAACAAGCCGTTACCCAC  >  NC_000913/616530‑616820

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