Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_vch_tpiA 482,230 G→A P725L (CCG→CTG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA482,2300GA97.2% 84.4 / NA 37P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (21/14);  minor base C (1/0);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAG  >  NC_000913_3_vch_tpiA/482090‑482367
                                                                                                                                            |                                                                                                                                         
aaaCTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCt                                                                                                                                   >  1:558323/1‑149 (MQ=255)
 aaCTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCtt                                                                                                                                  >  2:386199/1‑149 (MQ=255)
 aaCTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCtt                                                                                                                                  >  2:384047/1‑149 (MQ=255)
      tcttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGAGCGTCTTCCAGa                                                                                                                             >  1:407110/1‑149 (MQ=255)
             acGACCTCGGTCGATAAAGGCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGGGGTGTTAATGTCGTTGGTAGAAACACCCAGCGCCGGCGGTTTTTACTGGTTGAGATCAATCTTAAACTGCAGGGTATCTGCCAGAGCGGTTg                                                                                                                      >  2:19658/1‑149 (MQ=255)
             acGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTATTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTg                                                                                                                      >  1:92469/1‑149 (MQ=255)
                gtcgcgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCCTGCAGCGCCCATAGTGGTGTTAATGTCTTTGATAGAATCACCCAGCGTCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGAc                                                                                                                   <  1:214137/147‑1 (MQ=255)
                           tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCACTCTTAAACTCCCGGGTATCTTCCAGACCGCTTGCACGTACGCTGGTc                                                                                                        >  2:575742/1‑149 (MQ=255)
                                tCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACTCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGacaacaa                                                                                                   >  2:323751/1‑148 (MQ=255)
                                  gTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                 >  2:88629/1‑149 (MQ=255)
                                  ctttataTAGCGGCCGCCCCATGCATCTCCCAGAGTAGTGTTAATGTCGTTGATCGAAACACCCAGCTGCTGCGCTTTTTTCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACatat                                                                                                 <  2:493645/143‑1 (MQ=255)
                                      acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGTAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                             <  1:325992/149‑1 (MQ=255)
                                        aTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACa                                                                                                    >  1:30476/1‑140 (MQ=255)
                                        aTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACa                                                                                                    <  2:30476/140‑1 (MQ=255)
                                         tAGCTACCGCGCCCTACAGCAACAAGAGTGGGGTTAATGGCGTTGATAGAAACACCCAGCGCCTGCGCTGTTTCCTGGTCGCTATAAACCTTACACTGCAGGGCATCGTACAGCGCGTTTGGACGTACGCTGGTCAAAATATCAGGGTg                                                                                          >  2:577660/1‑149 (MQ=255)
                                            cTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGAATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCtt                                                                                       >  1:320496/1‑149 (MQ=255)
                                                          aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                         <  1:384047/149‑1 (MQ=255)
                                                          aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                         <  1:386199/149‑1 (MQ=255)
                                                          aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                         <  1:466669/149‑1 (MQ=255)
                                                               ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                    >  2:171436/1‑149 (MQ=255)
                                                                       gtgtTAATGTCGTTCATAGAAACACCCAGCCCCAGCGCTTTTTCCTGGTCGATATCAACCTTAAACTGCAGGGTATCTTCCAGATCGTTTAGACGTACGTTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcg                                                            <  1:323751/149‑1 (MQ=255)
                                                                        tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTTGCTGCTTCTGCGAGCAACTGGTTAcgc                                                           >  1:517748/1‑149 (MQ=255)
                                                                                  gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTACAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                 >  1:528611/1‑149 (MQ=255)
                                                                                                  aGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAg                                 <  1:88629/149‑1 (MQ=255)
                                                                                                       cTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTAGACGTACGCTAGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAATGGTTACGCGCCTGAGTCAGTTTTTAGTGACCAAGGCCAg                            >  1:50930/1‑149 (MQ=255)
                                                                                                        tGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGc                           >  1:385455/1‑149 (MQ=255)
                                                                                                        tGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGc                           <  1:150095/149‑1 (MQ=255)
                                                                                                           gcTTTTTCCTTGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                        <  2:77817/149‑1 (MQ=255)
                                                                                                             tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAAACCGTTTGAACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCCGCCTGGt                      >  1:167246/1‑149 (MQ=255)
                                                                                                                         gATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCCGTTTTTCGTGACCAAGGCCAGCCCGGTCAATCAGCtcat          >  2:426879/1‑148 (MQ=255)
                                                                                                                         gATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaa          <  2:136607/149‑1 (MQ=255)
                                                                                                                          atatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa         >  1:156889/1‑149 (MQ=255)
                                                                                                                           tatCAATCTTAACCTGCAGGGTATCTTGCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaag        <  2:92469/149‑1 (MQ=255)
                                                                                                                             tCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACAAAGGCCAGCCTGGTCAATCAGCtcaaagtc      >  1:90327/1‑149 (MQ=255)
                                                                                                                              cAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGAACAAAATAACAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtca     >  1:15265/1‑149 (MQ=255)
                                                                                                                                aTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTTAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaaa   <  1:171436/149‑1 (MQ=255)
                                                                                                                                 tCTTAAACTGCAGGGTATCTTCCAGACAGTAAGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCTCCTGAGTCAGTTTTTCGTGACCAAGTGCAGCCTGGACAATCAGCgtcaagtcaaag  >  2:31683/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                         
AAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAG  >  NC_000913_3_vch_tpiA/482090‑482367

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