Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,667,677 G→A 41.9% A124V (GCG→GTG)  yphB ← mutarotase superfamily protein, YphB family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,667,6770GA41.9% 8.5 / 17.0 23A124V (GCG→GTG) yphBmutarotase superfamily protein, YphB family
Reads supporting (aligned to +/- strand):  ref base G (6/7);  new base A (6/3);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

ACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCA  >  CP009273/2667531‑2667819
                                                                                                                                                  |                                                                                                                                              
aCTGCTCCCGCTCCAGCCAGTAACCGAAGGGCGGTCCCAAAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCCATGCCCCTTGATTGGTGGCACACAGCGTCACCGTCACCGTATccgcc                                                                                                                                              >  2:175082/1‑149 (MQ=255)
      cccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGCGAGCGTCACCGTCAGCGTAGCCGCCGTTaaa                                                                                                                                        >  1:307205/1‑149 (MQ=255)
       ccGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCACGGTTAAAt                                                                                                                                       >  2:182913/1‑149 (MQ=255)
        cGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGCATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGGCAGAGGGCGTCACCGTCAGCGTATCCGCCGTTAAATg                                                                                                                                      >  2:174588/1‑149 (MQ=255)
           tCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGaaa                                                                                                                                   <  2:94888/149‑1 (MQ=255)
                        cGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGGACCAAGTGGCAGCGTCTCTGCCCCTTGATGGGTGACAGAGAGCGTCACCGTCAGCGTATCCACCGTTGCATGAAACGCCTGACTGAcc                                                                                                                      >  2:146513/1‑149 (MQ=255)
                        cGCTCGCCTGTGCCCGAGTCCGCTTTGGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTAcc                                                                                                                      <  2:7147/149‑1 (MQ=255)
                                     cTGAACCCGCGTTTGCGACGGAGACCGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACATAGAGCGTCACCGTCAGCGTATCCGCTGTTAAATTAAACGCCTGACTTACCCGATAGTGATAGa                                                                                                         <  2:232651/149‑1 (MQ=255)
                                       gAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTATGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCACCGTTCAATAAAACCCCTGACTGACCCGCTaaatatataca                                                                                                       >  2:117409/1‑139 (MQ=255)
                                           ccGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGc                                                                                                   <  1:281274/149‑1 (MQ=255)
                                                             cGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATTGTGATAGACACCGCTGCGATGTTCATACACCa                                                                                 <  1:256567/149‑1 (MQ=255)
                                                                      aGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCACCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAAc                                                                        >  1:135384/1‑149 (MQ=255)
                                                                                ccGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAACGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCt                                                              <  2:226414/149‑1 (MQ=255)
                                                                                         aaTGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATATTCATCCCCAAAACACAAACTatcaccccccagttagac                                                     >  2:82078/1‑135 (MQ=255)
                                                                                           tGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCACCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCCTCGCTGTGCGAGacac                                                   <  2:53413/149‑1 (MQ=255)
                                                                                                 cGTCCCTGCCCCTTGATTGGTGACAGAGAGCGGCCCCGTCAGCGTATCCACCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGcc                                             <  1:182913/149‑1 (MQ=255)
                                                                                                    ctctGCCCCTTGATTGGTGACAGAGAGCGGCACCGTCAGCGTATCCACCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACt                                          <  1:167264/149‑1 (MQ=255)
                                                                                                     tctGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCACCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTc                                         >  2:231498/1‑149 (MQ=255)
                                                                                                       tGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCACCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGc                                       >  2:202359/1‑149 (MQ=255)
                                                                                                         ccccTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGGGCCACATTGCCACTCGccc                                     >  2:141698/1‑149 (MQ=255)
                                                                                                          cccTTGATTGGTGACAGAGAGCGACACCGTGACCGTAGCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCcca                                    >  1:345136/1‑149 (MQ=255)
                                                                                                              tGATTGGTGACGGAGAGCGTCACCGCCAGCGTATCCCCCGTTATATGAAACTCCTGCCTTACCCGATAGTGATAGACACCGCTGCCATGTTCATACACCAGACACAATCTATCATCGCTGTGCGAGACACATTGCCCCTCGCccagcca                                <  1:146513/149‑1 (MQ=255)
                                                                                                                                   aCCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCAGCGCCGGGTAGATAGTg           >  1:268530/1‑149 (MQ=255)
                                                                                                                                            gTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCa  <  2:247444/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
ACTGCTCCCGCTCCAGCCAGTAACCGCTCGCCTGTGCCTGAATCCGCGTTTGCGGCGACAACGGAAAATAAGGATGCCAGCCGGTACCAAATGGCAGCGTCTCTGCCCCTTGATTGGTGACAGAGAGCGTCACCGTCAGCGTATCCGCCGTTAAATGAAACGCCTGACTTACCCGATAGTGATAGACACCGCTGCGATGTTCATACACCAGACACAAACTATCATCGCTGTGCGAGACACATTGCCACTCGCCCAGCCAGCCATCGCCGTGTAGATAGTGCGCATCCCA  >  CP009273/2667531‑2667819

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