Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,499,480 2 bp→AC 100% 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,4800CA99.7% 114.6 / ‑7.3 43intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (22/20);  total (23/20)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0029471,499,4810AC100.0% 128.3 / NA 43intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (23/20);  total (23/20)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCT  >  NC_002947/1499339‑1499618
                                                                                                                                              ||                                                                                                                                             
cGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCtt                                                                                                                                         >  7:50957/1‑149 (MQ=255)
 gACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt                                                                                                                                        <  1:61944/149‑1 (MQ=255)
       ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                                  <  6:123915/149‑1 (MQ=255)
                gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                          >  8:23408‑M2/1‑141 (MQ=255)
                    aCGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAACAGGCGCCGCGGGGCATGGCACCCGCGTTGCCGGTGttcgcgggcacg                                                                                                                     <  2:204907‑M2/149‑13 (MQ=255)
                             cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctccc                                                                                                            >  1:204907‑M2/1‑128 (MQ=255)
                              ccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAGAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctccc                                                                                                            <  1:252945‑M2/149‑22 (MQ=255)
                                  gAGGCGGCCATTCAGAGCCGGAACGACGCCTTGCCACCACAGCAGCACTGAGGCTGCCAAACGGGGGGCGCGGCCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccacaga                                                                                                       >  1:233861‑M2/1‑123 (MQ=255)
                                   aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcaggcccgctcccacaga                                                                                                       >  8:54798‑M2/1‑122 (MQ=255)
                                   aTGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccacaga                                                                                                       >  1:231011‑M2/1‑122 (MQ=255)
                                      cGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccacagactga                                                                                                   >  8:204104‑M2/1‑119 (MQ=255)
                                                    cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                                   <  8:50957/149‑1 (MQ=21)
                                                    cTGAACGAAGCCTTGGCACCACAGCAGGAATGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGAAAAACGCGCCCCGGTTCATTGCACCGGATTTGCCGGTGTTAGAGGGCACGCCCGCGCCCACATCCTGAATGCCATTTATCgg                                                                                   >  8:157735/1‑149 (MQ=255)
                                                     tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCGTTGCCGGTGTTAGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCggg                                                                                  >  2:241822/1‑149 (MQ=25)
                                                     tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCCCGAAACAGGCGCCGCCGTGCCTGGCCCCCGCGTTGGCGGTGTTAGCGGGCACCCCCCCGCCCACAGCCGGATTGCCATTAACCggg                                                                                  >  7:276185/1‑149 (MQ=255)
                                                                          agcagGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACAAGCGCCGCCGTGCAAGGCACCCGCTTTGGCCGTGTTAGCGGGCACCCCCCCGCCCCCCGCCAGAATTCCAATTAGCGGGCCGATCGCGAAAGTggccggc                                                             >  6:283306/1‑143 (MQ=255)
                                                                                  tGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                     >  1:95204/1‑149 (MQ=21)
                                                                                  ggAGGCTGCCTATCGGTGGGTGCGGGGGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                      <  7:157735/147‑1 (MQ=21)
                                                                                    aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACTTGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                     <  1:237935/148‑1 (MQ=21)
                                                                                     ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                  <  2:252945/149‑1 (MQ=21)
                                                                                     ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                  <  5:283306/149‑1 (MQ=21)
                                                                                            aaaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCCCGGTGCATGGCACCGGCGTTGCCCGTGTTAGCGGGCCCGCCCCCTCCCACAGACTGATTTCCATTAAGCCGGCCGATCGCTACAGTGACGGCCCGGTTTTTGTTTCCTccc                                           >  6:246906/1‑147 (MQ=18)
                                                                                                   tGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGa                                     <  5:246906/149‑1 (MQ=21)
                                                                                                    gggAGCGGGCGTGCCCGGTAAACAGGCGCGGCTGTGCTTGGCACCGGCTTTGCCGTTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGaa                                    <  5:191783/149‑1 (MQ=25)
                                                                                                     ggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                                   <  3:284854/149‑1 (MQ=25)
                                                                                                      taGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                                  <  8:19517/148‑1 (MQ=25)
                                                                                                       aGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                                 >  3:271784/1‑149 (MQ=25)
                                                                                                        gcgggcggGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTTTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                <  7:121726/149‑1 (MQ=25)
                                                                                                        gCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                <  7:71391/149‑1 (MQ=25)
                                                                                                         cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                 >  6:242733/1‑99 (MQ=14)
                                                                                                         cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                 >  6:238461/1‑99 (MQ=14)
                                                                                                              gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                          <  5:135106/149‑1 (MQ=25)
                                                                                                               gccccgcGAAACCGGCGCCGCGGGGCATGGCACCGGCTTTGCCGGTGtgcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttg                           >  6:173830‑M2/3‑47 (MQ=255)
                                                                                                                gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                        >  8:160113/1‑149 (MQ=37)
                                                                                                                 cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctccccccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttg                           <  5:173830‑M2/147‑103 (MQ=255)
                                                                                                                  ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgt                          <  6:317540‑M2/147‑104 (MQ=255)
                                                                                                                  ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgggcacgcccgctcccaccgtttggcagcctcagtcctgctgtggtgccaaggcttcgttcaggctctgaatggccgcatcacgggcttcagcgttgt                          >  5:317540‑M2/1‑44 (MQ=255)
                                                                                                                  ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                        >  7:13180/1‑147 (MQ=37)
                                                                                                                  ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTa                      <  2:95204/149‑1 (MQ=37)
                                                                                                                    gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                     >  4:170815/1‑148 (MQ=37)
                                                                                                                      gAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAAc                  >  5:308662/1‑149 (MQ=37)
                                                                                                                                ccgcGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc        <  3:170815/149‑1 (MQ=37)
                                                                                                                                       gCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGTCCGGTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCt  >  1:134568/1‑148 (MQ=37)
                                                                                                                                              ||                                                                                                                                             
CGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCT  >  NC_002947/1499339‑1499618

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: