Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,328,157 C→G 12.1% V17L (GTG→CTG) ‡ udg ← UDP‑glucose 6‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,328,1570CG12.1% 110.9 / 5.2 40V17L (GTG→CTG) ‡udgUDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (15/20);  new base G (3/2);  total (18/22)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATGGCCGTTATTCC  >  NC_002947/3328018‑3328305
                                                                                                                                           |                                                                                                                                                    
cGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGt                                                                                                                                             <  2:157794/149‑1 (MQ=255)
 gACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGACCCAGGCAGTGTGCCTggccc                                                                                                                                            <  5:221880/149‑4 (MQ=255)
    cGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACAGCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAgg                                                                                                                                         <  5:156636/149‑1 (MQ=255)
           gAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACa                                                                                                                                   >  7:30470/1‑148 (MQ=255)
             ggTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAg                                                                                                                                >  1:45130/1‑149 (MQ=255)
              gTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCGCTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAg                                                                                                                                <  7:210869/148‑1 (MQ=255)
              gTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCTAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAg                                                                                                                                <  8:30470/148‑1 (MQ=255)
                    ctcCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCgg                                                                                                                            <  1:18156/146‑1 (MQ=255)
                       cagcagCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGc                                                                                                                               >  7:150481/1‑140 (MQ=255)
                       cagcagCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGc                                                                                                                               <  8:150481/140‑1 (MQ=255)
                       cagcagCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTc                                                                                                                         <  8:187220/146‑1 (MQ=255)
                             cggcgCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCg                                                                                                                >  6:192407/1‑149 (MQ=255)
                               gcgcAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGATGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGt                                                                                                               >  2:39732/1‑148 (MQ=255)
                                        cGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCAtt                                                                                                     >  3:188963/1‑149 (MQ=255)
                                         gggCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCAttt                                                                                                    <  5:98001/149‑1 (MQ=255)
                                                       ggCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGa                                                                                       >  4:30537/1‑148 (MQ=255)
                                                                          gCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTc                                                                                            <  4:174109/124‑1 (MQ=255)
                                                                          gCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTc                                                                                            >  3:174109/1‑124 (MQ=255)
                                                                              gCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGaa                                                                                 >  4:261223/1‑131 (MQ=255)
                                                                              gCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGaa                                                                                 <  3:261223/131‑1 (MQ=255)
                                                                                       tCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCgttgt                                                      <  6:19839/149‑1 (MQ=255)
                                                                                          gCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                                   <  3:30537/149‑1 (MQ=255)
                                                                                             tCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGa                                                <  2:45130/149‑1 (MQ=255)
                                                                                               gACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                                   >  4:256986/1‑144 (MQ=255)
                                                                                               gACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                                   <  3:256986/144‑1 (MQ=255)
                                                                                                    ccATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATtgtg                                          >  8:124774/1‑148 (MQ=255)
                                                                                                    ccATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGa                                         >  5:93490/1‑149 (MQ=255)
                                                                                                     cATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGAc                                        <  6:318532/149‑1 (MQ=255)
                                                                                                      aTGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACTGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGAcc                                       >  7:45153/1‑149 (MQ=255)
                                                                                                                       gCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                                   >  1:260815/1‑120 (MQ=255)
                                                                                                                       gCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                                   <  2:260815/120‑1 (MQ=255)
                                                                                                                       gCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTg                      <  1:39732/149‑1 (MQ=255)
                                                                                                                         cTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTggg                    <  4:188963/149‑1 (MQ=255)
                                                                                                                            ccTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTgg                     >  6:165346/1‑145 (MQ=255)
                                                                                                                            ccTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTgg                     <  5:165346/145‑1 (MQ=255)
                                                                                                                              tgacccAGGCAGTGTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                                   >  7:279044/5‑113 (MQ=255)
                                                                                                                              tgacccAGGCAGTGTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                                   <  8:279044/109‑1 (MQ=255)
                                                                                                                               gacccAGGCAGTGTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGGAATGCTCGTGGGCCatat              >  3:224513/4‑149 (MQ=255)
                                                                                                                               gacccAGGCAGTGTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGAGCCCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATAGTGACCTGAGGGAAATGCTCGTGGGCCatat              >  5:290594/4‑149 (MQ=255)
                                                                                                                                      gCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATGGCCGtt       <  6:93490/149‑1 (MQ=255)
                                                                                                                                           cTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATGGCCGTTATTcc  <  5:192407/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATGGCCGTTATTCC  >  NC_002947/3328018‑3328305

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: