Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,480 2 bp→AC intergenic (+57/‑106) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,499,4800CA97.4% 105.6 / ‑6.8 39intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (14/24);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4810AC97.4% 109.1 / NA 39intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (13/25);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACG‑‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCA  >  NC_002947/1499338‑1499613
                                                                                                                                              ||                                                                                                                                            
ccGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCt                                                                                                                                         >  1:47151/1‑149 (MQ=255)
      ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                   <  2:15566/149‑1 (MQ=255)
       gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCgg                                                                                                                                  >  8:9466/1‑149 (MQ=255)
        ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                                 <  4:47155/149‑1 (MQ=255)
        ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                                 <  6:212805/149‑1 (MQ=255)
          gCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt                                                                                                                               <  2:45015/149‑1 (MQ=255)
           cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtgcc                                                                                                                              <  8:232137/149‑3 (MQ=38)
                 gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                        <  8:203465/148‑6 (MQ=255)
                               ccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTCGCGGGCACCC‑‑‑CCGCTcccc                                                                                                     >  6:427297/1‑148 (MQ=25)
                                    aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGAc                                                                                               <  4:113173/149‑1 (MQ=21)
                                         gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGAtt                                                                                          <  6:307836/149‑1 (MQ=21)
                                             ttCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCa                                                                                      <  2:47151/149‑1 (MQ=21)
                                                    ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAACGCGCCGCGGTGCAAGGCACCGGCCTTGCCCGTGTTCGCGGGCACGC‑‑CCCCGCCCACCGACGGATTGCCATTAAGCg                                                                               >  3:322412/1‑149 (MQ=21)
                                                    ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACACGCGCCCCCGTGCATGGCACCGGCGTTGCCGGTGTTAGCGGGCCCGC‑‑CCGCGCCCACAAACTGATTGCCAATAACCg                                                                               >  8:382161/1‑149 (MQ=255)
                                                        aaCGAAGCCTTGGCACCACAGCAGGACTGCGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTGCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGcc                                                                           <  7:382161/149‑1 (MQ=21)
                                                           gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAAGGCGCCCCCGTGCAAGGCACCCGCATTGCCGGTGTTAGCGGGCACCC‑‑‑CCCCACCCACACACAGAATTCCAATAATCggggcaa                                                                         >  6:376745/1‑144 (MQ=16)
                                                            aaGCCTTGGCAACACAGCAGGACTGAGGCTGCCAAACGGTGGGAGAGGGCCTAACCCCGAAACACGCGCCCCCGTGAAATGCAACAGCATTGCCGGTGTTCGCGGGCACGC‑‑CCGCTACCAAAAACTGGTTGCCATTAATCGAGCCAATc                                                                       >  6:112311/1‑149 (MQ=255)
                                                             aGCCTGGGCACCACAGCAGGCCTGAGGCTGCCAACCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                      <  3:44334/149‑1 (MQ=18)
                                                              gCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGc                                                                     <  7:9466/149‑1 (MQ=21)
                                                                          cagcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                   <  5:149742/81‑1 (MQ=255)
                                                                          cagcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                   >  6:149742/1‑81 (MQ=255)
                                                                             cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                            <  2:66569/127‑1 (MQ=18)
                                                                             cagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                            >  1:66569/1‑127 (MQ=18)
                                                                                     aGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttg                                               <  5:376745/148‑1 (MQ=21)
                                                                                        cTGCCAAACGGTGGTTGCGGGCGTTCCCGCTAAACAGGCGCCGCTGTGCATGGCACGGGCTTTGCGGGTTTTCGCTGGAACGC‑‑CCTCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGc                                           <  1:371091/149‑1 (MQ=255)
                                                                                          gCCAAACGGTGGGTGGGGGCGTGCCCGCTAAACAGGCGCGGCGGTGCTTGGCACCGGCTTTGCCGGTGTTCGGGGGCACGC‑‑CCGCTCCCACAGCACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                <  1:97892/143‑1 (MQ=21)
                                                                                                 ggggggCGCGGGCGTGCCCGCGAAACAGGGGCCGCGGTGCATGGGACCGGGTTTGCCGGTGTTGGCGGGCACGCCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttt                                                  <  1:323251/132‑1 (MQ=21)
                                                                                                          gggCGTGCCCGCTAAACTGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                         <  3:349226/149‑1 (MQ=37)
                                                                                                          gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                         <  4:322412/149‑1 (MQ=37)
                                                                                                                 cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                             <  2:2673/138‑1 (MQ=25)
                                                                                                                 cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                             >  1:2673/1‑138 (MQ=25)
                                                                                                                 cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                           >  7:359184/1‑140 (MQ=25)
                                                                                                                 cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                           <  8:359184/140‑1 (MQ=25)
                                                                                                                  ccGCTAAACTGGCGCCGCGGTGCATGGCCCCGGCTTTGCGGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                 <  2:144332/149‑1 (MQ=37)
                                                                                                                         aCAGGCGCCGCGGTGCATGGCACCGGCTGTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt          <  5:427297/149‑1 (MQ=37)
                                                                                                                          cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa         >  5:427143/1‑149 (MQ=37)
                                                                                                                          cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa         >  1:347758/1‑149 (MQ=37)
                                                                                                                             gcgcCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACG‑‑‑CCAGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc       >  6:392407/1‑148 (MQ=38)
                                                                                                                                 cgcgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCa  <  6:427143/149‑1 (MQ=37)
                                                                                                                                              ||                                                                                                                                            
CCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACG‑‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCA  >  NC_002947/1499338‑1499613

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