Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,655,767 T→C T198A (ACT→GCT)  pepB ← aminopeptidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,655,7670TC98.0% 172.5 / NA 49T198A (ACT→GCT) pepBaminopeptidase B
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (22/26);  total (22/27)
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

AGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCTCA  >  NC_000913/2655631‑2655900
                                                                                                                                        |                                                                                                                                     
aGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGa                                                                                                                                     <  2:814880/139‑1 (MQ=255)
aGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGa                                                                                                                                     <  1:882552/139‑1 (MQ=255)
 gCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGAt                                                                                                                                    >  1:434822/1‑139 (MQ=255)
 gCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGAt                                                                                                                                    >  2:729129/1‑139 (MQ=255)
   ttCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATAc                                                                                                                                  <  2:232006/139‑1 (MQ=255)
      aCGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCtt                                                                                                                               <  1:846868/139‑1 (MQ=255)
        gcgcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                             <  2:621757/139‑1 (MQ=255)
        gcgcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTta                                                                                                                             <  2:308311/139‑1 (MQ=255)
          gcTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtacc                                                                                                                           <  2:1020687/139‑1 (MQ=255)
            ttGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTtaccta                                                                                                                         >  2:2674/1‑139 (MQ=255)
               ttCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCa                                                                                                                      <  2:669046/139‑1 (MQ=255)
                    tCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCAc                                                                                                                 >  2:946784/1‑139 (MQ=255)
                    tCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCAc                                                                                                                 >  2:483825/1‑139 (MQ=255)
                      cgcgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    >  1:710157/1‑134 (MQ=255)
                      cgcgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    <  2:710157/134‑1 (MQ=255)
                       gcgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCAcgcg                                                                                                              <  1:2674/139‑1 (MQ=255)
                        cgcgTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGt                                                                                                             >  1:453313/1‑139 (MQ=255)
                                   aaaTGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGct                                                                                                  <  2:537713/139‑1 (MQ=255)
                                           gcgcCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGa                                                                                                                                     <  2:587403/96‑1 (MQ=255)
                                           gcgcCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGa                                                                                                                                     >  1:587403/1‑96 (MQ=255)
                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAg                                                                                     <  1:729129/139‑1 (MQ=255)
                                                ccGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAg                                                                                     <  1:483825/139‑1 (MQ=255)
                                                   gTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTg                                                                                  <  2:453313/139‑1 (MQ=255)
                                                   gTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTg                                                                                  <  2:708515/139‑1 (MQ=255)
                                                                 gccgccCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAgc                                                                    <  1:946784/139‑1 (MQ=255)
                                                                             cGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATAcc                                                        >  1:254298/1‑139 (MQ=255)
                                                                               aCTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    >  2:635511/1‑77 (MQ=255)
                                                                               aCTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAgg                                                                                                                    <  1:635511/77‑1 (MQ=255)
                                                                                   cATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGgc                                                                                                     >  2:748200/1‑88 (MQ=255)
                                                                                   cATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGgc                                                                                                     <  1:748200/88‑1 (MQ=255)
                                                                                    aTCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAAc                                                 >  1:49511/1‑139 (MQ=255)
                                                                                     tCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACg                                                >  2:256761/1‑139 (MQ=255)
                                                                                          tCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGgaga                                                   >  1:423497/1‑131 (MQ=255)
                                                                                          tCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGgaga                                                   <  2:423497/131‑1 (MQ=255)
                                                                                            cATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGaa                                         <  1:715561/139‑1 (MQ=255)
                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGc                                                                                        >  2:663451/1‑82 (MQ=255)
                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGc                                                                                        <  1:663451/82‑1 (MQ=255)
                                                                                                      gTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGa                               >  1:623703/1‑139 (MQ=255)
                                                                                                       tCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGAc                              >  1:889561/1‑139 (MQ=255)
                                                                                                        cTGTTTGATGTTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACt                             >  2:135715/1‑139 (MQ=255)
                                                                                                           tttGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATa                                                          >  2:208143/1‑107 (MQ=255)
                                                                                                           tttGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATa                                                          <  1:208143/107‑1 (MQ=255)
                                                                                                             tGATGCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACtgtgtg                        <  2:397907/139‑1 (MQ=255)
                                                                                                                 gCTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGc                    >  1:529012/1‑139 (MQ=255)
                                                                                                                  cTGTAGCCGCCGGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGcc                   <  2:889561/139‑1 (MQ=255)
                                                                                                                          gccgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATaa           >  1:118352/1‑139 (MQ=255)
                                                                                                                           ccgGAGTCAAAGGCGATACCTTTACCTACCAGGAACGCGTACATTGGCGCATCTGTATCGCCAGTTGGGGTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGATCCGCGTCCGCCTGTGTGCAGCCCCAGATAAc          <  2:323943/139‑1 (MQ=255)
                                                                                                                            cgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGaa                                         >  2:614599/1‑107 (MQ=255)
                                                                                                                            cgGAGTCAAAAGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGaa                                         <  1:614599/107‑1 (MQ=255)
                                                                                                                                   aaaaGCGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCTCa  <  2:254298/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
AGCTTCACGCGCTTGTTCAGTCCGCGCGTAATGGCAAATGCCAGCGCCCCGGTAACCGTTGCCGCGCCGCCCATGTCCGACTTCATCGAGTCCATAAACGCAGTCTGTTTGATGCTGTAGCCGCCGGAGTCAAAAGTGATACCTTTACCTACCAGGCACGCGTACACTGGCGCTTCTTTATCGCCAGTTGGGTTGTAATCCAGCGCCAGCAATACCGGAGAACGTTCTGAACCGCGTCCGACTGTGTGCAGCCCCATATAACCTTGCTCA  >  NC_000913/2655631‑2655900

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