Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,583,110 G→A 7.3% G123G (GGC→GGT)  nudK ← GDP‑mannose pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,583,1100GA 7.3% 171.9 / 2.1 55G123G (GGC→GGT) nudKGDP‑mannose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (23/27);  new base A (2/2);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCAGCTGCCCGCTTTC  >  NC_000913/2582969‑2583253
                                                                                                                                             |                                                                                                                                               
gCCGGTTTTGATCATCTACAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAACGTGGATTAGCTCAGTCACACCGCCTgg                                                                                                                                            <  4:220723/147‑1 (MQ=255)
     ttttGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACat                                                                                                                                      <  2:77032/148‑1 (MQ=255)
       ttGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACata                                                                                                                                     >  6:17432/1‑147 (MQ=255)
          atcatcTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATAc                                                                                                                                  <  4:174365/147‑1 (MQ=255)
            catcTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGt                                                                                                                               >  8:234635/1‑148 (MQ=255)
                    gcgcCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAaataa                                                                                                                       <  5:17432/148‑1 (MQ=255)
                          ggCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTgcg                                                                                                                 >  6:379106/1‑148 (MQ=255)
                           gCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTTTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTgcgc                                                                                                                <  8:122267/148‑1 (MQ=255)
                           gCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTgcgc                                                                                                                <  3:38832/148‑1 (MQ=255)
                           gCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCACCTGGCGACATATACAGTTCAAATAATTTgcg                                                                                                                 >  4:23022/1‑147 (MQ=255)
                            cTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCACCTGGCGACATATACAGTTCAAATAATTTGCGCa                                                                                                               <  3:23022/148‑1 (MQ=255)
                             tGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCAc                                                                                                              >  8:160048/1‑148 (MQ=255)
                                       cTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACt                                                                                                    >  7:347774/1‑148 (MQ=255)
                                          gAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTc                                                                                                  <  7:160048/147‑1 (MQ=255)
                                          gAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCa                                                                                                 <  8:38017/148‑1 (MQ=255)
                                                 tCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATAc                                                                                                                                  <  3:338527/108‑1 (MQ=255)
                                                 tCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATAc                                                                                                                                  >  4:338527/1‑108 (MQ=255)
                                                 tCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTgc                                                                                                                  <  1:103481/124‑1 (MQ=255)
                                                 tCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTgc                                                                                                                  >  2:103481/1‑124 (MQ=255)
                                                      atcttcatcttcGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGtctcc                                                                                     >  5:224041/1‑147 (MQ=255)
                                                       tcttcatcttcGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCtt                                                                                    >  8:293470/1‑148 (MQ=255)
                                                       tcttcatcttcGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCtt                                                                                    >  4:29073/1‑148 (MQ=255)
                                                          tcatcttcGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGAGACATATACAGTTCAAATAATTTGCTAACTTCACCCACTTCATAGCCCGTCTCTTCa                                                                                  >  1:341467/1‑147 (MQ=255)
                                                                 cGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGctt                                                                           <  8:347774/147‑1 (MQ=255)
                                                                       cGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCa                                                                                                 >  2:336523/1‑119 (MQ=255)
                                                                       cGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCa                                                                                                 <  1:336523/119‑1 (MQ=255)
                                                                        gcaaCCGGCGTTAGCGATCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCTCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTATTTGCGa                                                                   <  2:341467/145‑1 (MQ=255)
                                                                        gCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCCTCAATCGCTTCTTTGCGa                                                                   <  6:224041/148‑1 (MQ=255)
                                                                           aCCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATa                                                                <  3:29073/148‑1 (MQ=255)
                                                                              ggCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATacac                                                             <  5:379106/148‑1 (MQ=255)
                                                                                   tAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTc                                                        >  4:132641/1‑148 (MQ=255)
                                                                                   tAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCAACCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTc                                                        >  2:149798/1‑148 (MQ=255)
                                                                                     gcgcgcTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCg                                                      <  7:293470/148‑1 (MQ=255)
                                                                                          cTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTtcg                                                 >  7:140031/1‑148 (MQ=255)
                                                                                            gATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCACCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTg                                                                      >  5:143634/1‑125 (MQ=255)
                                                                                            gATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCACCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTg                                                                      <  6:143634/125‑1 (MQ=255)
                                                                                                 tCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCg                                                                    <  1:85316/122‑1 (MQ=255)
                                                                                                      gTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCcag                                      >  4:58340/1‑147 (MQ=255)
                                                                                                       tATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCcagc                                     <  1:149798/147‑1 (MQ=255)
                                                                                                        aTTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCcagcag                                   >  2:283765/1‑148 (MQ=255)
                                                                                                          tCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAAt                                                                 >  2:25897/1‑116 (MQ=255)
                                                                                                          tCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAAt                                                                 <  1:25897/116‑1 (MQ=255)
                                                                                                          tCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGcc                                 >  8:24519/1‑148 (MQ=255)
                                                                                                                aTAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGgcg                            >  7:387743/1‑147 (MQ=255)
                                                                                                                        gTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCa                   <  7:234635/148‑1 (MQ=255)
                                                                                                                            aTTAGCTCAGTCACACCGCCTGGCTACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCa               <  8:387743/148‑1 (MQ=255)
                                                                                                                            aTTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCa               <  7:24519/148‑1 (MQ=255)
                                                                                                                                      tCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATcc                                        <  7:206379/113‑1 (MQ=255)
                                                                                                                                       cacaCCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTc                    <  2:97939/132‑1 (MQ=255)
                                                                                                                                       cacaCCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTc                    >  1:97939/1‑132 (MQ=255)
                                                                                                                                       cacaCCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCAGCTGCCCGCt     >  7:173948/1‑147 (MQ=255)
                                                                                                                                        acacCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCAGCTGc          >  4:174567/1‑141 (MQ=255)
                                                                                                                                        acacCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCAGCTGc          <  3:174567/141‑1 (MQ=255)
                                                                                                                                        acacCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCAGCTGCCCGCttt   <  8:140031/148‑1 (MQ=255)
                                                                                                                                         cacCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCAGCTGCCCGCTTTc  <  3:58340/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCTGGCGACATATACAGTTCAAATAATTTGCGCACTTCACCCACTTCATAGCCCGTCTCTTCAATCGCTTCTTTGCGAATACACACTTCCGGTTCGTCGTTATCCAGCAGCCCGGCGCAGCTTTCAATCAGCTGCCCGCTTTC  >  NC_000913/2582969‑2583253

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