Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,388,041	T→G	T50P (ACG→CCG)	aaeB ←	aromatic carboxylic acid efflux pump subunit AaeB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,388,041	0	T	G	97.1%	100.6 / NA	34	T50P (ACG→CCG)	aaeB	aromatic carboxylic acid efflux pump subunit AaeB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (11/22);  total (11/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCA  >  NC_000913/3387806‑3388279                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               gcgGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAAt                                                                                                                                                                                                                                 <  1:80680/251‑1 (MQ=255)            ggCCAGCCCCCACGCATACCAGTTTTCTATTCGTACCAGCGATGATATCCAGTTTCAAAAACCGGCCCAGATACAGCACACCAGAATCAGCAATAATGGTGCGCGGATCATCGCAATGATGTTCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCa                                                                                                                                                                                                                      <  1:14486/251‑1 (MQ=255)                aGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTAc                                                                                                                                                                                                                  >  1:395046/1‑251 (MQ=255)                           aTACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGAGAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCGGTCAGTACCGCCCAGCGTg                                                                                                                                                                                                       <  2:87348/251‑1 (MQ=255)                             aCGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAACCCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAAGGATGATCACCAGTCCGTCAATACAGCCAATAAAGGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGc                                                                                                                                                                                                     <  2:142111/251‑1 (MQ=255)                                        aTTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCAAACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCGCCTCCCGCAGCAAATGCCGGACCGGCGGCAACAGTCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAg                                                                                                                                                                                             <  2:85349/248‑1 (MQ=255)                                                             cAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAg                                                                                                                                                                     <  2:104755/251‑1 (MQ=255)                                                              aGGTACAAACACCGGCCCAAATACAGCACACCAGAATCATCAATAATGGAGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAgc                                                                                                                                                                    >  2:138167/1‑251 (MQ=255)                                                                    aaaaaCCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGt                                                                                                                                                              <  1:215274/251‑1 (MQ=255)                                                                        cccGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGAGCATCGCAATGATGATCCCCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTCCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACaa                                                                                                                                                          <  1:174732/250‑1 (MQ=255)                                                                              ccAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCgg                                                                                                                                                    <  1:199035/251‑1 (MQ=255)                                                                                    aCAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGcc                                                                                                                                               <  1:377/250‑1 (MQ=255)                                                                                              aGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCTCCCGGTTTACCGc                                                                                                                                    >  2:13863/1‑251 (MQ=255)                                                                                                      caATAATGGTGCGCGGATCATCGCAATGATGATCAGCAGTCCGGCAATACAGCCAATAAATGTGGCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCGCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGACCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTCCAATGGCGGTCGCCAGTTTTACCGCAAAGCGaa                                                                                                                            <  1:52388/251‑1 (MQ=255)                                                                                                         taatGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCCGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGGatat                                                                                                                         >  1:343912/1‑251 (MQ=255)                                                                                                             ggTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGtt                                                                                                                      <  1:316122/250‑1 (MQ=255)                                                                                                                         atagCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGAGAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGga                                                                                                         <  2:86736/248‑1 (MQ=255)                                                                                                                                       cACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCg                                                                                           <  1:364192/251‑1 (MQ=255)                                                                                                                                        aCCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGt                                                                                          <  1:129546/251‑1 (MQ=255)                                                                                                                                            gTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGTCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCCAAGCGAATATGTTGGTTTGCAATGGGGAAAATACCCATCGTGAtt                                                                                      >  1:73930/1‑251 (MQ=255)                                                                                                                                              ccGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTaa                                                                                    <  1:81752/251‑1 (MQ=255)                                                                                                                                               cGGCAATACAGCCAATAAATGTGCCGATGATGCGCTAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCCAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGAtttac                                                                                   >  1:108686/1‑248 (MQ=255)                                                                                                                                                  cAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCCGTTTTACCGCCAAGCGAATATGTTGGTTATCAATGGAGGAAAAACCCATCGTGATTAACCa                                                                                 >  2:53694/1‑250 (MQ=255)                                                                                                                                                      aCAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGAGTAACCAAACt                                                                             <  1:114018/250‑1 (MQ=255)                                                                                                                                                         gCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCg                                                                                                                                     <  1:169366/191‑1 (MQ=255)                                                                                                                                                           cAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCa                                                                       <  1:53093/251‑1 (MQ=255)                                                                                                                                                                aaTGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGTTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGa                                                                  <  2:60672/251‑1 (MQ=255)                                                                                                                                                                     gCCGATGATTCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCAGCGTGATTAACCAAACTCACGCAGGCGATGGGc                                                             <  2:215669/251‑1 (MQ=255)                                                                                                                                                                                            gATAGCGAATAGCGCCAGAAGACGGTTCACCGCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCTGGCGATGGGCCATTGGACGTAAGAACGAATCCt                                      <  2:268136/251‑1 (MQ=255)                                                                                                                                                                                                          ccAGAATACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGACGCCAGTTTTACCGCAAAGCGAAAATGTTGGTTAGCAATGGAGAAAATACCCATCGTGAGTAAACAAAATCACCCAGGCGGTGGGCCATTTTACGGAAGGACGAATCCTGGCTTTCCTcgcgg                        >  2:290385/1‑250 (MQ=255)                                                                                                                                                                                                            cgtagtcGGTTGACCTCACGCAGGAAAGGCCAGACGGGCGGCAATGATCTGCGCGGTCAGTACCGCCCGGCGTGGAGTTGCCAGCTGGAAGTGAAAGCCAACAAAAATGGCCAGGCAAAGGGCGGGCGCCAGGTTTACCGCAAAGCGAATATGTTGGTTAGCAAGGGAGAGAATACCCATCGTGATTAACCAAACTCACGCAGGCGATTGGCCATTTTACGGAAGGACGAATCCTGGCGTGCGTCGCGATc                      <  2:98891/245‑1 (MQ=255)                                                                                                                                                                                                              aaTACGGTTCACCTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCtt                    >  1:120900/1‑251 (MQ=255)                                                                                                                                                                                                                          cTCCCGCAGCAAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGCAAGTGAAAGCCAACACACAGCACAAGTACAATGGCGGTCGCCCGTTTTACCGCACAGGGAATATGGTGGTTATCAACGGGGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCCATGGGCAATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGGTCTTGTATGCCAggg         >  1:132615/1‑248 (MQ=255)                                                                                                                                                                                                                                cagcaAAGGCCGGACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGaccaccc  >  1:154535/1‑250 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               GCGGTATAACCGGCCAGCCCCCACGCATACGAGTTTTCTATTCGTACCAGCGAGGATATCCAGGTACAAAAACCGGCCCAGATACAGCACACCAGAATCATCAATAATGGTGCGCGGATCATCGCAATGATGATCACCAGTCCGGCAATACAGCCAATAAATGTGCCGATGATGCGCAAAAAGCCACGATAGCGAATAGCGCCAGAATACGGTTCACCTCCCGCAGCAAAGGCCGTACCGGCGGCAACAATCGCCGCTGTCAGTACCGCCCAGCGTGGCGTTTCCAGCTGGAAGTGAAAGCCAACAAACAGCGCCAGTACAATGGCGGTCGCCAGTTTTACCGCAAAGCGAATATGTTGGTTAGCAATGGAGAAAATACCCATCGTGATTAACCAAACTCACGCAGGCGATGGGCCATTTTACGGAAGAACGAATCCTGGCTTTCGTCGCGATCTTGTTTGCCAGTGACCACCA  >  NC_000913/3387806‑3388279

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