Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,160,541 G→A 9.1% A93A (GCG→GCA)  pabC → 4‑amino‑4‑deoxychorismate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,160,5410GA 9.1% 122.3 / 3.0 44A93A (GCG→GCA) pabC4‑amino‑4‑deoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base G (16/24);  new base A (2/2);  total (18/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTGGCCGGCCTCAAACACC  >  NC_002947/2160395‑2160686
                                                                                                                                                   |                                                                                                                                                 
gCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATgcgc                                                                                                                                                  <  2:159061/148‑1 (MQ=255)
   cGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCGCATAGCCACGCTGGCTAT‑CGCCacg                                                                                                                                              <  2:167796/149‑3 (MQ=255)
       tGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTc                                                                                                                                          >  7:248769/1‑149 (MQ=255)
         gCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTc                                                                                                                                          <  8:35129/147‑1 (MQ=255)
         gCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCg                                                                                                                                         <  3:84328/148‑1 (MQ=255)
         gCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGa                                                                                                                                        <  6:226341/149‑1 (MQ=255)
           aCTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATg                                                                                                                                      <  8:149591/149‑1 (MQ=255)
               ggCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCa                                                                                                                                    <  8:153112/147‑1 (MQ=255)
               ggCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAgc                                                                                                                                  >  6:141048/1‑149 (MQ=255)
                    ccAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGcgccgc                                                                                                                             >  6:57474/1‑149 (MQ=255)
                         gcCTGGCGATCAAGGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTTCGCCTGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCGAGCGTGGGTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGt                                                                                                                        <  6:114218/149‑1 (MQ=255)
                                    aaCGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAgg                                                                                                             <  3:105949/149‑1 (MQ=255)
                                                gcgcTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGTGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTa                                                                                                   <  1:31539/147‑1 (MQ=255)
                                                     ggTGCGCGACGAGCTGCTGCTCTTGGCCAGCCAGCTGGGTGAGGGTCTCGTCAATCTGATGCTTCCCGGTGGAG‑ATATCCAGCGGGGCTTTGCGCCAGCTGTTGTTGCAGCGCCTCGGTGTATCCTGCAAGGCAGGCCATTatgcagt                                                                                             <  8:186220/148‑7 (MQ=255)
                                                         cgcgACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGGCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATc                                                                                         <  5:26214/148‑1 (MQ=255)
                                                            gacgaGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGt                                                                                     >  8:46475/1‑149 (MQ=255)
                                                              cgaGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCg                                                                                   <  5:141048/149‑1 (MQ=255)
                                                                 gctgctGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGc                                                                                <  5:57474/149‑1 (MQ=255)
                                                                    gctgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGcc                                                                                                                                                 >  2:233414/1‑81 (MQ=255)
                                                                    gctgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGcc                                                                                                                                                 <  1:233414/81‑1 (MQ=255)
                                                                    gctgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATg                                                                             >  3:62710/1‑149 (MQ=255)
                                                                     ctgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATg                                                                             <  6:210881/148‑1 (MQ=255)
                                                                      tgcGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCt                                                                           >  3:221209/1‑149 (MQ=255)
                                                                        cgcTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGa                                                                         <  6:171266/149‑1 (MQ=255)
                                                                          cTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATGGCCAGCGTGGCTATGCACCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTAccct                                                                                               >  4:2448/1‑124 (MQ=255)
                                                                          cTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATGGCCAGCGTGGCTATGCACCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTAccct                                                                                               <  3:6279/125‑2 (MQ=255)
                                                                          cTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATGGCCAGCGTGGCTATGCACCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTAccct                                                                                               <  3:2448/125‑2 (MQ=255)
                                                                          cTACGCCAGCCAGCTGGGTCAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATGGCCAGCGTGGCTATGCACCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTAccct                                                                                               >  4:6279/1‑124 (MQ=255)
                                                                           tcacccagctggCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCa                                                                      >  4:214098/5‑149 (MQ=255)
                                                                                               gggTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCtt                                                  <  3:214098/149‑1 (MQ=255)
                                                                                                      gtcAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATAc                                           >  6:106551/1‑149 (MQ=255)
                                                                                                          aaCTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCgg                                       <  1:154501/149‑1 (MQ=255)
                                                                                                             tCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCtt                                    >  3:149982/1‑149 (MQ=255)
                                                                                                               aTCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTgg                                  <  2:156999/149‑1 (MQ=255)
                                                                                                                     aCCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTTGGGaa                              <  4:62710/147‑1 (MQ=255)
                                                                                                                        cGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGaacaac                          >  3:186544/1‑148 (MQ=255)
                                                                                                                                    ccAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTggccg              >  8:215314/1‑149 (MQ=255)
                                                                                                                                     cAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCtt                                                  >  1:145472/1‑112 (MQ=255)
                                                                                                                                     cAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCtt                                                  <  2:145472/112‑1 (MQ=255)
                                                                                                                                             cTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTGGCCGGCCTCAAac     >  5:170234/1‑149 (MQ=255)
                                                                                                                                              tATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTggc                >  2:103603/1‑137 (MQ=255)
                                                                                                                                              tATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTggc                <  1:103603/137‑1 (MQ=255)
                                                                                                                                              tATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTGGCCGGCCTCAAaca    >  5:47120/1‑149 (MQ=255)
                                                                                                                                                tGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTGGCCGGCCTCAAACAcc  <  6:47120/149‑1 (MQ=255)
                                                                                                                                                tGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTGGCCGGCCTCAAACAcc  <  4:221209/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCACGCCTGGCACTGGGCTGCCAGCGCCTGGCGATCAACGCCGATCTGGCGCTGGTGCGCGACGAGCTGCTGCGCTACGCCAGCCAGCTGGGTGAGGGTGTCGTCAAACTCATCCTTACCCGTGGCG‑ATAGCCAGCGTGGCTATGCGCCAGCTGTCGATGCAGCGCCGCGGCGTATCCTGCAAGGCAGCCCATTACCCCGTTATCCTGTCGAGCATGCTGAGCATGGCGTGCGTTTGTTCCCTTGCCATACCCGGCTTGGGGAACAACCACTGCTGGCCGGCCTCAAACACC  >  NC_002947/2160395‑2160686

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