Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,744,751 Δ1 bp coding (1295/1374 nt) mdtK → multidrug efflux pump MdtK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,744,7510C.95.7% 85.6 / ‑3.3 23coding (1295/1374 nt)mdtKmultidrug efflux pump MdtK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (12/10);  minor base A (1/0);  total (13/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 = 9.39e-01

TGATGTTGCTGGCGGCGGTATATCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGCAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTATCTCAACGTGTATATCCCCGGTAACTTCCCCT  >  NC_000913/1744526‑1744991
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
tGATGTTGCTGGCGGCGGTATATCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGc                                                                                                                                                                                                                        >  2:279078/1‑251 (MQ=255)
tGATGTTGCTGGCGGCGGTATATCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCGTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTTCGTATg                                                                                                                                                                                                                         <  2:264329/250‑1 (MQ=255)
              gcggTATATCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGt                                                                                                                                                                                                          >  1:373470/1‑251 (MQ=255)
                   atatCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGcc                                                                                                                                                                                                     >  1:205710/1‑251 (MQ=255)
                                  cTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTAt                                                                                                                                                                                                                          >  2:245088/1‑215 (MQ=255)
                                  cTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTAt                                                                                                                                                                                                                          <  1:245088/215‑1 (MQ=255)
                                        ccAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGTGCTTGCCAAGCGGCTATAGTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGGGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACg                                                                                                                                                                                <  1:183559/251‑1 (MQ=255)
                                           ggTGATTGGCAGTGGGATTTTGCGCGGTTATACAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGc                                                                                                                                                                             >  2:383856/1‑251 (MQ=255)
                                              gATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATc                                                                                                                                                                          <  2:8778/251‑1 (MQ=255)
                                                      gTGGGATCTTGCGTGGTTATAACGATCAGCGTTCCCTTTTCTATATTACCTTTACGGCATACTGGCTGCTGGGATTGCCCAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGACCCTATGGGGCAAGCAGGCTTCTGGCTAGGCATTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCATATGCGGTTCCTGCAACGTATGCCGTCAGCCATCATTCTGCAAAGAGCATCCCGAtaat                                                                                                                                                                  >  2:366788/1‑251 (MQ=255)
                                                           aTTTTGCGTGGTTATAAAGATAAGAGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGAACTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCAGTCAGCCATAATTCTGCAACGAGCATCCCGCtactaa                                                                                                                                                                >  2:192957/1‑248 (MQ=255)
                                                            ttttGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGAc                                                                                                                                                            >  2:371068/1‑251 (MQ=255)
                                                                         taAAGATACGCGTTCCATTTTCTATATTACCTTTGCGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCtt                                                                                                                                               >  1:379092/1‑251 (MQ=255)
                                                                                   tgTTCCATTATGTATCTTACCTTTACGGCTTACTGGGTGCCGGGCTTGCCATGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTTGCCTGACGTCGGAGCCATTTTGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCCGCCATCCATCTGCAACGAGCATCCCGCTAATACAGACAATGCGCTACCTTCACGGGtagc                                                                                                                                     >  2:58902/2‑248 (MQ=255)
                                                                                                      ccTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGAACTGACCGATCTGGACGTTGAACATATGGGGCCAGCAGGCTTCTGGATAGGCGTTATTATTGGCAGGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGGTCCTGCAACGTCTGCCGTTAGCCATCATTCTGCAACGAGCAACCCGCTAATAAAGACAAGGCGCAACCTTCAAGGGTTGCGACGGTAGTTTTACGCAgg                                                                                                                  >  1:391278/1‑251 (MQ=255)
                                                                                                                          cTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCcatcat                                                                                                                                                                                         <  2:354054/160‑1 (MQ=255)
                                                                                                                          cTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCcatcat                                                                                                                                                                                         >  1:354054/1‑160 (MQ=255)
                                                                                                                                                                   ggTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACag                                                         <  1:221462/247‑1 (MQ=255)
                                                                                                                                                                                     gCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGtt                                   <  1:405785/251‑1 (MQ=255)
                                                                                                                                                                                          cagGCTTATGGATAGGATTTATTACTCCCCTGACGTCGAAACCATTATGATGATGCTGCGTATGCGATTCCTGCAACGTCTGCCGTCAACCATCATTCTGAAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTAAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTTtctc                              >  2:48496/1‑251 (MQ=255)
                                                                                                                                                                                                   ggATAGGCTTTATTATTGGCCTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTATCTCAACgtgt                       <  2:63669/249‑1 (MQ=255)
                                                                                                                                                                                                                   gggCCTGACGTCGGAGCCATTATGATGATGTTGCGTATTCGGTTCCTGCAACGTCTGCCTGCAGCCATCATTCTGGAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCGAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGTGGTTTCCGACTGGCAGCTGTGATCTCAACGTGGATATCCCCGGTAGCTTcc     <  2:176056/250‑1 (MQ=255)
                                                                                                                                                                                                                      ccTGACGTCGGAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACGGAGTTTGCCGACTGGCATCTGTTATCTCAACGTGTATATCCCCGGTAACTTCCCCt  <  2:379092/251‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
TGATGTTGCTGGCGGCGGTATATCAGATTTCTGACTCAATCCAGGTGATTGGCAGTGGGATTTTGCGTGGTTATAAAGATACGCGTTCCATTTTCTATATTACCTTTACGGCTTACTGGGTGCTGGGCTTGCCAAGCGGCTATATTCTGGCACTGACCGATCTGGTCGTTGAACCTATGGGGCCAGCAGGCTTCTGGATAGGCTTTATTATTGGCCTGACGTCGGCAGCCATTATGATGATGTTGCGTATGCGGTTCCTGCAACGTCTGCCGTCAGCCATCATTCTGCAACGAGCATCCCGCTAATAAAGACAAGGCGCAACCTTCACGGGTTGCGCCTGTATTTTTACGCAGGCTGGAGCGTTGCGCCAATCCCGTCTTCGTCTGGCTGTAATTTCAGAGCGTTACACAGAGTTTGCCGACTGGCATCTGTTATCTCAACGTGTATATCCCCGGTAACTTCCCCT  >  NC_000913/1744526‑1744991

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