Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 485,272 C→T Q116Q (CAG→CAA acrA ← multidrug efflux pump membrane fusion lipoprotein AcrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913485,2720CT93.8% 38.2 / NA 16Q116Q (CAG→CAAacrAmultidrug efflux pump membrane fusion lipoprotein AcrA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/6);  minor base G (0/1);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTAC  >  NC_000913/485058‑485493
                                                                                                                                                                                                                      |                                                                                                                                                                                                                             
tCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTGTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTGCTCTTGCGTACTGATGTACGGAGTACCGAGCAGTTTCTTATAACGATTCACCGTCAATTGCGCTATATTGGCTGCATCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTa                                                                                                                                                                                            <  1:194576/250‑1 (MQ=255)
                       aTCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATc                                                                                                                                                                    <  1:148434/251‑1 (MQ=255)
                             gagGTGACTTTGGTGGAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCCGGTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCGTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTg                                                                                                                                                              <  1:82407/251‑1 (MQ=255)
                                         gtgtAAGACAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAgagagaga                                                                                                                                                   >  1:171009/1‑250 (MQ=255)
                                                   gATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGGTACCGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCTAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGGTGCAGCTTGGGCTTTCGCCAGAGCACCTTTCGCACTATCGTATGTCGCCTGAAAGGTCGCAGGATCAATCTGATAGAGAGAGACTCCTGCTTCg                                                                                                                                        <  1:158634/251‑1 (MQ=255)
                                                       gATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGt                                                                                                                                    >  1:12310/1‑251 (MQ=255)
                                                                                    ccgcAGTTACCGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCt                                                                                                                                           >  1:190500/1‑215 (MQ=255)
                                                                                             ccGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGAGCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCg                                                                                              <  1:176316/251‑1 (MQ=255)
                                                                                                               gcgcATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGc                                                                                                                                                                                                                    <  1:199726/115‑1 (MQ=255)
                                                                                                                 gcATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTggg                                                                                                                                                                                                                           <  1:50844/106‑1 (MQ=255)
                                                                                                                        ccAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTCCCGAGCAGTTTCTGATAACGATTCACCGTCACTTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGAACGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGCACTTCTGCGATc                                                                   >  1:27636/1‑251 (MQ=255)
                                                                                                                             gCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGt                                                               >  1:177902/1‑250 (MQ=255)
                                                                                                                                        ctctTGCTTACTGATGTACTGAGTACCGAGCAGTTTCAGATAACGATTCCCAGTCAATTGCGCGAAATTGGCCGCCGCTTGGGATTTCGCCAGATCACCTTACGCCCTGTCGTATGACGCCTGATAGGACGAAGGCTCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTAAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGc                                                   >  1:144480/1‑251 (MQ=255)
                                                                                                                                                            gAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGANACANNNNNtt                                                                                                                                          >  2:128237/1‑144 (MQ=255)
                                                                                                                                                                             tGATAACGATTAACCGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCAGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCt                           >  1:72233/1‑238 (MQ=255)
                                                                                                                                                                                         nccGTCAATTGCGCGATATTGGCTGCAGCTTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTtac  >  1:19563/2‑251 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                             
TCGACTTACCAATGCGACCGCTAATCGGAGAGGTGACTTTGGTGTAAGCCAGATTGATCCGCGCAGTTTCAACGGCAGCTTTCGCCGCAGTTACCGCAGCATTCGCCTGTTGCGCATCAGCCAGAGCCTGATCGTACTCTTGCTTACTGATGTACTGAGTACCGAGCAGTTTCTGATAACGATTCACCGTCAATTGCGCGATATTGGCTGCAGCCTGGGCTTTCGCCAGATCACCTTTCGCACTGTCGTATGTCGCCTGATAGGTCGCAGGATCAATCTGATAGAGAGAGACACCTGCTTCGATGTCGCTACCTTCTTTGAAATTACGCTTCAGGATAATCCCGCTAACTTGAGGACGAACTTCTGCGATCCGGTAGGCACTGGTGCGACCCGGAAGCTCGGTTGTGATCTGCAGAGGTTCAGTTTTGACTGTTAC  >  NC_000913/485058‑485493

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