Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,359,909	C→T	W41* (TGG→TGA)	cadB ←	lysine:cadaverine antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,359,909	0	C	T	91.7%	55.5 / NA	24	W41* (TGG→TGA)	cadB	lysine:cadaverine antiporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/12);  minor base G (0/2);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTG  >  NC_000913/4359673‑4360154                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      ggAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGGAATACCAATCGCCAGGTTACCAATCCAGTTAGCATTGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCTACATTGCACCAATAATAGAGATAATCCAACCTCAGTTAGCTATacc                                                                                                                                                                                                                                         <  2:86143/251‑1 (MQ=255) ggAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGAAATACCAATCGCCCGGTTACCAATCCAGTTAGCATGGTAATAAAGCACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATacc                                                                                                                                                                                                                                         >  1:436262/1‑251 (MQ=255)         aTCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACCCCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGAGAATTCAACCGCAGATAGTAATACCACCGATAc                                                                                                                                                                                                                                 <  1:326749/251‑1 (MQ=255)            aTTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTg                                                                                                                                                                                                                              >  1:293063/1‑251 (MQ=255)                  tACTGGGAAGAAGGAGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATCAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTACTCTGGCATAAGCCATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACCTACGCCAGCGACATTGCACCTATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGt                                                                                                                                                                                                                        >  1:450817/1‑251 (MQ=255)                                                    gTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCcataat                                                                                                                                                                                      <  1:430336/251‑4 (MQ=255)                                                     aaaTACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGATATAATCCAACCTCAGATATCAATACCACCGATACTTGCTAGTTTCGCAGGTAATAATGCAATACCGCTCCCCATCata                                                                                                                                                                                     <  2:181131/250‑1 (MQ=255)                                                                    ggTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAacaaca                                                                                                                                                                      <  1:321773/251‑1 (MQ=255)                                                                     gTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAacac                                                                                                                                                                     >  2:70518/1‑251 (MQ=255)                                                                              tcAGTTAGCATTGTAATAAACAACACCTATATGAAAGCCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGACCACCTTGTTGCGGGTTTTTTGTTGCCAGTCAGGCATATACATACGCCCGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATCGCAATACCACAGATACTTGCGAGGTTCGCAGGTAATAATGCAATACAGCTCCCCATCATCTTACCGGCAACAACACAGGAACAgg                                                                                                                                                            >  2:57857/2‑251 (MQ=255)                                                                               cagttcgaacggTAGTAAAGAACCCCCGTCTGAAAACCAAATCCAGGGGAAATTTCTCCGGCCTAAGCACTTGGGCCACCTTGTTCCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATCCCGCTCCCCATCATATTACCGGCAACAAAACCGGTACAGGc                                                                                                                                                           >  2:413651/11‑251 (MQ=255)                                                                                 gTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCTCCTTGTTGCGGTTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCaa                                                                                                                                                         <  1:424431/251‑1 (MQ=255)                                                                                  ttAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAGTGCAGGGGAAATTTCTCCGGCCTAAGCAATTGGGCAACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCCATCCCGCTCCCCATCATATTACCGGCCACACAACCGGTACAGGCaa                                                                                                                                                         >  2:426276/1‑250 (MQ=255)                                                                                            aataaAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAAGAGCCCGAtc                                                                                                                                              <  2:305800/251‑1 (MQ=255)                                                                                                                  aCCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGAAAGTCGGGCATATACATACGCCAGCGACATTGAACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACAGCTCCCCATCATATTACCGGCAAAAACACCGGTACAGGCAAATAGCCCGATATTCTTGGCAGACCTCATGctct                                                                                                                        >  1:223600/1‑251 (MQ=255)                                                                                                                          cAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTaa                                                                                                                <  1:250330/251‑1 (MQ=255)                                                                                                                                            cATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTgg                                                                                              >  2:129212/1‑251 (MQ=255)                                                                                                                                                              cTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGgtgt                                                                             >  1:109514/1‑250 (MQ=255)                                                                                                                                                                       ggTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTAttt                                                                   <  2:89837/251‑1 (MQ=255)                                                                                                                                                                                                  catacGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGa                                        <  2:353364/251‑1 (MQ=255)                                                                                                                                                                                                                aTTGCACCAATAATAGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGaaaaa                          <  1:189976/251‑1 (MQ=255)                                                                                                                                                                                                                      ccAATAATAGAGATAATCCAACCGCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTGTGAATTTGGAGTCCGTGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGAtt                    <  2:242732/251‑1 (MQ=255)                                                                                                                                                                                                                            ataGAGATAATCCAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCt              <  1:228638/251‑1 (MQ=255)                                                                                                                                                                                                                                        cAACCTCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTg  <  2:293063/251‑1 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      GGAACAGGATCATTTAATACTGGGAAGAAGGTGGAAAGATAAGATACAGCGGTAATACCAATCGCCAGGTTACCAATCCAGTTAGCATGGTAATAAAGAACACCTGTCTGAAAACCAAATGCAGGGGAAATTTCTCCGGCATAAGCAATTGGGCCACCTTGTTGCGGGTTTTTTGTTGCCAGTCGGGCATATACATACGCCAGCGACATTGCACCAATAATAGAGATAATCCAACCCCAGATAGCAATACCACCGATACTTGCTAGGTTCGCAGGTAATAATGCAATACCGCTCCCCATCATATTACCGGCAACAACACCGGTACAGGCAAATAGCCCGATCTTCTTGGCAGAACTCATGCTCTTCTCCTAATTTCATTTTTGAATTTGGAGTCCGGGTCATGATGTATAACTATTTCCTGACCAGACCAAACTGGCGATAAGATTACTCACGAAAAAAGGATTAATCCTAAAGATTAGGTG  >  NC_000913/4359673‑4360154

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