Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 6,672 G→A L430F (CTC→TTC)  yaaJ ← sodium:alanine symporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092736,6720GA94.3% 184.0 / 0.9 53L430F (CTC→TTC) yaaJsodium:alanine symporter family protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (27/23);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGC  >  NZ_CP009273/6552‑6807
                                                                                                                        |                                                                                                                                       
aCGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCaa                                                                                                                       <  2:332929/139‑1 (MQ=255)
 cgcgTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCaaa                                                                                                                      >  1:548191/1‑139 (MQ=255)
       cGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCaaa                                                                                                                      >  2:920160/1‑133 (MQ=255)
       cGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCaaa                                                                                                                      <  1:920160/133‑1 (MQ=255)
       cGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTTGt                                                                                                                 >  1:962507/1‑138 (MQ=255)
           gaAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAg                                                                                                            <  2:757250/139‑1 (MQ=255)
                cTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATg                                                                                                       <  2:250188/139‑1 (MQ=255)
                    cGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAgg                                                                                                   <  1:72198/139‑1 (MQ=255)
                          atatCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTa                                                                                             <  2:1285867/139‑1 (MQ=255)
                          atatCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTa                                                                                             <  1:1072550/139‑1 (MQ=255)
                                 gATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGcc                                                                                                 >  2:43307/1‑128 (MQ=255)
                                 gATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGcc                                                                                                 <  1:43307/128‑1 (MQ=255)
                                     aCGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGc                                                                                   >  1:886544/1‑138 (MQ=255)
                                                cAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACAt                                                                       >  1:1286802/1‑139 (MQ=255)
                                                      ccGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCgg                                                                 <  1:899968/139‑1 (MQ=255)
                                                       cGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCgg                                                                 >  1:968825/1‑138 (MQ=255)
                                                       cGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGa                                                                <  1:84732/139‑1 (MQ=255)
                                                        gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACAt                                                                       >  2:874871/1‑131 (MQ=255)
                                                        gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACAt                                                                       <  1:874871/131‑1 (MQ=255)
                                                        gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCgg                                                                 >  2:34608/1‑137 (MQ=255)
                                                        gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCgg                                                                 >  1:1067159/1‑137 (MQ=255)
                                                        gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCgg                                                                 >  1:733151/1‑137 (MQ=255)
                                                        gggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCgg                                                                 >  1:402255/1‑137 (MQ=255)
                                                         ggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGt                                                                                                                          >  2:477705/1‑79 (MQ=255)
                                                         ggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGt                                                                                                                          <  1:477705/79‑1 (MQ=255)
                                                         ggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGc                                                                                                  >  2:576182/1‑103 (MQ=255)
                                                         ggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGc                                                                                                  <  1:576182/103‑1 (MQ=255)
                                                         ggCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAg                                                              >  2:451585/1‑139 (MQ=255)
                                                                    aGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCa                                                                                                     <  1:957721/89‑1 (MQ=255)
                                                                    aGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCa                                                                                                     >  2:957721/1‑89 (MQ=255)
                                                                    aGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAgg                                                   >  1:170188/1‑139 (MQ=255)
                                                                              ggCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGAt                                         >  2:236456/1‑139 (MQ=255)
                                                                               tggTAGAGAATCACTGGCAATGGTATGAACCACAGGCGAGAAGAGTAAGATGTCAGTCAAATTGGGAATGGCCATGCAGGCCATTATGATGTCTGCCAGTGGCCACATCAGCGGATGGCTTAGCAAGGTGCCGCCGAt                                         <  2:310502/136‑1 (MQ=255)
                                                                                         tCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAg                                                                                                            <  2:810408/61‑1 (MQ=255)
                                                                                         tCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAg                                                                                                            >  1:810408/1‑61 (MQ=255)
                                                                                         tCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCg                               >  2:170943/1‑138 (MQ=255)
                                                                                         tCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCg                               >  2:144768/1‑138 (MQ=255)
                                                                                         tCACTGGCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAACAAGGTGCCGCCGATGACCGTTGCg                               >  2:1055918/1‑138 (MQ=255)
                                                                                               gCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTg                         <  2:606874/138‑1 (MQ=255)
                                                                                               gCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTg                         >  1:1281730/1‑138 (MQ=255)
                                                                                               gCAATGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTg                         >  1:748072/1‑138 (MQ=255)
                                                                                                  aTGGTATGAACCACAGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGa                     <  1:236499/139‑1 (MQ=255)
                                                                                                   tGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTg                         >  2:532996/1‑134 (MQ=255)
                                                                                                               caGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCcaca                                                                        >  2:533996/1‑75 (MQ=255)
                                                                                                               caGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCcaca                                                                        <  1:533996/75‑1 (MQ=255)
                                                                                                               caGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCa               <  1:667777/132‑1 (MQ=255)
                                                                                                               caGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCa               >  2:667777/1‑132 (MQ=255)
                                                                                                                aGGCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAacac          >  1:754301/1‑136 (MQ=255)
                                                                                                                  gCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGAt     <  2:754301/139‑1 (MQ=255)
                                                                                                                  gCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGAt     <  2:402255/139‑1 (MQ=255)
                                                                                                                  gCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGAt     <  1:34608/139‑1 (MQ=255)
                                                                                                                  gCGAGAACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGAt     <  1:451585/139‑1 (MQ=255)
                                                                                                                    gagaACAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCg   >  1:1002294/1‑139 (MQ=255)
                                                                                                                     agaaCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGc  >  2:705830/1‑139 (MQ=255)
                                                                                                                        |                                                                                                                                       
ACGCGTCCGGAGAAAGCTGGCGACCGATATCCGGATAACGCAATGGATCAAACACCGGGCGCACGCCGAGTTTACGCTGGCGTAGATAATCACTGGCAATGGTATGAACCACAGGCGAGAGCAGTAAAATGGCGGTCAAATTGGTAATAGCCATGCAGGCCATTATGATATCTGCCAGTTGCCACATCAGCGGAAGGCTTAGCAAGGTGCCGCCGATGACCGTTGCGAAGGTGCAGATCCGCAAACACCAGATCGC  >  NZ_CP009273/6552‑6807

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