Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,510,928	A→T	64.4%	L177Q (CTG→CAG)	fecE ←	ferric citrate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,510,928	0	A	T	64.4%	23.8 / 34.9	35	L177Q (CTG→CAG)	fecE	ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (5/7);  new base T (12/9);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

TTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACT  >  NC_000913/4510674‑4511182                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                ttACGGCCTGTGCAATCTACCTCAGTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCATGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGtgatt                                                                                                                                                                                                                                                 <  1:2568/271‑1 (MQ=255)                        aTTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCCCAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCCACCGCCCCTGGTCCACCCGTGGAATGAAATCAAGACAggggggttgggg                                                                                                                                                                                                                         >  2:16458/1‑265 (MQ=255)                             gCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGTGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGACGTTCAGCACAGCAACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGCTGATTGATATCAAGATAGGTGGTAGGCtcatcac                                                                                                                                                                                                                    >  1:153285/1‑271 (MQ=255)                                        ggccatgatccGGTTCGGGGGGTTTTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACg                                                                                                                                                                                                         <  2:153285/263‑1 (MQ=255)                                              aTACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTa                                                                                                                                                                                                   <  1:16458/272‑1 (MQ=255)                                                       cGGGGTGTATTTCTGCTTCCACGCTGAATACTGTTCTCAGCAATCATGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCACAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCCTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGGTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTGATCAAGTAATACAAc                                                                                                                                                                                                          >  1:353840/1‑256 (MQ=255)                                                                          cACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTa                                                                                                                                                                                                   <  1:222135/244‑1 (MQ=255)                                                                                              gCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGt                                                                                                                                                                                                    >  2:209508/1‑223 (MQ=255)                                                                                              gCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTT‑‑TCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGtgat                                                                                                                                                                                                                                                  >  2:414075/1‑175 (MQ=255)                                                                                               cAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTa                                                                                                                                                                                                   <  1:209508/223‑1 (MQ=255)                                                                                               cAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTT‑‑TCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGtgatt                                                                                                                                                                                                                                                 <  1:414075/175‑1 (MQ=255)                                                                                                aaTCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACTCCCAGTTGATCGCAGTACCGGCGAGCCTGATTAAGGTCGTGCAGCACAGCGACAACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTCGTACACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCATATACTACAACGAGCGTAGTCTTGGCCCGGACCATT‑CCACAAACTGCACGCCGACCCTGACCGCAGGa                                                                                                                                                >  1:133523/1‑272 (MQ=255)                                                                                                               aCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTc                                                                                                                                                                                                                                                             >  2:452834/1‑149 (MQ=255)                                                                                                               aCCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGtgattg                                                                                                                                                                                                                                                >  2:57607/1‑162 (MQ=255)                                                                                                                ccTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTcc                                                                                                                                                                                                                                                            <  1:452834/149‑1 (MQ=255)                                                                                                                ccTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGtgattga                                                                                                                                                                                                                                               <  1:57607/162‑1 (MQ=255)                                                                                                                         gtgtGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGtgattga                                                                                                                                                                                                                                               >  1:190173/1‑153 (MQ=255)                                                                                                                          tgtgCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGtgattgat                                                                                                                                                                                                                                              <  2:190173/153‑1 (MQ=255)                                                                                                                                                 cGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAAc                                                                                                                                                                                                          >  2:113645/1‑166 (MQ=255)                                                                                                                                                  gTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACg                                                                                                                                                                                                         <  1:113645/166‑1 (MQ=255)                                                                                                                                                                                     tAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGt                                                                             <  1:259321/254‑1 (MQ=255)                                                                                                                                                                                           gATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAAt                                                                                                                                                                                                               >  1:363879/1‑119 (MQ=255)                                                                                                                                                                                            aTTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATa                                                                                                                                                                                                              <  2:363879/119‑1 (MQ=255)                                                                                                                                                                                                     gTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTc                                                                                               >  2:383267/1‑220 (MQ=255)                                                                                                                                                                                                      tGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCa                                                                                              <  1:383267/220‑1 (MQ=255)                                                                                                                                                                                                         agcaCAGCGACCAACGATTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTACACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGCAAGCTCGGTTAACCGACGCACGGAAAGCTGATTGATCCGGGTCTGGCTCATGGCGACATTGACTCGTGCATTGTCTTCTGGGGAGAGACACCCCCAGAGTGACAg                                       >  1:95368/1‑272 (MQ=255)                                                                                                                                                                                                                   ccaccGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGCTTGATATCAAGATAGGTAGTTGGCTCATCAAGTAATACAACAGGCGTATTCTGGGCCAGGACAATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCTGAAAGCCCGGTTACTCGACGAACGGCAAGAAGATTGATCCTGGTCTGTTTc                                                                                               >  1:14296/1‑206 (MQ=255)                                                                                                                                                                                                                   ccaccGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCgcgctgg                                                                                                                                                              >  1:400942/1‑143 (MQ=255)                                                                                                                                                                                                                   ccaccGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGAtt                              >  1:162394/1‑271 (MQ=255)                                                                                                                                                                                                                    caccGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCgcgctggc                                                                                                                                                             <  2:400942/143‑1 (MQ=255)                                                                                                                                                                                                                        gTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCTGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACCCGTGCATTTTCCTTAGCGGAAAGACGCCCCCCAAGGGAGAGCCAgggggtgtggac                        >  2:421501/1‑263 (MQ=255)                                                                                                                                                                                                                               ccTGAGTCCGGAGTTCGCCCATCAACCGCATCGGGTCCACCGGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATGGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGaa                 <  1:421501/272‑1 (MQ=255)                                                                                                                                                                                                                                    gcacgggtgagcgccAACAACCGTATGTGGTCCACCTGGTGATGGACATCAAGATAGGTGGTTGGTTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑GGCGCGCGGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGAAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAg            <  2:95368/259‑1 (MQ=255)                                                                                                                                                                                                                                       aacagccagCCCATCAACCGCATCTGGTCCACCTGGTGATGGATATGAAGATAGGTGGTGGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGGCGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTc         <  1:424108/264‑1 (MQ=255)                                                                                                                                                                                                                                              cccccATAAACCGCTTCGGGTCCACCTGGTGTTTGTTATAAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACt  <  2:162394/270‑1 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                TTACGGCCTGTGCAATCTACCTCATTAGGCACATCGGCCTGCCAGATACCGGCTCGGGGTGTATTTCCGCTTCCACGCTGAATACTGTTCTCAGCAATCCTGGGGTCATCACCTCTTCTGGTGTGCCTTGCGCCATAACATGTCCGTTTGCCATTACCACCAGTTGATCGCAGTACCGGCTAGCCTGATTAAGGTCGTGCAGCACAGCGACCACCGTTTTCCCCTGAGTCCGGAGTTCGCCCATCAACCGCATCAGGTCCACCTGGTGATTGATATCAAGATAGGTGGTTGGCTCATCAAGTAATACAACGGGCGTATTCTGGGCCAGGACCATTGCCAGAAA‑TGCGCGCTGGCGCTGACCGCCGGAAAGCTCGGTTAACCGACGAACGGCAAGATGATTGATCCGGGTCTGGTTCATGGCGACATTAACTCGTGCATTGTCTTCAGCGGAGAGACGCCCCCAGAGTGACAGCCAGGGATTACGACCATACGAAACCAGCTCCTGGACT  >  NC_000913/4510674‑4511182

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑