Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,567,272 C→A intergenic (‑26/+15) yjiN ← / ← mdtM DUF445 domain‑containing protein YjiN/multidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,567,2720CA94.7% 46.9 / ‑4.2 21intergenic (‑26/+15)yjiN/mdtMDUF445 domain‑containing protein YjiN/multidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (10/8);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

CGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATT  >  NC_000913/4567025‑4567472
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                        
cGTATGGCGGGAAATGATCGGAACCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTGCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGa                                                                                                                                                                                                       <  2:169605/251‑1 (MQ=255)
 gTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGAt                                                                                                                                                                                                      <  1:219569/251‑1 (MQ=255)
                aTCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACt                                                                                                                                                                                       <  1:259118/251‑1 (MQ=255)
                           cACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAAGTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACt                                                                                                                                                                            <  1:236763/251‑1 (MQ=255)
                                        aGCGCCACCACCGCAAACCAGCCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCGAAAATCACTGCTCCTCCACTGGCGCGGCTGCCt                                                                                                                                                               <  1:86117/251‑1 (MQ=255)
                                                        cccAGTCCGCCATCGCGCCAACCACCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCATCCTTTGGGCGCGTCTGAGTTCGAGGAGTTTATTCAGAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTATCTCGGCTGCATGATGCTGGCGTAcgcg                                                                                                                                               <  2:247842/250‑1 (MQ=255)
                                                          cAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAAGCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTGGGCGCGTCTGAGTTCGATGAGGTTATTCATAACGTTAAGAATAGAGCCTGACGGCTGGATGGCGCAAAAATCTCTGCTCCTCCACTATCTCGGCTGCCTGATGCTGGCGCACGCGAt                                                                                                                                             <  2:263397/251‑1 (MQ=255)
                                                                      gcgcCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGAGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCACAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCg                                                                                                                                 >  2:140064/1‑251 (MQ=255)
                                                                        gcctaCCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGACCGCGCAAAAATCACTGCTCCTCCACTAGCTCTGCTGCCTGATGCTGGCGCACGCGCTTGAGCAATCCCGcc                                                                                                                               >  2:80187/5‑251 (MQ=255)
                                                                          cGACCATCGCCGCTCCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAg                                                                                                                             <  1:191792/251‑1 (MQ=255)
                                                                                                cGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAAcgccc                                                                                                       >  1:248076/1‑251 (MQ=255)
                                                                                                  ccTTCACGCCGCTCACACAAAAATTGGGCGGCTAAAAACGCGTAACGACACAGGTATCAGCGGCGCTAAGCAGTAAAGATAGCTCCACCTTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTGACTCGGCTGTCTGACGCTGCCGCCCGCCATTGCGCACTCTCGCCAGGGTGATCACGACAGTAAcgccccc                                                                                                     >  2:198904/1‑249 (MQ=255)
                                                                                                              tCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGACAGGGTGAAAACGACGATAACGCCCGCCCCAACGGCTa                                                                                          >  1:92641/1‑250 (MQ=255)
                                                                                                                            gggcgcAAAAACAACCGTACCGACAAAGGTAGCAGCGGCATAAATAATAAAAGAGAGCGCAACCCTTTGGGCGCTTCTGATTCCGATGTGTTTATTAATAAGTTTAAAAATAGGGCCGGCCGGCGGAATCGCACAAAAATCCCTGCTCCGCCACAAGCTCGTCTGCCTAATGCTGGCGCCCGTGATTGAGCAATCCCGCGAGGGTAAAAACGACGAAACCGCCCGCCACAACGGCTAACAGATGAAACGGc                                                                            <  1:252036/247‑1 (MQ=255)
                                                                                                                               ggcAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGc                                                                        >  1:11575/1‑251 (MQ=255)
                                                                                                                                 cAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCATTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCg                                                                       >  1:112364/1‑250 (MQ=255)
                                                                                                                                     aaaCAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCGATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGAccg                                                                   >  2:151886/1‑251 (MQ=255)
                                                                                                                                                    cAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAgg                                                                                                                            <  1:209856/180‑1 (MQ=255)
                                                                                                                                                    cAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAgg                                                                                                                            >  2:209856/1‑180 (MQ=255)
                                                                                                                                                                                                       cgTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCAGCCACAACGGCGAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGAGGCCGATTTCGACCGAGACCGTCATCAACCTCAGCAACACCatat   >  1:150651/1‑251 (MQ=255)
                                                                                                                                                                                                        gTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCAGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATAtt  >  1:147570/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                        
CGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAGCGCCAACCTTTTGGCGCGTCTGAGTTCGATGAGTTTATTCATAACGTTAAGAATAGAGCCTGACGGCCGGATCGCGCAAAAATCACTGCTCCTCCACTAGCTCGGCTGCCTGATGCTGGCGCACGCGATTGAGCAATCCCGCCAGGGTGAAAACGACGATAACGCCCGCCACAACGGCTAACAGATGAAACGGCAAGCGACCGCCGTTAAACCATAGCCAGCGGCCGATTTCGACCGAGACCGACATCACCATCAGGATCACCATATT  >  NC_000913/4567025‑4567472

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