Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,388,726 A→G V134A (GTT→GCT)  aaeA ← aromatic carboxylic acid efflux pump membrane fusion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,388,7260AG98.2% 191.8 / ‑5.0 55V134A (GTT→GCT) aaeAaromatic carboxylic acid efflux pump membrane fusion protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (22/32);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCCTCGAGTAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAACCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAGCAACATCGGCTTGCGCT  >  NC_000913/3388590‑3388860
                                                                                                                                        |                                                                                                                                      
tCCTCGAGTAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCg                                                                                                                                     <  2:554702/140‑1 (MQ=255)
tCCTCGAGTAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCg                                                                                                                                     <  2:663890/140‑1 (MQ=255)
 ccTCGAGTAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGt                                                                                                                                    <  2:852053/140‑1 (MQ=255)
   tCGAGTAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGttt                                                                                                                                  >  1:13678/1‑140 (MQ=255)
     gAGTAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgt                                                                                                                                >  2:814910/1‑140 (MQ=255)
       gTAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtag                                                                                                                              <  1:601469/140‑1 (MQ=255)
         aataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                                            <  2:13678/140‑1 (MQ=255)
         aataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                                            <  1:837484/140‑1 (MQ=255)
         aataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                                            <  1:79217/140‑1 (MQ=255)
         aataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTgtagta                                                                                                                            <  1:200299/140‑1 (MQ=255)
          ataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTAc                                                                                                                           >  1:818832/1‑140 (MQ=255)
          ataaACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTAc                                                                                                                           >  2:646486/1‑140 (MQ=255)
                    cGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGcc                                                                                                                 <  1:547828/140‑1 (MQ=255)
                    cGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGcc                                                                                                                 <  1:77834/140‑1 (MQ=255)
                                 gAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCtt                                                                                                    >  2:594786/1‑140 (MQ=255)
                                 gAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCtt                                                                                                    >  1:590107/1‑140 (MQ=255)
                                 gAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCtt                                                                                                    >  1:414233/1‑140 (MQ=255)
                                       ggTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCgag                                                                                              >  1:761402/1‑140 (MQ=255)
                                         tCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGAc                                                                                            <  2:467772/140‑1 (MQ=255)
                                         tCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGAc                                                                                            <  1:814910/140‑1 (MQ=255)
                                         tCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACa                                                                                           >  1:744579/1‑140 (MQ=255)
                                           aCCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTTGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACAt                                                                                          <  2:761402/140‑1 (MQ=255)
                                           aCCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACAt                                                                                          <  2:679696/140‑1 (MQ=255)
                                           aCCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACAt                                                                                          <  1:594786/140‑1 (MQ=255)
                                            cccAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGTCCTGGTCg                                                                                                           >  2:237384/1‑122 (MQ=255)
                                            cccAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGTCCTGGTCg                                                                                                           <  1:237384/122‑1 (MQ=255)
                                                  cATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACAt                                                                                          >  2:121038/1‑133 (MQ=255)
                                                  cATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACAt                                                                                          <  1:121038/133‑1 (MQ=255)
                                                        ctgGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCg                                                                             >  2:6925/1‑140 (MQ=255)
                                                         tgGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCga                                                                            >  1:801535/1‑140 (MQ=255)
                                                          ggCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGca                                                                                  <  1:234086/133‑1 (MQ=255)
                                                          ggCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGca                                                                                  >  2:234086/1‑133 (MQ=255)
                                                              cgcgGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACgg                                                                       >  2:69985/1‑140 (MQ=255)
                                                                cgGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGTCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGtt                                                                     >  1:832798/1‑140 (MQ=255)
                                                                cgGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGtt                                                                     >  2:637134/1‑140 (MQ=255)
                                                                           gCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCgg                                                          >  2:174117/1‑140 (MQ=255)
                                                                             gTTCAAGATCCAGTTTTGCCAGATCGCGTGTCGCCTGCGCTTTCGCTAACTGATGCAGAACCGTTTGTAGTACGTTGTTGACCTGGGCGATCTCTTCGCGAGACATCGCTTGGACACTGAGACAGTCACGCCGTCCGGcc                                                        >  2:618649/1‑140 (MQ=255)
                                                                               tCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGcctc                                                      <  1:782134/140‑1 (MQ=255)
                                                                                        aGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGttt                                             >  2:174202/1‑140 (MQ=255)
                                                                                         gTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTtc                                            <  2:829040/140‑1 (MQ=255)
                                                                                         gTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTtc                                            >  2:292038/1‑140 (MQ=255)
                                                                                            ttGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTcc                                         <  1:358570/140‑1 (MQ=255)
                                                                                                 aGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGc                                    <  2:818832/140‑1 (MQ=255)
                                                                                                 aGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGc                                    <  2:801535/140‑1 (MQ=255)
                                                                                                   aTCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCa                                  <  1:65090/140‑1 (MQ=255)
                                                                                                    tCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTcc                                         >  2:859077/1‑132 (MQ=255)
                                                                                                    tCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTcc                                         <  1:859077/132‑1 (MQ=255)
                                                                                                       gcgggcGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTAc                              <  2:744579/135‑1 (MQ=255)
                                                                                                        gggTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCTTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTAcc                             <  1:699578/140‑1 (MQ=255)
                                                                                                        gggTCGCCTGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTAcc                             <  1:69985/140‑1 (MQ=255)
                                                                                                                tGCGCTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGataata                     <  2:147140/140‑1 (MQ=255)
                                                                                                                  cgcTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAg                   <  1:637134/140‑1 (MQ=255)
                                                                                                                  cgcTTTCGCTAACTGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAg                   <  1:174202/140‑1 (MQ=255)
                                                                                                                               tGATGCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAGCAACATCGGCTTg      >  2:155311/1‑140 (MQ=255)
                                                                                                                                   gCAGAGCCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAGCAACATCGGCTTGCGCt  <  2:138497/140‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TCCTCGAGTAATAAACTCACCGGTATAGACGTTGAGGTTGGTCACCCAGCCATCTGCTGGCGCGCGGATCACCGTGCGTTCAAGATCCAGTTTTGCCAGATCGCGGGTCGCCTGCGCTTTCGCTAACTGATGCAGAACCGTTTGTAGTACGTTGTTGGCCTGGTCGATCTCTTCGCGAGACATCGCCTGCACACCGAGACGGTTACGACGTCCGGCCTCCTGGCGTTTCTCCTGTGCCAGTACCTGATAATAAGCAACATCGGCTTGCGCT  >  NC_000913/3388590‑3388860

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