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,4800CA92.9% 66.8 / ‑6.0 28intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (10/16);  total (10/18)
Fisher's exact test for biased strand distribution p-value = 5.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4810AC96.4% 64.1 / ‑6.4 28intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (10/17);  total (10/18)
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

CGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑‑TCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCT  >  NC_002947/1499347‑1499618
                                                                                                                                     ||                                                                                                                                                
cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                                     <  5:45049/148‑1 (MQ=255)
cGCAGGGCGACAACGCTGAAGCCCGTGATGCGGCCCTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                                     <  5:146706/148‑1 (MQ=38)
   aGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGtt                                                                                                                                 <  5:152829/148‑1 (MQ=255)
        gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcg                                                                                                                             <  6:98158/147‑5 (MQ=255)
              gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGACCGACGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgc                                                                                                                  <  3:2072/148‑1 (MQ=25)
                ggAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgccc                                                                                                                <  7:401/147‑2 (MQ=37)
                    gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcgcccgctc                                                                                                            >  7:151903/1‑143 (MQ=25)
                      ccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCGCGGGGCATGTCACCCGCGTTGCCGGTGTT‑‑TGCGGGCCcgcccgcgccc                                                                                                          >  1:87849/1‑141 (MQ=25)
                      ccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCCcgcccgcgccc                                                                                                          >  3:140760/1‑141 (MQ=25)
                                            cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCCGCTTTGCCGGTGTTCGCCGGCACGCCCGCTCCCAAAGACTGATTTCCATTAAGcc                                                                                     >  8:153621/1‑147 (MQ=21)
                                                aCGAAGCCTGGGCACCACAGCAGGACTGAGGCGGCCAAACGGGGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCGGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                                  <  1:155464/147‑1 (MQ=21)
                                                     gCCTTGCCACCCCCGCCGCCCTGAGGCTGCAAAAGGGTGGGAGCGGGCGGCCCCGGTAAACTGGCGCCGCGGGGATGGGCCCCGGCTTTGCGGGTGTG‑‑GGGGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                            <  8:65936/148‑1 (MQ=255)
                                                                          tGAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                      <  5:89337/148‑1 (MQ=21)
                                                                             ggCTGCCAAACGGTGGGTTGGGGCGTGCACGCGAAAACGGCGCGGCGGTTCATGGAACGGGGTTTGCCGGGGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAATCGTGCCTATCGCTACAGTGACGTCCCGtttttt                                                     <  4:163586/146‑1 (MQ=255)
                                                                             ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                     <  6:84121/146‑1 (MQ=21)
                                                                              gCTGCCTAAGGGTGGGCGCGGGCGTGCCCGCGAAACCGGCGCCGCGGGGAATGGCACCGGCTTGGCCGGTGTGC‑‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                  <  7:6161/148‑1 (MQ=18)
                                                                              gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                  <  8:151903/148‑1 (MQ=21)
                                                                                    aaaCGGTGGGAGCGGGCGTGCCCCCCAAACACGCGCCCCGGTGAATGGCACCCGCTTTGGCGGTGTTAGCGGGCCCCCCCGCGCCCCCACACTGATTGCCAATAAGCCGGCCGATCGCTAAAGTGACCGCCCGGTTTTTTTTTGcgc                                              >  2:113712/1‑145 (MQ=255)
                                                                                     aaCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATTGCACCCGCTTTGCCGGTGTTCGCGGGCCCCCCCGCTCCCACAGACTGATTGCCATAAACCGGGCCTATCGCTACACTGACGGCCCGGTTTTTGTTGCCTcc                                            >  6:94508/1‑146 (MQ=18)
                                                                                             gAGCGGGCGTGCCCGCGAAAAACGCGCCCCGGTGCATGGCACCCGCGTTTCCGGTGTTAGCGGGCACGCCCGCGACCACAGACAGATTTCCATTAAGCGGGCCTATAGCGACAGTGACGGCCCGGTTTTTGTTTCCTCGCTAACGaa                                    >  8:136342/1‑147 (MQ=16)
                                                                                                ggggggggCCGGCTAAACTGGCGCGGGGTGGCTTGCCACGGGTTTGGCGGTTTTGCGCGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                <  2:87849/141‑1 (MQ=255)
                                                                                                cGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                                <  7:153621/148‑1 (MQ=25)
                                                                                                         ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                                         >  7:64919/1‑99 (MQ=17)
                                                                                                          cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTAc                                                                        <  8:64919/99‑1 (MQ=17)
                                                                                                             gAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGaa                   >  1:67322/1‑148 (MQ=37)
                                                                                                                 cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt               >  7:152587/1‑148 (MQ=37)
                                                                                                                          cgGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAg      <  1:113712/148‑1 (MQ=37)
                                                                                                                              gCATGGCACGGGCTTGGCCGGTGTTCGCGGGAACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATTACTCCCTTAGAACGCTATCAAGCAGATCt  <  8:137010/148‑1 (MQ=255)
                                                                                                                                     ||                                                                                                                                                
CGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTT‑‑TCGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGATCT  >  NC_002947/1499347‑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: