Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 181,707 C→G G243A (GGC→GCC)  dapD ← 2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑succinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273181,7070CG95.5% 66.3 / ‑3.1 22G243A (GGC→GCC) dapD2,3,4,5‑tetrahydropyridine‑2,6‑dicarboxylate N‑succinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (11/10);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

GTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCACGGTCG  >  NZ_CP009273/181620‑181790
                                                                                       |                                                                                   
gTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGGCa                                                                                   <  2:475212/90‑1 (MQ=255)
       gcagTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGGCATCTTTTg                                                                            >  1:166992/1‑90 (MQ=255)
           ttCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGGCATCTTTTGACgg                                                                        <  2:197978/90‑1 (MQ=255)
                    gCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGGCATCTTTTGACGGCAGATTAcc                                                               >  2:140914/1‑90 (MQ=255)
                           ttGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCCCAGTAGAGGCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAaaca                                                        <  1:496993/90‑1 (MQ=255)
                                                    tCTTAACGATAACCGCACAGTAGAGGCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGAc                               <  2:410361/90‑1 (MQ=255)
                                                          cGATAACCGCACAGTAGAGGCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGt                         <  2:539945/90‑1 (MQ=255)
                                                                    acaGTAGAGGCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGAc                               <  2:232198/74‑1 (MQ=255)
                                                                    acaGTAGAGGCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGAc                               >  1:232198/1‑74 (MQ=255)
                                                                       gTAGAGGCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGt            >  1:251012/1‑90 (MQ=255)
                                                                          gagGCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGAc                               <  1:316315/68‑1 (MQ=255)
                                                                          gagGCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGAc                               >  2:316315/1‑68 (MQ=255)
                                                                            ggCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCAc       >  1:50209/1‑90 (MQ=255)
                                                                            ggCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCAc       <  1:554904/90‑1 (MQ=255)
                                                                            ggCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCAc       >  1:324087/1‑90 (MQ=255)
                                                                            ggCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCAc       >  2:14641/1‑90 (MQ=255)
                                                                            ggCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCAc       <  2:300331/90‑1 (MQ=255)
                                                                            ggCTGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCAc       <  2:509436/90‑1 (MQ=255)
                                                                               tGTATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCACGGt    >  2:506588/1‑90 (MQ=255)
                                                                                 tATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCACGGTCg  >  2:12053/1‑90 (MQ=255)
                                                                                 tATTTGGCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGGAGTGGATTTCGCCGGTTTCACGGTCg  >  1:299686/1‑90 (MQ=255)
                                                                                 tATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCGCAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCACGGTCg  >  1:468282/1‑90 (MQ=255)
                                                                                       |                                                                                   
GTACGCAGCAGTTCGTTAATGCCGACTTTGCCGCGAGTTTTCGCGTCAACTTTCTTAACGATAACCGCACAGTAGAGGCTGTATTTGCCATCTTTTGACGGCAGATTACCTGAAACAACCACAGACCCCGCCGGAACGCGACCGTAGTGGATTTCGCCGGTTTCACGGTCG  >  NZ_CP009273/181620‑181790

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