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,4800CA96.0% 134.2 / ‑3.9 50intergenic (+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 (24/24);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4810AC96.0% 140.0 / ‑6.1 50intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (24/24);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCCA‑C‑GGCTTTGCCGGTGTT‑TCGGG‑CACG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAT  >  NC_002947/1499340‑1499616
                                                                                                                                             ||                                                                                                                                             
gACGGGCCGCAGTGCGACTACGCTGAACCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCttt                                                                                                                                      <  8:267858/149‑1 (MQ=38)
  cGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGc                                                                                                                                    >  7:15702/1‑149 (MQ=255)
  cGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGc                                                                                                                                    <  5:191155/149‑1 (MQ=255)
   gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGcc                                                                                                                                   >  1:390897/1‑149 (MQ=255)
    ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCg                                                                                                                                  <  1:251250/149‑1 (MQ=255)
     gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGc                                                                                                                                         <  8:111127/141‑1 (MQ=255)
     gCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGc                                                                                                                                         >  7:111127/1‑141 (MQ=255)
      ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGgt                                                                                                                                <  6:292568/149‑1 (MQ=255)
       cGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGgtg                                                                                                                               <  4:9330/149‑1 (MQ=255)
        gCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCGTTGCCGgtgt                                                                                                                              >  2:299499/1‑149 (MQ=255)
         cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGACAC‑C‑GGCAAAGCCGgtgcc                                                                                                                             <  1:302228/149‑3 (MQ=38)
              cGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGttcgcg                                                                                                                         <  2:390897/148‑5 (MQ=255)
               gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGttcgcgg                                                                                                                        <  8:15702/148‑6 (MQ=255)
                     gCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTAGCGGG‑CAcg‑cc                                                                                                               >  8:4744/1‑149 (MQ=37)
                                tGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCTCGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCcaca                                                                                                     <  1:299499/149‑1 (MQ=21)
                                     cGGCCATTCAGAGCCTCAACGAAGCCTTGGCACCACAGCAGGACTGTGGCTGCCTAACGGTGGGAGGGGGCGTGCCCGCGAAACAGGCGG‑CGCGGTGAATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGa                                                                                               <  1:274379/149‑1 (MQ=21)
                                       gCCATTCAGCGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGAtt                                                                                             <  3:31486/149‑1 (MQ=21)
                                              agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGtt                                                                                                                             <  1:343720/112‑1 (MQ=255)
                                              agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGtt                                                                                                                             >  2:343720/1‑112 (MQ=255)
                                              agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑ACTCCCACAGACTGATTGCCAtt                                                                                       <  3:330634/148‑1 (MQ=21)
                                                     gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCGCGAAACAGCCGC‑CCCGGTGCATGGCAC‑C‑GGCTTTGACGGTGTTCGCGGG‑CACGCCC‑GCTCCCACACACAGATTGCCAATAACCGGGc                                                                               >  5:336191/1‑149 (MQ=18)
                                                     gAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACAGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCGTTGCCCGTGTTAGCGGG‑CACGCCC‑GCTCCCACCGACTGATTGCCCTTTAGCGGGc                                                                               >  1:273026/1‑149 (MQ=21)
                                                          aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                          <  1:310546/149‑1 (MQ=21)
                                                                 ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTAGCGGG‑GACGCCC‑GCTCCCCCAGACTGATTGCCATTAAGCGGGCCTATCGCTACAg                                                                   >  8:285146/1‑149 (MQ=25)
                                                                    accacAGCAGGACTTAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                <  3:40447/149‑1 (MQ=21)
                                                                    accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAAAACGCGC‑CCCGGTGCCTTGCAC‑C‑CGCTTTTGCGGTGTTCGCGGGCCCCC‑CC‑GCTCCCCCAGACTGATTGGCCTTAATCGGGCCAATCGGTACAGTgg                                                                >  8:232790/1‑148 (MQ=16)
                                                                              gccTGAGGCTCCCAAAGGTTGGGAGCGGGCGTGCCCGCTAAACAGGCGC‑CGCGGTGCATGGAAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGa                                                                                               <  1:120108/106‑1 (MQ=14)
                                                                              gACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGGATGGCAC‑C‑CGCATTGCCGGTGTTGGCGGG‑CCCGCCC‑CCGCCCACAGACTGATTGCCAATAAGCGGGGCTATAGCGACAGTTACGGCCCGgtt                                                      >  4:84180/1‑149 (MQ=18)
                                                                                cTGCGGCTGCAAAAGGGTGGGAGGGGGCGTGCCCGCTAAACCGGCGC‑CGCGG‑GGATGGTCACC‑GGTTTTGCGGGTGTTCGCGGG‑AACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGttttt                                                    <  8:241208/149‑1 (MQ=255)
                                                                                          cAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCCC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CCCGCCC‑GCGCCCCCAGACTGCTTGCCATTAAGCGGGCCTATGGCGACAGTGACGGCCCGGTTTTTGTTGCcgcg                                          >  3:61986/1‑149 (MQ=18)
                                                                                          cAAACGGTGGGAGCGGGCGTGACCGCGAAACAGGCGC‑CGCCGTGCATGGCAC‑C‑CGCGTTGCCCGTGTTACCGCG‑CACGCCC‑GCTCCCACAGACTGAATGCCATTAAGCGGGGCTAAGGCTACAGTGAAGGCCCGTTTTTTGTTGCCACg                                          >  8:106976/1‑149 (MQ=255)
                                                                                             aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                                <  3:284964/124‑1 (MQ=18)
                                                                                             aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCCGTGCATGGCAC‑C‑CGCATTGCCCGTGTTAGCCGG‑CACGCCC‑CCGCCCCCAGACAGATTGCCATTAAGCGGGCCGATAGCGACAGTgg                                                                >  4:284964/1‑123 (MQ=12)
                                                                                                     aGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                               >  5:125995/1‑149 (MQ=25)
                                                                                                     aGCGGGCGTCCCGGCTAACCCGGCGC‑CGCGGTGCA‑GGCCA‑CGGGCTTTGCCGGTGTT‑CGGGGCCACG‑CCGGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                               <  6:220837/149‑1 (MQ=255)
                                                                                                      gCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCGGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                              <  3:84180/149‑1 (MQ=25)
                                                                                                      gCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                              <  7:4744/149‑1 (MQ=25)
                                                                                                             tGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGCCGGCCCGTTTTTTGTTGCCTCGCTAACGACTCATGACTccc                       >  4:209523/1‑149 (MQ=25)
                                                                                                              gCCCGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                      >  1:260380/1‑149 (MQ=37)
                                                                                                               cccGCGAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                      >  8:104788/1‑148 (MQ=37)
                                                                                                                  gcgAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                   >  8:297421/1‑148 (MQ=37)
                                                                                                                    gAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                       >  5:223380/1‑142 (MQ=37)
                                                                                                                    gAAACAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTccc                       <  6:223380/142‑1 (MQ=37)
                                                                                                                     aaaCAGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg               <  6:125995/149‑1 (MQ=37)
                                                                                                                         aGGTGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                >  6:307118/1‑112 (MQ=21)
                                                                                                                         aGGCGC‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                                <  5:307118/112‑1 (MQ=21)
                                                                                                                           gcgc‑CGCGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc          >  7:287077/1‑148 (MQ=37)
                                                                                                                             gc‑cgcGGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg       >  8:73905/1‑149 (MQ=37)
                                                                                                                                  cgGTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg       >  7:186512/1‑145 (MQ=37)
                                                                                                                                   ggTGCATGGCAC‑C‑GGCTTTGCCGGTGTTCGCGGG‑CACGCCC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAt  <  7:104788/149‑1 (MQ=38)
                                                                                                                                             ||                                                                                                                                             
GACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGC‑CGCGGTGCATGGCCA‑C‑GGCTTTGCCGGTGTT‑TCGGG‑CACG‑CC‑GCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAGAT  >  NC_002947/1499340‑1499616

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: ‑