Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,328,147 G→C 7.2% A20G (GCC→GGC)  udg ← UDP‑glucose 6‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,328,1470GC 7.2% 164.0 / 4.2 56A20G (GCC→GGC) udgUDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base C (2/2);  total (28/28)
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.18e-01

GTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATG  >  NC_002947/3328003‑3328294
                                                                                                                                                |                                                                                                                                                   
gtggtgAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAg                                                                                                                                                 <  6:106174/149‑1 (MQ=255)
gtggtgAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAg                                                                                                                                                 <  7:114552/149‑1 (MQ=255)
 tggtgAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAgg                                                                                                                                                <  6:15592/149‑1 (MQ=255)
    tgAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGcac                                                                                                                                             <  6:312720/149‑1 (MQ=255)
        gcgcAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGACCCAGGcagtgtg                                                                                                                                         >  3:179943/1‑144 (MQ=255)
           cAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCAATGCACAACACCGAATGGCCTACCTCGGCCAGGCACACTGCCt                                                                                                                                      >  4:129986/1‑149 (MQ=255)
             gCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTgg                                                                                                                                    >  5:345036/1‑149 (MQ=255)
             gCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTgg                                                                                                                                    <  8:27248/149‑1 (MQ=255)
                gACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGt                                                                                                                                  >  2:292131/1‑148 (MQ=255)
                 aCCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGACCCAGGCAGTGTGCCTggccca                                                                                                                                <  4:179943/149‑5 (MQ=255)
                        gCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGa                                                                                                                          >  6:278710/1‑148 (MQ=255)
                                cttctCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCgg                                                                                                                 >  7:68798/1‑149 (MQ=255)
                                  tctcCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGtt                                                                                                               >  7:158656/1‑149 (MQ=255)
                                  tctcCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGtt                                                                                                               >  1:142277/1‑149 (MQ=255)
                                       agcagcGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCg                                                                                                            <  2:163364/147‑1 (MQ=255)
                                       agcagcGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCg                                                                                                            >  1:163364/1‑147 (MQ=255)
                                          agcGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGCCCAGGCAGTGTGCCTggc                                                                                                                                   >  4:182177/1‑120 (MQ=255)
                                          agcGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGCCCAGGCAGTGTGCCTggc                                                                                                                                   <  3:182177/121‑2 (MQ=255)
                                             ggcgCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGt                                                                                                    <  5:278710/149‑1 (MQ=255)
                                               cgcAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGTCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGa                                                                                                  <  3:189482/149‑1 (MQ=255)
                                                gcAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGAc                                                                                                 >  1:117137/1‑149 (MQ=255)
                                                gcAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGAc                                                                                                 <  2:334911/149‑1 (MQ=255)
                                                     gCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCa                                                                                            >  5:60089/1‑149 (MQ=255)
                                                           cTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGAc                                                                                      <  2:142277/149‑1 (MQ=255)
                                                                     gggCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGa                                                                            >  1:116255/1‑149 (MQ=255)
                                                                               cTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCt                                                                  <  8:68798/149‑1 (MQ=255)
                                                                                   aCTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTg                                                                                                   >  2:266229/1‑112 (MQ=255)
                                                                                   aCTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTg                                                                                                   <  1:266229/112‑1 (MQ=255)
                                                                                     tCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCTAATGGCCTACCTGTGCCAGGAACACTGCCTTGGTAAGGCAGACATAGCCGGTTCCGAAAACAGTGACCTTAATTTTACAGCTCCTTGATAGGAAATCTTaaaa                                                             >  6:270484/1‑148 (MQ=255)
                                                                                     tCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGaaaa                                                             >  6:270486/1‑148 (MQ=255)
                                                                                     tCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGaaaa                                                             >  1:20084/1‑148 (MQ=255)
                                                                                             gCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGt                                                      <  3:129986/147‑1 (MQ=255)
                                                                                             gCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTca                                                    <  6:52888/149‑1 (MQ=255)
                                                                                              ccACCCGATCGACATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTca                                                    >  8:268388/1‑148 (MQ=255)
                                                                                                 cccGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATc                                                >  6:96528/1‑149 (MQ=255)
                                                                                                  ccGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCg                                               >  6:249024/1‑149 (MQ=255)
                                                                                                   cGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCgt                                              >  4:179640/1‑149 (MQ=255)
                                                                                                      tCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCgttgt                                           >  8:235095/1‑149 (MQ=255)
                                                                                                         gCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGCAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCgttgtt                                          >  4:200941/1‑147 (MQ=255)
                                                                                                              gACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAAt                                   >  5:127182/1‑149 (MQ=255)
                                                                                                               aCATCCATGCACACCACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGTGGCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAAtt                                  <  8:77405/149‑1 (MQ=255)
                                                                                                               aCATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAAtt                                  <  4:7848/149‑1 (MQ=255)
                                                                                                               aCATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAAtt                                  <  6:60089/149‑1 (MQ=255)
                                                                                                                 aTCCTTGCACAACTCCGAATGTCCTACCTGGGCAAGGGACTCTGCCTTGGTCAGGCCGACATAGCCGGTTCCGAATACCGTGACCTTCATTTGAGATCTCCTTGATAGGAAAGCTGAATTTAAATGTCACATCGTTGTTTAGGAATtgt                                <  5:122596/149‑1 (MQ=255)
                                                                                                                    cATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGAc                             >  3:70182/1‑149 (MQ=255)
                                                                                                                        cacaACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGa                         <  3:99107/149‑1 (MQ=255)
                                                                                                                         acaacaCCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAg                        <  7:235095/149‑1 (MQ=255)
                                                                                                                         acaacaCCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAg                        <  1:68552/149‑1 (MQ=255)
                                                                                                                              acCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGaaa                   <  1:292131/149‑1 (MQ=255)
                                                                                                                               ccGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAAt                  <  3:179640/149‑1 (MQ=255)
                                                                                                                                      gCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTg           <  5:73141/149‑1 (MQ=255)
                                                                                                                                       ccTACCTGGGCCAGGCACACTGCCTGGGTCTGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAg                                                                    >  1:28874/1‑91 (MQ=255)
                                                                                                                                       ccTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAg                                                                    <  2:28874/91‑1 (MQ=255)
                                                                                                                                       ccTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTg           >  4:29257/1‑148 (MQ=255)
                                                                                                                                       ccTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTg           <  3:29257/148‑1 (MQ=255)
                                                                                                                                         tACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGc        >  4:40956/1‑149 (MQ=255)
                                                                                                                                               ggCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATg  >  4:162783/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATG  >  NC_002947/3328003‑3328294

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