Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,708,538 A→G intergenic (+10/‑96) dtpA → / → gstA dipeptide/tripeptide permease DtpA/glutathione transferase GstA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,708,5380AG94.7% 47.1 / ‑3.2 19intergenic (+10/‑96)dtpA/gstAdipeptide/tripeptide permease DtpA/glutathione transferase GstA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (7/11);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CGTATTCTTGCAGATTGGTGTCGCTACTGCCGTTATTGCAGTACTGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGT  >  NZ_CP009273/1708402‑1708657
                                                                                                                                        |                                                                                                                       
cgTATTCTTGCAGATTGGTGTCGCTACTGCCGTTATTGCAGTACTGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGgaa                                                                                                                       >  2:766550/1‑139 (MQ=255)
cgTATTCTTGCAGATTGGTGTCGCTACTGCCGTTATTGCAGTACTGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGgaa                                                                                                                       >  2:371854/1‑139 (MQ=255)
               tGGTGTCGCTACTGCCGTTATTGCAGTACTGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCg                                                                                                        >  2:392825/1‑139 (MQ=255)
                       cTACTGCCGTTATTGCAGTACTGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACttt                                                                                                <  1:638466/139‑1 (MQ=255)
                              cGTTATTGCAGTACTGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAg                                                                                         <  1:766550/139‑1 (MQ=255)
                                     gCAGTACTGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCttt                                                                                  >  2:210577/1‑139 (MQ=255)
                                              aTGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGtt                                                                         >  2:654869/1‑139 (MQ=255)
                                                      gACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACg                                                                 >  1:186340/1‑139 (MQ=255)
                                                        ccGCGCCGAAACTGCCCCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCa                                                               <  1:210577/139‑1 (MQ=255)
                                                         cgcgCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTTAGCAACGccc                                                              >  1:62471/1‑137 (MQ=255)
                                                             ccGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTAcc                                                          <  1:751516/139‑1 (MQ=255)
                                                                     gCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAACCAGCCGTGGCGTAAATTTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTa                                                  <  2:88771/139‑1 (MQ=255)
                                                                          gCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACca                                             <  1:721906/139‑1 (MQ=255)
                                                                           cATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACca                                             >  2:310707/1‑138 (MQ=255)
                                                                            aTGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACcaca                                           <  2:1272/139‑1 (MQ=255)
                                                                                       ggACGCTGCAGAAAAAGGGGCGAAAGCAGCCGGAGCGTAAAATTCGGGGGACCTCTTTTCAAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAg                                <  2:62471/138‑1 (MQ=255)
                                                                                          cGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACg                                                                 >  2:386876/1‑103 (MQ=255)
                                                                                          cGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACg                                                                 <  1:386876/103‑1 (MQ=255)
                                                                                               cAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGAt                        >  2:439971/1‑139 (MQ=255)
                                                                                                                  agcCGTAGCGTAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCg     <  1:145165/139‑1 (MQ=255)
                                                                                                                     cGTAGCGAAAAATTCAGGGGAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGt  <  1:310707/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                       
CGTATTCTTGCAGATTGGTGTCGCTACTGCCGTTATTGCAGTACTGATGCTGCTGACCGCGCCGAAACTGCACCGCATGACGCAGGATGACGCTGCAGACAAAGCGGCGAAAGCAGCCGTAGCGTAAAATTCAGGGAAACTCTTTTACAAGCCGCTAACTTTTCGTTAGCGGCTTTTTTTTTGTTCAGCAACGCACTACCATACTTTAAACCACAGACAAAAGGAGTTACCGATGAAATTGTTCTACAAACCGGGT  >  NZ_CP009273/1708402‑1708657

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