Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 4,009,285 T→C E82G (GAG→GGG)  ysgA ← dienelactone hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,009,2850TC92.3% 83.8 / 1.2 26coding (245/816 nt)ysgAdienelactone hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (7/17);  minor base . (1/1);  total (8/18)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATC  >  NZ_CP009273/4009160‑4009412
                                                                                                                             |                                                                                                                               
gaTCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACg                                                                                                                    <  2:228403/139‑1 (MQ=255)
      ccAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGAc                                                                                                               >  1:293614/1‑138 (MQ=255)
             cTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTcg                                                                                                       <  1:14939/139‑1 (MQ=255)
             cTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTcg                                                                                                       <  1:348009/139‑1 (MQ=255)
                        gCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCa                                                                                           >  2:120418/1‑139 (MQ=255)
                           cTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATg                                                                                         <  2:233576/139‑1 (MQ=255)
                                        aCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTcgc                                                                                                      <  1:104629/113‑1 (MQ=255)
                                        aCCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTcgc                                                                                                      >  2:104629/1‑113 (MQ=255)
                                          cGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTc                                                                          <  1:210166/139‑1 (MQ=255)
                                               aGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCa                                                                     <  2:386883/139‑1 (MQ=255)
                                                          gATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCa                                                         <  1:120418/139‑1 (MQ=255)
                                                           atatCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCa                                                         >  1:245081/1‑139 (MQ=255)
                                                           atatCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCa                                                         >  2:210530/1‑139 (MQ=255)
                                                              tCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACt                                                       <  2:445399/138‑1 (MQ=255)
                                                               cGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTgg                                                     >  2:401744/1‑139 (MQ=255)
                                                                                        gcgGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGt                            <  1:353915/139‑1 (MQ=255)
                                                                                        gcgGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGt                            <  2:293614/139‑1 (MQ=255)
                                                                                              gTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGctttgctttgtt                            >  2:144883/1‑131 (MQ=255)
                                                                                              gTAAAGTTCAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGt                            <  1:144883/133‑1 (MQ=255)
                                                                                                      cAGGTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAg               >  2:81357/1‑138 (MQ=255)
                                                                                                        ggTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAg            <  1:468584/139‑1 (MQ=255)
                                                                                                        ggTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAg            <  2:237539/139‑1 (MQ=255)
                                                                                                         gTGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAgg           <  1:401744/139‑1 (MQ=255)
                                                                                                          tGCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGc          <  1:210530/139‑1 (MQ=255)
                                                                                                           gCGATAGCCAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCa         <  2:225712/139‑1 (MQ=255)
                                                                                                                  ccAGATACCCCCCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATc  <  1:208225/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                               
GATCGGCCAGCACCTGCGAGTCAGGCACTTTTGCTACCAGACCGCTAAGCAACGTGGGGATATCGGCAAAATCATTCGGATCGCCTTCGCGGAAGTAAAGTTCAGGTGCGATAGCCAGATACCCCTCCAGCGCCAGACGGCGACAAATGTCGCGGATATGTTCATGCACGCCAAAAATTTCCTGCACTACAATGACCACTGGCAGTGGGCCATCGCTTTGCTTTGGTCTGGCATGGTAAGCAGGCATGTTATC  >  NZ_CP009273/4009160‑4009412

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