Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,831,436 G→A 7.8% R230C (CGC→TGC)  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,4360GA 7.8% 141.8 / 2.0 51R230C (CGC→TGC) PP_3384gluconate 2‑dehydrogenase gamma subunit
Reads supporting (aligned to +/- strand):  ref base G (25/22);  new base A (2/2);  total (27/24)
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.99e-01

CCCGGCCTCGGTCAGCTCTTTGGCCATGATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGC  >  NC_002947/3831290‑3831580
                                                                                                                                                  |                                                                                                                                                
cccGGCCTCGGTCAGCTCTTTGGCCATGATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCt                                                                                                                                                >  5:326708/1‑149 (MQ=255)
  cGGCCTCGGTCAGCTCTTTGGCCATGATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTc                                                                                                                                               >  3:292778/1‑148 (MQ=255)
   ggCCTCGGTCAGCTCTTTGGCCATGATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGa                                                                                                                                             <  8:25153/149‑1 (MQ=255)
     ccTCGGTCAGCTCTTTGGCCATGATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGAcc                                                                                                                                           <  4:292778/149‑1 (MQ=255)
                  tttGGCCATGATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGaa                                                                                                                                <  2:74146/147‑1 (MQ=255)
                        cATGATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGgcgc                                                                                                                        >  6:240850/1‑149 (MQ=255)
                           gATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCggg                                                                                                                     >  8:12468/1‑149 (MQ=255)
                            aTCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGgc                                                                                                                    >  1:130135/1‑149 (MQ=255)
                              cGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGgcgc                                                                                                                  <  2:113486/149‑1 (MQ=255)
                                  ccGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCggg                                                                                                              <  8:217438/149‑1 (MQ=255)
                                        cAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGcc                                                                                                        <  5:240850/149‑1 (MQ=255)
                                           ccgaagccgaCGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGAt                                                                                                     >  4:98110/1‑149 (MQ=255)
                                                   gacgaTCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTc                                                                                              >  7:81022/1‑148 (MQ=255)
                                                        tCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGACCGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGACGCCGATCTGAGTCCAGCCg                                                                                        >  2:265976/1‑149 (MQ=255)
                                                              cATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCAtt                                                                                  <  7:8151/149‑1 (MQ=255)
                                                                  cACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCt                                                                               <  6:326708/148‑1 (MQ=255)
                                                                   aCTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCttt                                                                             <  8:81022/149‑1 (MQ=255)
                                                                    cTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTttg                                                                            >  3:21487/1‑149 (MQ=255)
                                                                    cTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTttg                                                                            >  3:240100/1‑149 (MQ=255)
                                                                     tttGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCTGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTttgt                                                                           >  8:256826/1‑149 (MQ=255)
                                                                         ttCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCttt                                                                             >  1:282408/1‑143 (MQ=255)
                                                                         ttCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCttt                                                                             <  2:282408/143‑1 (MQ=255)
                                                                           cAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCACTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAg                                                                                           >  3:110819/1‑127 (MQ=255)
                                                                           cAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCACTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAg                                                                                           >  3:106151/1‑127 (MQ=255)
                                                                           cAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCACTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAg                                                                                           <  4:106151/127‑1 (MQ=255)
                                                                           cAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCACTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAg                                                                                           <  4:110819/127‑1 (MQ=255)
                                                                                 cGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGa                                                                                                 <  7:296765/115‑1 (MQ=255)
                                                                                 cGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGa                                                                                                 >  8:296765/1‑115 (MQ=255)
                                                                                   tCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATc                                                             <  1:47005/149‑1 (MQ=255)
                                                                                            cTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTttgt                                                                           <  8:248791/126‑1 (MQ=255)
                                                                                            cTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTttgt                                                                           >  7:248791/1‑126 (MQ=255)
                                                                                               aGGCCCTCTCACCACGGATCGATACAGCTGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCagaa                                                  <  7:256826/148‑1 (MQ=255)
                                                                                               aGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCAt                                                                  >  5:259784/1‑132 (MQ=255)
                                                                                               aGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCAt                                                                  <  6:259784/132‑1 (MQ=255)
                                                                                                   ccTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAAcc                                             >  8:183578/1‑149 (MQ=255)
                                                                                                       tcACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCtt                                          >  1:301913/1‑148 (MQ=255)
                                                                                                         accacGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCa                                        <  7:12468/148‑1 (MQ=255)
                                                                                                                   gATCCAGCCGGGATCGGGTAAGGCTCATTGCGCTCGACCCTGTCCATGAATTCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGtttt                              <  4:240100/148‑1 (MQ=255)
                                                                                                                        aGCCGGGAACGGGTAAGGCTCATTCCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACca                        >  3:180412/1‑149 (MQ=255)
                                                                                                                          cccggAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGCTCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACcag                       >  6:124817/4‑148 (MQ=255)
                                                                                                                          ccGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACcagc                      <  2:307674/149‑1 (MQ=255)
                                                                                                                            gggAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGCGTCCAGCCGACCAATCCTTTGTTGCCGCCATGCATCGGGTCGCAGAAGAACCCATCACGCGTGTTTTGCACcagcag                    >  5:73286/1‑149 (MQ=255)
                                                                                                                             ggAACGGGAAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCAATCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCCCGCCTCTTTTGCACAAGCAgg                   >  2:132816/1‑149 (MQ=255)
                                                                                                                                  gggTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGc                            <  7:135763/135‑1 (MQ=255)
                                                                                                                                  gggTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGc                            >  8:135763/1‑135 (MQ=255)
                                                                                                                                     tAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTgaagaa           >  8:56584/1‑149 (MQ=255)
                                                                                                                                         gCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCagaa                                                  >  3:241962/1‑106 (MQ=255)
                                                                                                                                         gCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCagaa                                                  <  4:241962/106‑1 (MQ=255)
                                                                                                                                         gCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACAcc       >  3:302924/1‑149 (MQ=255)
                                                                                                                                             aTTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTg    <  3:98110/148‑1 (MQ=255)
                                                                                                                                              ttGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGc  <  7:56584/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CCCGGCCTCGGTCAGCTCTTTGGCCATGATCGCACCGGCCCAGCCGAAGCCGACGATCACCGCATCCACTTTGTTCAACACCGTCGCCATGCCTCAGGCCCTCTCACCACGGATCGATACAGCCGGGAACGGGTAAGGCTCATTGCGCTCGACCCAGTCCATGAAATCGGCGCGGGCGCCGGGGAAGCCGATCTGAGTCCAGCCGACCATTCCTTTGTTGCCGCCATGAATCGGGTCGCAGAAGAACCCTTCACGCGTGTTTTGCACCAGCAGGCTGAAGAACACCTTGGC  >  NC_002947/3831290‑3831580

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