Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 195,495 A→G E198G (GAG→GGG)  ispU → undecaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913195,4950AG82.5% 97.6 / 14.3 40E198G (GAG→GGG) ispUundecaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (2/5);  new base G (16/17);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

TATGCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGAGCATCGCATTAGTAACTTTTTGCTTTGGCAAATTGCCTATGCC  >  NC_000913/195448‑195538
                                               |                                           
tATGCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGc                                           <  1:762343/50‑1 (MQ=255)
  tGCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCAt                                         <  1:1007681/50‑1 (MQ=255)
  tGCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCAt                                         >  1:1498757/1‑50 (MQ=255)
  tGCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCAt                                         >  1:3206001/1‑50 (MQ=255)
  tGCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCAt                                         <  1:1773354/50‑1 (MQ=255)
  tGCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCAt                                         >  1:2529115/1‑50 (MQ=255)
   gCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATc                                        >  1:2796124/1‑50 (MQ=255)
    cATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATc                                        <  1:1186253/49‑1 (MQ=255)
    cATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCg                                       >  1:1084811/1‑50 (MQ=255)
        aaCTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCAtt                                   <  1:970204/50‑1 (MQ=255)
        aaCTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCAtt                                   <  1:124151/50‑1 (MQ=255)
        aaCTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCAt                                    >  1:3929534/1‑49 (MQ=255)
         aCTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTa                                  <  1:3080563/50‑1 (MQ=255)
         aCTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGAGCATCGCATTa                                  <  1:3691383/50‑1 (MQ=255)
          cTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAg                                 >  1:732112/1‑50 (MQ=255)
          cTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAg                                 >  1:3299404/1‑50 (MQ=255)
             gCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAGTa                               <  1:2793048/49‑1 (MQ=255)
               cccTGTAGATTTAGTAATTAGGACTGGGGGGGAGCATCGCATTAGTAACt                            <  1:878140/50‑1 (MQ=255)
                 cTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACttt                          <  1:2671002/50‑1 (MQ=255)
                 cTGTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACttt                          >  1:3781795/1‑50 (MQ=255)
                   gTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACttttt                        >  1:334138/1‑50 (MQ=255)
                   gTAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACttttt                        >  1:1610452/1‑50 (MQ=255)
                    tAGATTTTGTAATTAGGGCTGGGGGGGGGCATCGCATTAGTAACttttt                        <  1:33125/49‑1 (MQ=255)
                    tAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTg                       <  1:2219481/50‑1 (MQ=255)
                    tAGATTTAGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTg                       >  1:2000117/1‑50 (MQ=255)
                      gATTTAGTAATTAGGACTGGGGGGGAGCATCGCATTAGTAACTTTTTGc                      <  1:3193451/49‑1 (MQ=255)
                        tttAGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCttt                   >  1:2524297/1‑50 (MQ=255)
                           aGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTgg                 <  1:2608352/49‑1 (MQ=255)
                           aGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTgg                 <  1:3587496/49‑1 (MQ=255)
                           aGTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTgg                 <  1:536446/49‑1 (MQ=255)
                            gTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTGGc                >  1:3868480/1‑49 (MQ=255)
                            gTAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTGGc                <  1:2731404/49‑1 (MQ=255)
                             tAATTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTGGCaa              <  1:654797/50‑1 (MQ=255)
                              aaTTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTGGCaaa             >  1:3756964/1‑50 (MQ=255)
                               aTTAGGACTGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTGGCAAAt            <  1:2847476/50‑1 (MQ=255)
                                  aGGACTGGGGGGGAGCATCGCATTAGTAACTTTTTGCTTTGGCAAATTg          <  1:3071015/49‑1 (MQ=255)
                                  aGGACTGGGGGGGAGCATCGCATTAGTAACTTTTTGCTTTGGCAAATTGc         >  1:3672015/1‑50 (MQ=255)
                                      cTGGGGGGGAGCATCGCATTAGTAACTTTTTGCTTTGGCAAATTGCCTAt     >  1:2116672/1‑50 (MQ=255)
                                       tGGGGGGGGGCATCGCATTAGTAACTTTTTGCTTTGGCAAATTGCCTATg    >  1:2984191/1‑50 (MQ=255)
                                         ggggggggCATCGCATTAGTAACTTTTTGCTTTGGCAAATTGCCTATGcc  <  1:400624/50‑1 (MQ=255)
                                         ggggggAGCATCGCATTAGTAACTTTTTGCTTTGGCAAATTGCCTATGc   <  1:1308206/49‑1 (MQ=255)
                                               |                                           
TATGCATGAACTGGCCCCTGTAGATTTAGTAATTAGGACTGGGGGGGAGCATCGCATTAGTAACTTTTTGCTTTGGCAAATTGCCTATGCC  >  NC_000913/195448‑195538

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