Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,591,637 G→T 11.6% A82D (GCC→GAC)  nudF ← ADP‑ribose/sugar pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,591,6370GT11.6% 118.9 / 6.9 43A82D (GCC→GAC) nudFADP‑ribose/sugar pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (22/16);  new base T (3/2);  total (25/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 = 6.49e-01

GAAGTAGCGGGTCATCGGCCACAGGGCGCTGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAGCTCGTGGCGCAGGCGTAAC  >  NC_002947/5591489‑5591783
                                                                                                                                                    |                                                                                                                                                  
gAAGTAGCGGGTCATCGGCCACAGGGCGCTGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGAAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAggg                                                                                                                                                    >  2:51830/1‑149 (MQ=255)
 aaGTAGCGGGTCATCGGCCACAGGGCGCTGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGgc                                                                                                                                                   <  2:259444/149‑1 (MQ=255)
      gCGGGTCATCGGCCACAGGGCGCTGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCa                                                                                                                                              >  7:211406/1‑149 (MQ=255)
      gCGGGTCATCGGCCACAGGGCGCTGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCACGGGTTGGCGACCTTGTCCAGGGCGCCCa                                                                                                                                              >  1:137374/1‑149 (MQ=255)
       cGGGTCATCGGCCACAGGGCGCTGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGCCCTGGACAAGGTCGccaa                                                                                                                                              <  4:231528/148‑3 (MQ=255)
         ggTCATCGGCCACAGGGCGCTGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCAcgc                                                                                                                                           <  6:73932/149‑1 (MQ=255)
                             tGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGcac                                                                                                                        >  8:278155/1‑148 (MQ=255)
                                                ttcttcAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGt                                                                                                    >  7:77791/1‑149 (MQ=255)
                                                          tCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGc                                                                                                                          <  1:285303/117‑1 (MQ=255)
                                                          tCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGc                                                                                                                          >  2:285303/1‑117 (MQ=255)
                                                                        ttCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAggg                                                                                                                                                    <  4:95184/77‑1 (MQ=255)
                                                                        ttCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAggg                                                                                                                                                    >  3:95184/1‑77 (MQ=255)
                                                                           cTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGt                                                                         >  4:199897/1‑149 (MQ=255)
                                                                                  tgctCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGca                                                                   >  4:182262/1‑148 (MQ=255)
                                                                                    ctCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGcaca                                                                 >  4:301646/1‑148 (MQ=255)
                                                                                        tCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAg                                                            >  1:127428/1‑149 (MQ=255)
                                                                                        tCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCACCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAg                                                            >  3:22598/1‑149 (MQ=255)
                                                                                              tCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTcgc                                                       >  6:24870/1‑148 (MQ=255)
                                                                                                gATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGAc                                                    >  7:158252/1‑149 (MQ=255)
                                                                                                 aTCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACt                                                   <  2:127428/149‑1 (MQ=255)
                                                                                                 aTCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCACCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACt                                                   <  4:22598/149‑1 (MQ=255)
                                                                                                       ccGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGTGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGAtctc                                             >  5:230477/1‑149 (MQ=255)
                                                                                                          gCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTcgcg                                          <  8:158252/149‑1 (MQ=255)
                                                                                                           cAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCgcgc                                         <  4:91499/149‑1 (MQ=255)
                                                                                                            aaCCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGcc                                        <  3:199897/149‑1 (MQ=255)
                                                                                                              ccATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGccc                                       >  4:197547/1‑148 (MQ=255)
                                                                                                                  ctcGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGc                                  >  3:21941/1‑149 (MQ=255)
                                                                                                                     gATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTcg                                            >  2:147596/1‑136 (MQ=255)
                                                                                                                     gATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTcg                                            <  1:147596/136‑1 (MQ=255)
                                                                                                                     gATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGcccc                               <  4:21941/149‑1 (MQ=255)
                                                                                                                      aTCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGccccc                              <  7:278155/149‑1 (MQ=255)
                                                                                                                        cAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAACACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCgg                            >  1:104333/1‑149 (MQ=255)
                                                                                                                         aGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCgg                            <  8:211406/148‑1 (MQ=255)
                                                                                                                          gCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGc                           >  5:135026/1‑148 (MQ=255)
                                                                                                                              ggggTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCaaa                        >  2:293135/1‑147 (MQ=255)
                                                                                                                               gggTTGGCGACCTTGTCCAGGGTGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACACTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAg                     <  6:230477/149‑1 (MQ=255)
                                                                                                                                ggTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCgg                            <  5:167714/141‑1 (MQ=255)
                                                                                                                                ggTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCgg                            >  6:167714/1‑141 (MQ=255)
                                                                                                                                 gtgggCGCCCTGGACAAGGTCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAGCt                   >  1:2923/4‑149 (MQ=255)
                                                                                                                                    ggCGCCCTGGACAAGGTCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAGCTCGt                >  1:126728/1‑149 (MQ=255)
                                                                                                                                    ggCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCCCGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAGCTCGt                >  5:67639/1‑149 (MQ=255)
                                                                                                                                     gCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTgg                                                                                                      >  4:23260/1‑62 (MQ=255)
                                                                                                                                     gCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTgg                                                                                                      <  3:23260/62‑1 (MQ=255)
                                                                                                                                     gCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAGCTCGTg               <  8:77791/149‑1 (MQ=255)
                                                                                                                                      cGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAGCTCGTgg              <  5:24870/149‑1 (MQ=255)
                                                                                                                                            ggacaaGGTCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTcgc                                                       <  6:266034/97‑1 (MQ=255)
                                                                                                                                            ggacaaGGTCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTcgc                                                       >  5:266034/6‑102 (MQ=255)
                                                                                                                                                  gggCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACAGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAGCTCGTGGCGCAGGCGTAAc  <  6:195448/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAAGTAGCGGGTCATCGGCCACAGGGCGCTGAAGGTCAGGCCGGCTTCTTCTTCAGCCTCGCGATGGGCGACTTCCTCGGGCTGCTCGTCTTTGTCGATCAGCCCGGCAACCATCTCGATCAGCCAGGGGTTGGCGACCTTGTCCAGGGCGCCCACGCGGAACTGTTCGATCAGCACCACTTCATCGCGCGCTGGGTCGTACGGCAGCACGCATACCGCATCGTGGCGCACAAACAGTTCGCGACTGATCTCGCGCCCCATGCCCCCGGCAAATAGCTCGTGGCGCAGGCGTAAC  >  NC_002947/5591489‑5591783

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