Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_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_pae_tpiA482,2300GA96.8% 85.6 / NA 31P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (14/16);  total (14/17)
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

TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTT  >  NC_000913_3_pae_tpiA/482084‑482373
                                                                                                                                                  |                                                                                                                                               
tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                               >  2:552913/1‑149 (MQ=255)
tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                               >  2:546792/1‑149 (MQ=255)
tgacatAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                               <  2:458023/149‑1 (MQ=255)
tgaattgaaCTTTCTGAACACGGCCGCGGTCGAGAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGGCGTTGATAGTAACACCCGGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCagg                                                                                                                                               <  1:217756/142‑1 (MQ=255)
     tAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGACTGGTGGTAATGTCGTTGATCGAACCCCCCAGCTCCTGCCCTTTTTCCTGCTCGATACCACTCCAACACTGCAGGCTATc                                                                                                                                          >  2:277588/1‑149 (MQ=255)
                  cacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                             <  2:380489/149‑1 (MQ=255)
                     gACCGCGGTCGATAAAGTCGTTCACATAGCGGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGa                                                                                                                          <  2:258953/149‑1 (MQ=255)
                        cgcgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCATCGTCGACATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGt                                                                                                                       >  2:229540/1‑149 (MQ=255)
                        cgcgGTCGATAAAGTAGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGCGTGGTGCTACTGTCGGTGCTAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATCTCAACCTTACACTACAGCGGCTCTTCCGGACCGTTTgcccga                                                                                                                       >  2:116750/1‑144 (MQ=255)
                                aTAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGt                                                                                                               <  2:635491/149‑1 (MQ=255)
                                 tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                              <  1:546792/149‑1 (MQ=255)
                                 tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTACCTCTTTGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                              <  1:344644/149‑1 (MQ=255)
                                   aaGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCaa                                                                                                            >  2:468902/1‑149 (MQ=255)
                                    aGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAAc                                                                                                           <  2:282111/149‑1 (MQ=255)
                                            acaTAGCTGCCGCCGCATGCCGCTCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCCTTTCCCCGGGCGATATCAATCTTAAACTGCGGGGAATAGTCGAGACCGTTTGGATGTACGCTGGTCAACATATCAgg                                                                                                   <  2:272466/149‑1 (MQ=255)
                                            acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCTCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTTCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                   <  2:522346/149‑1 (MQ=255)
                                            acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                   <  1:552913/149‑1 (MQ=255)
                                              atatCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATCTCCGGGt                                                                                                 >  2:409384/1‑149 (MQ=255)
                                                              gCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgc                                                                                 <  2:266866/149‑1 (MQ=255)
                                                                     ccAGAGTGGTGTTAATGTCGTTGATCGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAAc                                                                          >  1:318361/1‑149 (MQ=255)
                                                                     ccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGAGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCATCTGGCAGTAGCCTGATCAACAAAACATGGTGCTTCGCTCCTTCTGCAAGCAAc                                                                          >  2:252018/1‑149 (MQ=255)
                                                                            ggtgTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAc                                                                   <  1:335944/149‑1 (MQ=255)
                                                                                      tCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAAGCGCCTGAGt                                                         >  2:452471/1‑149 (MQ=255)
                                                                                      gcggtGATAGAAACACCCAGCGCCGGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATGTCAGGGTGCTTCGCTGCTGCTGCAAGCAACTGGTTACGCGCCTGAGt                                                         <  2:119916/145‑1 (MQ=255)
                                                                                               gAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                <  1:136354/149‑1 (MQ=255)
                                                                                                aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTAAGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                               >  1:594353/1‑149 (MQ=255)
                                                                                                  acacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGCGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTg                                             >  2:327611/1‑149 (MQ=255)
                                                                                                       cAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCCGTTTTTCGTGCCCaa                                        >  2:96648/1‑149 (MQ=255)
                                                                                                               gcgcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGcc                                >  2:635387/1‑149 (MQ=255)
                                                                                                                   tttttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGt                            <  1:318935/149‑1 (MQ=255)
                                                                                                                           ggTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACAGGTTACGCGCCTGCTTCAGTTTTTCGTAACCAAGCCCAGCCTGGTCAACCacc                    >  2:230277/1‑147 (MQ=255)
                                                                                                                                atatCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa               <  1:271541/149‑1 (MQ=255)
                                                                                                                                  atCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa               <  1:58275/147‑1 (MQ=255)
                                                                                                                                  atCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaa               >  2:58275/1‑147 (MQ=255)
                                                                                                                                             aCTGCAGGGTATCTACCAGACCGTTTGGACGTACGCTGGTCAAAATATAAGCTCTCGTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGtt  <  1:409384/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGACATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTT  >  NC_000913_3_pae_tpiA/482084‑482373

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