Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,521,249 A→C K315N (AAA→AAC yfeH → bile acid:sodium symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,521,2490AC97.7% 156.6 / ‑3.3 44K315N (AAA→AACyfeHbile acid:sodium symporter
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (27/16);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

GGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCA  >  NZ_CP009273/2521123‑2521378
                                                                                                                              |                                                                                                                                 
ggCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGaa                                                                                                                       <  2:837096/139‑1 (MQ=255)
 gCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGaa                                                                                                                       <  1:923030/138‑1 (MQ=255)
 gCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGa                                                                                                                        <  1:518654/137‑1 (MQ=255)
             ccGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGgcg                                                                                                           <  2:1050987/138‑1 (MQ=255)
                               ttCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGcc                                                                                         >  1:1032052/1‑138 (MQ=255)
                                    cacaTCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATa                                                                                     >  1:636575/1‑137 (MQ=255)
                                         cGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCtt                                                                               >  2:1195253/1‑138 (MQ=255)
                                           gTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGaa                                                                                                                       <  2:1122478/96‑1 (MQ=255)
                                           gTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGaa                                                                                                                       >  1:1122478/1‑96 (MQ=255)
                                                 ggTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       <  1:985134/138‑1 (MQ=255)
                                                 ggTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       <  2:619520/138‑1 (MQ=255)
                                                             ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAgg                                                                                                             >  2:700159/1‑88 (MQ=255)
                                                             ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAgg                                                                                                             >  2:948864/1‑88 (MQ=255)
                                                             ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAgg                                                                                                             <  1:948864/88‑1 (MQ=255)
                                                             ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAgg                                                                                                             <  1:700159/88‑1 (MQ=255)
                                                             ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       >  1:607855/1‑126 (MQ=255)
                                                             ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       <  2:607855/126‑1 (MQ=255)
                                                             ctgcCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCt                                                          >  1:673680/1‑139 (MQ=255)
                                                                    tgatgaTTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGAACAg                                                   >  2:732308/1‑139 (MQ=255)
                                                                             tCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGgcgc                                                                                                          >  2:523354/1‑75 (MQ=255)
                                                                             tCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGgcgc                                                                                                          <  1:523354/75‑1 (MQ=255)
                                                                                                ggTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATcag                        >  1:446791/1‑138 (MQ=255)
                                                                                                 gTCTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATcag                        >  1:322834/1‑137 (MQ=255)
                                                                                                   cTGTGCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCagcagc                    >  2:97415/1‑139 (MQ=255)
                                                                                                      tgCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCg                   <  1:955120/137‑1 (MQ=255)
                                                                                                      tgCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGt                  >  1:285184/1‑138 (MQ=255)
                                                                                                      tgCGGTGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGt                  >  2:1210907/1‑138 (MQ=255)
                                                                                                           tGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       >  2:716153/1‑80 (MQ=255)
                                                                                                           tGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGa                                                                       <  1:716153/80‑1 (MQ=255)
                                                                                                           tGCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTg            >  2:414044/1‑139 (MQ=255)
                                                                                                           tGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTg            >  1:746185/1‑139 (MQ=255)
                                                                                                            gCTGGCGCGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTg            >  1:405736/1‑138 (MQ=255)
                                                                                                                gcgcGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCt                                     <  2:353356/109‑1 (MQ=255)
                                                                                                                gcgcGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCt                                     >  1:353356/1‑109 (MQ=255)
                                                                                                                gcgcGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAtt        <  2:3090/138‑1 (MQ=255)
                                                                                                                gcgcGTCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCAt         >  2:425448/1‑137 (MQ=255)
                                                                                                                     tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAATCGGACGCTTCAGGGTCTGGACCAGTTGCGTCATCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCa  >  1:1052929/1‑139 (MQ=255)
                                                                                                                     tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGCTTTGCATTAGCTCa  >  2:1180702/1‑139 (MQ=255)
                                                                                                                     tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGAt               >  1:686671/1‑126 (MQ=255)
                                                                                                                     tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGAt               <  2:686671/126‑1 (MQ=255)
                                                                                                                     tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCa  >  2:480399/1‑139 (MQ=255)
                                                                                                                     tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCa  >  2:1005258/1‑139 (MQ=255)
                                                                                                                     tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCa  >  1:890928/1‑139 (MQ=255)
                                                                                                                     tCGATACAACCGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCa  >  1:1125617/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                 
GGCAAATGGCATCCCGATGGCAAACATTCTGTTCCCCACATCGGTGATCGGTATGATGGTGCTGCCCCTGATGATTTTCCATCAGATCCAATTGATGGTCTGTGCGGTGCTGGCGCGTCGATACAAACGCCAGACCGAACAGTTACAGGCGCAGCAGGAAAGCAGCGCCGATAAAGCTTAAAGCGGACGCTTCAGGGGCTGGACCAGTTGCGTCAGCCCCTCGGTTTTAATCAGCAGCGTGATTTGCATTAGCTCA  >  NZ_CP009273/2521123‑2521378

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