Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	481,228	C→A	13.2%	intergenic (‑520/+26)	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,228	0	C	A	13.2%	120.4 / 9.0	38	intergenic (‑520/+26)	tomB/acrB	Hha toxicity attenuator; conjugation‑related protein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (19/14);  new base A (2/3);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCC  >  NC_000913/481085‑481375                                                                                                                                                |                                                                                                                                                    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)                                                                                                                                               aCGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATcc  <  2:383863/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCGCTACTTAATATAAATTAACTAACAATCAGTATTTTTTATGAATTTTATCCGTGGTTAATACTGGTTTTCGTATGAGATCCTGAGTTGGTGGTTCAATTACTCCTTAATGTTCGTAGGTTATGCATAAAAAAGGCCGCTTACGCGGCCTTAGTGATTACACGTTGTATCAATGATGATCGACAGTATGGCTGTGCTCGATATCTTCATTCTTGCGGCTAAAGCGGCGGCGAACCACCACAAAGAATACCGGAACGAAGAAGATTGCCAGTACCGTTGCGGTCACCATCC  >  NC_000913/481085‑481375

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: ‑