Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,831,777 A→G 8.1% F116S (TTC→TCC)  PP_3384 ← gluconate 2‑dehydrogenase gamma subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,831,7770AG 8.1% 136.8 / 2.5 50F116S (TTC→TCC) PP_3384gluconate 2‑dehydrogenase gamma subunit
Reads supporting (aligned to +/- strand):  ref base A (21/24);  new base G (2/2);  total (24/26)
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.13e-01

GGTCCCGGGTAGCGCTATCTTGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATGCGCGACACGGCAGCCTGGACAAAAGCCCA  >  NC_002947/3831631‑3831925
                                                                                                                                                  |                                                                                                                                                    
ggTCCCGGGTAGCGCTATCTTGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAAc                                                                                                                                                    >  3:134290/1‑149 (MQ=255)
    ccGGGTAGCGCTATCTTGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGAcc                                                                                                                                                <  8:41038/149‑1 (MQ=255)
      gggTAGCGCTATCTTGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCt                                                                                                                                              >  7:183795/1‑149 (MQ=255)
           gcgcTATCTTGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCAt                                                                                                                                         <  5:50493/149‑1 (MQ=255)
                 tCTTGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCa                                                                                                                                   >  1:290751/1‑149 (MQ=255)
                   ttGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTAc                                                                                                                                     <  8:37085/145‑1 (MQ=255)
                   ttGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTAc                                                                                                                                     >  7:37085/1‑145 (MQ=255)
                      cgcAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACtgt                                                                                                                               <  3:221889/148‑1 (MQ=255)
                         agcaAAATCCTGTCGAGCATTGGGCTTGCACCAGCCATGCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTTTgcgc                                                                                                                           <  3:2999/149‑1 (MQ=255)
                         agcaAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTgcgc                                                                                                                           <  6:224575/149‑1 (MQ=255)
                               aaCCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTcg                                                                                                                     <  6:81254/149‑1 (MQ=255)
                                aCCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTcgc                                                                                                                    <  6:197958/149‑1 (MQ=255)
                                       caccaTTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAgggg                                                                                                             <  2:290751/149‑1 (MQ=255)
                                       caccaTTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAggg                                                                                                              <  3:148478/148‑1 (MQ=255)
                                         ccaTTGGCCTTGCACCAGCCAACCACGGCCGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCCGATCGGGCGGGGAGTCGGCGTGGCACGGCCCGTGCATGGACACCTGTGCGCCAGTCGAGTAggggt                                                                                                            >  4:182233/1‑147 (MQ=255)
                                               gCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCa                                                                                                      >  6:269293/1‑148 (MQ=255)
                                                 cTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGGTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCt                                                                                                   >  6:89650/1‑149 (MQ=255)
                                                       ccAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAgg                                                                                                               <  6:217453/131‑1 (MQ=255)
                                                       ccAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAgg                                                                                                               >  5:217453/1‑131 (MQ=255)
                                                         aGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGGACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTc                                                                                            >  5:245558/1‑148 (MQ=255)
                                                         aGCCATCCACGGCAGCGATGACCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGGACGGCCCTTGCATGTACCACTGTGCGCAAGTCGCGTAGGGGGTATTCATCTGGCgg                                                                                              >  5:241585/1‑146 (MQ=255)
                                                          gCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGa                                                                                          >  3:184043/1‑149 (MQ=255)
                                                             aTCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATg                                                                                        >  3:84715/1‑148 (MQ=255)
                                                               ccACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGa                                                                                          >  2:130307/1‑144 (MQ=255)
                                                               ccACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGa                                                                                          <  1:130307/144‑1 (MQ=255)
                                                                     cAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTa                                                                                                                                      >  8:203135/1‑94 (MQ=255)
                                                                     cAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTa                                                                                                                                      <  7:203135/94‑1 (MQ=255)
                                                                     cAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTcggcg                                                                               <  5:269293/149‑1 (MQ=255)
                                                                      aGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATACATCTGGCGGTCGATGAACTcggcgg                                                                              <  4:134290/149‑1 (MQ=255)
                                                                             cccAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGc                                                                       >  8:20844/1‑149 (MQ=255)
                                                                                     gTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTcggc                                                                                >  4:137501/1‑132 (MQ=255)
                                                                                     gTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTcggc                                                                                <  3:137501/132‑1 (MQ=255)
                                                                                        gATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGccc                                                            <  8:183795/149‑1 (MQ=255)
                                                                                          tttGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCt                                                          >  4:251678/1‑149 (MQ=255)
                                                                                              ttgCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTCGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGc                                                      >  3:274402/1‑149 (MQ=255)
                                                                                                 cGGGCTGAGCTGCAGCTGGAAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGcc                                                     >  2:293630/1‑147 (MQ=255)
                                                                                                  gggCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGGACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAg                                                  <  6:241585/149‑1 (MQ=255)
                                                                                                  gggCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGGACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAg                                                  <  6:245558/149‑1 (MQ=255)
                                                                                                         gCTGCAGCTGGTAGCCCAGTTCTGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCTATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATc                                            <  1:19881/148‑1 (MQ=255)
                                                                                                         gCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCg                                           >  7:240119/1‑149 (MQ=255)
                                                                                                             cAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGc                                         <  7:20844/147‑1 (MQ=255)
                                                                                                                cTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGga                                    <  1:232961/149‑1 (MQ=255)
                                                                                                                   gTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGgatg                                  <  5:144176/148‑1 (MQ=255)
                                                                                                                      gCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCg                                       >  5:213538/1‑140 (MQ=255)
                                                                                                                       cccAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATgcg                             >  2:193664/1‑149 (MQ=255)
                                                                                                                        ccAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATgcgc                            >  7:275333/1‑149 (MQ=255)
                                                                                                                        ccAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATgcgc                            <  8:275333/149‑1 (MQ=255)
                                                                                                                                   ggggcgTCGGCGTTGAACGGCCCTCGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATGCGCGACACGGCAGc                 <  4:274402/149‑1 (MQ=255)
                                                                                                                                        gTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGAGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATGCGCGACACGGCAGCCTGGa            <  2:302147/149‑1 (MQ=255)
                                                                                                                                         tCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATGCGCGACACGGCAGCCTGGa            >  7:36325/1‑148 (MQ=255)
                                                                                                                                                  aaCGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATGCGCGACACGGCAGCCTGGACAAAAGCCCa  <  1:288069/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGTCCCGGGTAGCGCTATCTTGCGCAGCAAAAACCTGTCCACCATTGGCCTTGCACCAGCCATCCACGGCAGCGATGCCCAGGCGGTAGATTTGTTGCGGGCTGAGCTGCAGCTGGTAGCCCAGTTCGGGCGGGGCGTCGGCGTTGAACGGCCCTTGCATGTACCACTGTGCGCCAGTCGCGTAGGGGGTATTCATCTGGCGGTCGATGAACTCGGCGGCACCCGCTTCGAGTGCCCCTGGGCCAAGCTCATCGGCCGGGATGATGCGCGACACGGCAGCCTGGACAAAAGCCCA  >  NC_002947/3831631‑3831925

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