Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	481,225	T→G	13.5%	intergenic (‑517/+29)	tomB ← / ← acrB	Hha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	481,225	0	T	G	13.5%	95.2 / 11.7	37	intergenic (‑517/+29)	tomB/acrB	Hha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (19/13);  new base G (2/3);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

TCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGC  >  NC_000913/481078‑481364                                                                                                                                                    |                                                                                                                                            tCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCt                                                                                                                                             >  3:219314/1‑148 (MQ=255)        gcgaTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCgg                                                                                                                                      <  3:363936/148‑1 (MQ=255)          gcTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCt                                                                                                                                   <  1:179140/149‑1 (MQ=255)           cTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCGTAAGCGGCCtt                                                                                                                                  <  1:156468/149‑1 (MQ=255)               ttAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTg                                                                                                                              >  1:253726/1‑149 (MQ=255)                           aactaacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGt                                                                                                                  <  2:103817/149‑1 (MQ=255)                              taacAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATc                                                                                                               >  2:272923/1‑149 (MQ=255)                                   aTCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatga                                                                                                          >  2:216025/1‑149 (MQ=255)                                    tCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAatgat                                                                                                         <  1:272923/149‑1 (MQ=255)                                      aGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAAtgatga                                                                                                       >  2:374323/1‑149 (MQ=255)                                        tATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGa                                                                                                                             <  1:314908/125‑1 (MQ=255)                                        tATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGa                                                                                                                             >  2:314908/1‑125 (MQ=255)                                        tATTTTTTACGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATc                                                                                                     >  1:3511/1‑149 (MQ=255)                                          ttttttATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGa                                                                                                   <  1:29005/149‑1 (MQ=255)                                          ttttttACGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGa                                                                                                   <  2:3511/149‑1 (MQ=255)                                                   aaTTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                                                          >  2:229038/1‑149 (MQ=255)                                                   aaTTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGc                                                                                          <  2:253726/149‑1 (MQ=255)                                                      tttATCCGTGGTTAATACTGGTTTTAGTATGAGATCCTGAGTTGGTGGTTAAATTACTCCTTAATGTACGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCtg                                                                                        <  4:18628/148‑1 (MQ=255)                                                            cGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCg                                                                                  <  3:66709/148‑1 (MQ=255)                                                                 ttAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCt                                                                            >  1:383863/1‑149 (MQ=255)                                                                 ttAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCt                                                                            <  1:216025/149‑1 (MQ=255)                                                                      aCTGGATTTCGTATGAGATCCTGAGTTGGTGGTACAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCAtt                                                                       >  1:116575/1‑149 (MQ=255)                                                                                gTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTa                                                             >  2:262722/1‑149 (MQ=255)                                                                                  aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAg                                                                                                                                >  3:221453/1‑80 (MQ=255)                                                                                  aTGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAg                                                                                                                                <  4:221453/80‑1 (MQ=255)                                                                                      gaTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCtt                                                                    <  3:393547/136‑1 (MQ=255)                                                                                      gaTCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCtt                                                                    >  4:393547/1‑136 (MQ=255)                                                                                                    ggtTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACcaccac                                         <  2:116575/149‑1 (MQ=255)                                                                                                               tCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGACCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCg                              >  1:250416/1‑149 (MQ=255)                                                                                                                ccTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCg                              >  3:371551/1‑148 (MQ=255)                                                                                                                  ttAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGa                            >  3:257128/1‑148 (MQ=255)                                                                                                                         tCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACgaaga                     <  3:393951/148‑1 (MQ=255)                                                                                                                             aGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGAtt                >  1:173132/1‑149 (MQ=255)                                                                                                                                        aaaaaaGGCCTCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGTCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg      >  4:115203/1‑148 (MQ=255)                                                                                                                                        aaaaaaGGCCGCTTACGCGGCCTTAGTGATTACACATTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCg      >  4:108056/1‑148 (MQ=255)                                                                                                                                          aaaaGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTg   >  1:202751/1‑149 (MQ=255)                                                                                                                                           taaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc  >  1:337713/2‑149 (MQ=255)                                                                                                                                            aaGGCCGCGTAAGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGc  >  3:98923/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                            TCTATTAGCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGC  >  NC_000913/481078‑481364

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