Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,328,023	G→A	7.8%	C61C (TGC→TGT)	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,023	0	G	A	7.8%	143.2 / 2.6	51	C61C (TGC→TGT)	udg	UDP‑glucose 6‑dehydrogenase
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.00e+00

GCGTGTTCGAGGATGCTGTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGG  >  NC_002947/3327877‑3328170                                                                                                                                                   |                                                                                                                                                      gCGTGTTCGAGGATGCTGTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCa                                                                                                                                                     <  1:38308/149‑1 (MQ=255)      ttCGAGGATGCTGTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGc                                                                                                                                                  >  4:295519/1‑147 (MQ=255)            cgatcagTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGttc                                                                                                                                           <  5:240002/142‑1 (MQ=255)             aTGCTGTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAg                                                                                                                                                    >  4:96780/1‑138 (MQ=255)             aTGCTGTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAg                                                                                                                                                    <  3:96780/138‑1 (MQ=255)              tGCTGTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGttcttc                                                                                                                                        <  1:136783/149‑1 (MQ=255)                cTGTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTtcttct                                                                                                                                       <  3:302487/148‑1 (MQ=255)                   tCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGttc                                                                                                                                           >  6:240002/1‑141 (MQ=255)                       ccACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCagcagc                                                                                                                               >  5:171887/1‑149 (MQ=255)                       ccACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCAagcagc                                                                                                                               >  5:197021/1‑149 (MQ=255)                         aCGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGttc                                                                                                                                           >  7:70079/1‑135 (MQ=255)                         aCGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGttc                                                                                                                                           <  8:70079/135‑1 (MQ=255)                         aCGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAgcgg                                                                                                                             <  6:197021/149‑1 (MQ=255)                          cGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAgcggc                                                                                                                            >  8:87405/1‑149 (MQ=255)                                   cGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCACAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCg                                                                                                                 >  5:34276/1‑149 (MQ=255)                                   cGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCACAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCg                                                                                                                 >  5:29754/1‑149 (MQ=255)                                        ttCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGttcttc                                                                                                                                        <  7:97866/123‑1 (MQ=255)                                        ttCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGttcttc                                                                                                                                        >  8:97866/1‑123 (MQ=255)                                          cAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCg                                                                                                           >  3:259100/1‑148 (MQ=255)                                                 gcgcTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCgg                                                                                                   <  5:4984/149‑1 (MQ=255)                                                       ccATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCaa                                                                                               <  5:304647/147‑1 (MQ=255)                                                             gCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCtt                                                                                        >  1:139705/1‑148 (MQ=255)                                                                   ggcggAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGTCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAg                                                                                 >  3:5545/1‑149 (MQ=255)                                                                      ggAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCt                                                                              >  6:93005/1‑149 (MQ=255)                                                                        aGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTgg                                                                            <  7:87405/149‑1 (MQ=255)                                                                          tACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGcc                                                                          >  6:116152/1‑149 (MQ=255)                                                                          tACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGc                                                                           >  2:95766/1‑148 (MQ=255)                                                                             cGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCAcc                                                                       >  4:72252/1‑149 (MQ=255)                                                                                 ggCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGAt                                                                   <  3:295519/149‑1 (MQ=255)                                                                                      tGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGAGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCa                                                              >  1:216122/1‑149 (MQ=255)                                                                                        aaCTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATc                                                            >  8:124009/1‑149 (MQ=255)                                                                                              aaCCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATc                                                      >  3:143436/1‑149 (MQ=255)                                                                                                     aaTAGTGCGCAGCCGCTGCGTCAGCGGTGGAGACGCGCAGACGACAGGAGTCGATGTCCTTTTCCCGCCGCCGCCAGCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATg                                                  <  1:95766/147‑1 (MQ=255)                                                                                                         gTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCacaacc                                           >  8:161519/1‑148 (MQ=255)                                                                                                             gaaGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCg                                        >  5:274133/3‑148 (MQ=255)                                                                                                                 ccGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTg                                                                             <  8:19636/107‑1 (MQ=255)                                                                                                                 ccGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTg                                                                             >  7:19636/1‑107 (MQ=255)                                                                                                                   gCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCcagcag                                                                                                                                >  1:2626/1‑56 (MQ=255)                                                                                                                   gCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCcagcag                                                                                                                                <  2:2626/56‑1 (MQ=255)                                                                                                                     tGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTa                               <  6:274133/149‑1 (MQ=255)                                                                                                                       cGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTAcc                             >  7:178454/1‑149 (MQ=255)                                                                                                                           aTCCGTGGTGAAGCGCAGCCGACCACAGGCGTGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTggg                         <  6:34276/149‑1 (MQ=255)                                                                                                                           aTCCGTGGTGAAGCGCAGCCGACCACAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCAGGCGTTCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTggt                         <  6:29754/149‑2 (MQ=255)                                                                                                                              cGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGcc                       <  5:93005/148‑1 (MQ=255)                                                                                                                                tggtgAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAgg                    <  4:259100/149‑1 (MQ=255)                                                                                                                                  gtgAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGaa                                      >  8:54943/1‑129 (MQ=255)                                                                                                                                  gtgAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGaa                                      <  7:54943/129‑1 (MQ=255)                                                                                                                                        cgcAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGc            <  7:161519/149‑1 (MQ=255)                                                                                                                                         gcAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGcc           >  6:276651/1‑149 (MQ=255)                                                                                                                                                aCCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCa    <  5:116152/149‑1 (MQ=255)                                                                                                                                                  cGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAgg  <  8:178454/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GCGTGTTCGAGGATGCTGTCGACCACGGCGAAAACGTGCTTCAGGTCGGCGCTACCATCGGCTTGCGGCGGAGTACCGACGGCAATGAACTGCAACCGGGCATAGTTGGCAGCCGCTGCGGCATCCGTGGTGAAGCGCAGCCGACCGCAGGCGAGGTTCTTCTCCAGCAGCGGC‑‑GCAAGGCCGGGCTCGAAGATCGGGCAATGGCCTTCACTCAGGCTGGCCACCCGATCGGCATCGACATCCATGCACAACACCGAATGGCCTACCTGGGCCAGGCACACTGCCTGGGTCAGG  >  NC_002947/3327877‑3328170

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