Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,438,552 Δ1 bp coding (488/810 nt) yfcA ← sulfite exporter TauE/SafE family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,438,5510G.81.6% 105.6 / 27.1 38coding (489/810 nt)yfcAsulfite exporter TauE/SafE family protein
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base . (16/15);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

TAAACCGCCGATGTTTGACGTTGCGTTGAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACA  >  NZ_CP009273/2438446‑2438685
                                                                                                         |                                                                                                                                      
tAAACCGCCGATGTTTGACGTTGCGTTGAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGa                                                                                                      >  1:217769/1‑139 (MQ=255)
       ccgATGTTTGACGTTGCGTTGAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGc                                                                                                >  2:617454/1‑138 (MQ=255)
       ccgATGTTTGACGTTGCGTTGAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGc                                                                                                >  1:209215/1‑138 (MQ=255)
                       cgttgAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAAcgc                                                                                <  2:82693/139‑1 (MQ=255)
                        gttgAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGa                                                                              <  1:192948/139‑1 (MQ=255)
                              tAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTa                                                                        <  1:364322/139‑1 (MQ=255)
                              tAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTa                                                                        <  2:542366/139‑1 (MQ=255)
                                    ggCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCg                                                                   >  2:363772/1‑138 (MQ=255)
                                        tGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCgcag                                                                                                                                     <  2:801922/69‑4 (MQ=255)
                                        tGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCgcag                                                                                                                                     >  1:801922/1‑66 (MQ=255)
                                          aGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATc                                                             >  2:158406/1‑139 (MQ=255)
                                               tGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATa                                                                                                              >  1:206750/1‑84 (MQ=255)
                                               tGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATa                                                                                                              <  2:206750/84‑1 (MQ=255)
                                               tGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGgcgc                                                       >  2:780464/1‑139 (MQ=255)
                                               tGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGgcgc                                                       >  1:375576/1‑139 (MQ=255)
                                                  cTTTGGCGAGGTTAAAGCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTgg                                                    <  2:499305/139‑1 (MQ=255)
                                                          aGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGtcttc                                            >  2:689143/1‑139 (MQ=255)
                                                                 aTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCAcc                                      <  2:82936/139‑1 (MQ=255)
                                                                 aTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCAcc                                      >  2:743982/1‑138 (MQ=255)
                                                                 aTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCt                                     <  1:795089/139‑1 (MQ=255)
                                                                      cacaGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCATAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGa                                                                              >  1:629221/1‑93 (MQ=255)
                                                                      cacaGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGa                                                                              <  2:629221/93‑1 (MQ=255)
                                                                      cacaGCGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCtt                                >  1:84620/1‑139 (MQ=255)
                                                                           cGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCgg                                                          >  2:273422/1‑108 (MQ=255)
                                                                           cGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCgg                                                          <  1:273422/108‑1 (MQ=255)
                                                                           cGTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCa                           >  1:350759/1‑139 (MQ=255)
                                                                            gTAACGAAGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCAt                          <  1:689143/139‑1 (MQ=255)
                                                                                  aaGGCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCaa                    <  2:217769/139‑1 (MQ=255)
                                                                                   aGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGc                      >  1:672258/1‑137 (MQ=255)
                                                                                     gCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATc                         <  1:581215/132‑1 (MQ=255)
                                                                                     gCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATc                         >  2:581215/1‑132 (MQ=255)
                                                                                     gCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCAt                          <  2:264664/130‑1 (MQ=255)
                                                                                     gCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCAt                          >  1:264664/1‑130 (MQ=255)
                                                                                     gCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGa                  >  1:175003/1‑138 (MQ=255)
                                                                                     gCCAGCGCGTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGa                  >  1:225953/1‑138 (MQ=255)
                                                                                             gTAAAACGACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGc         <  2:225953/139‑1 (MQ=255)
                                                                                                 aaCGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGa      <  1:158406/139‑1 (MQ=255)
                                                                                                    gACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACa  <  1:780464/139‑1 (MQ=255)
                                                                                                    gACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACa  <  2:209215/139‑1 (MQ=255)
                                                                                                    gACCCGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACa  <  2:350759/139‑1 (MQ=255)
                                                                                                         |                                                                                                                                      
TAAACCGCCGATGTTTGACGTTGCGTTGAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCATAAAAACCGACGCAGCCACCGGCTATTAACGCGAAGGGTAAACCGTACATCCGGCGCTGGCGGTCTTCTTCACCTAGCTTTGGCATCAGCAAGAAATAGAGGCCGATACA  >  NZ_CP009273/2438446‑2438685

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