Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,328,156	A→T	11.9%	V17E (GTG→GAG) ‡	udg ←	UDP‑glucose 6‑dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,328,156	0	A	T	11.9%	95.7 / 4.3	40	V17E (GTG→GAG) ‡	udg	UDP‑glucose 6‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (15/20);  new base T (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.36e-01

GAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATGGCCGTT  >  NC_002947/3328008‑3328300                                                                                                                                                     |                                                                                                                                                 gAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGcaca                                                                                                                                                  >  5:318532/1‑149 (MQ=255)           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                                                                            <  3:261223/131‑1 (MQ=255)                                                                                         gCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGaa                                                                            >  4:261223/1‑131 (MQ=255)                                                                                                  tCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCgttgt                                                 <  6:19839/149‑1 (MQ=255)                                                                                                     gCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                              <  3:30537/149‑1 (MQ=255)                                                                                                        tCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGa                                           <  2:45130/149‑1 (MQ=255)                                                                                                          gACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                              <  3:256986/144‑1 (MQ=255)                                                                                                          gACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTa                                              >  4:256986/1‑144 (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)                                                                                                                                                     |                                                                                                                                                 GAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGCGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGGCCGACATAGCCGGTTCCGAAAACCGTGACCTTCATTTGACAGCTCCTTGATAGGAAAGCTGAAAAAAAATGTCACATCGTTGTTTAGGAATTGTGACCTGAGGGAAATGCTCGTGGGCCATATGGCCGTT  >  NC_002947/3328008‑3328300

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