Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,236,197	G→A	9.5%	I192I (ATC→ATT)	catA‑I ←	catechol 1,2‑dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,236,197	0	G	A	9.5%	113.6 / 3.2	42	I192I (ATC→ATT)	catA‑I	catechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (19/19);  new base A (2/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CAAAGTTGATCTGCGTGGTCAGGTGGCGGTGCCCCGGTGCCGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCCACAGGTCGACGGTGGCACCGGCCAA  >  NC_002947/4236052‑4236334                                                                                                                                                  |                                                                                                                                               cAAAGTTGATCTGCGTGGTCAGGTGGCGGTGCCCCGGTGCCGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATc                                                                                                                                             >  6:256618/1‑149 (MQ=255)  aaaGTTGATCTGCGTGGTCAGGTGGCGGTGCCCCGGTGCCGAGATGAAGAAGTGCACGTGCGACGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGTGTCGCAGCCATACCCGGCCGGCACGATc                                                                                                                                             >  5:22035/1‑148 (MQ=255)     gTTGATCTGCGTGGTCAGGTGGCGGTGCCCCGGTGCCGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAgc                                                                                                                                         <  1:111536/149‑1 (MQ=255)                  gTCAGGTGGCGGTGCCCCGGTGCCGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCg                                                                                                                         <  2:44962/149‑1 (MQ=255)                                      tGCCGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGtgatga                                                                                                   <  7:176990/149‑1 (MQ=255)                                         cGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGtgatga                                                                                                   <  8:230873/146‑1 (MQ=255)                                         cGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGtgatga                                                                                                   >  7:230873/1‑146 (MQ=255)                                         cGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATAc                                                                                                <  4:19945/149‑1 (MQ=255)                                          gagaTGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACg                                                                                               >  6:229128/1‑149 (MQ=255)                                                aagaagTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACg                                                                                         >  5:70457/1‑149 (MQ=255)                                                aagaagTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACAATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACg                                                                                         >  8:20155/1‑149 (MQ=255)                                                aagaagTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACAATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACg                                                                                         >  8:15410/1‑149 (MQ=255)                                                 agaagTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGtga                                                                                                      <  4:14337/135‑1 (MQ=255)                                                 agaagTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGtga                                                                                                      >  3:14337/1‑135 (MQ=255)                                                        cACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTg                                                                                  <  6:70457/148‑1 (MQ=255)                                                          cGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAAc                                                                               <  1:292611/149‑1 (MQ=255)                                                            tGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCGTGCGGGTCGCAGCCTTACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTc                                                                             <  8:26298/149‑1 (MQ=255)                                                            tGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTTCGGGTAGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTc                                                                             <  8:21646/149‑1 (MQ=255)                                                                          gCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCgg                                                                           >  4:280301/1‑137 (MQ=255)                                                                          gCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCgg                                                                           <  3:280301/137‑1 (MQ=255)                                                                            cGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGtcgaatc                                                              >  5:115628/1‑148 (MQ=255)                                                                                         cAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGCGTGGAGAGCGGCCCTCGGGACATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACg                                                      >  1:95565/1‑148 (MQ=21)                                                                                            gTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTgatgat                                                                                                  >  3:272816/1‑96 (MQ=255)                                                                                            gTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTgatgat                                                                                                  <  4:272816/96‑1 (MQ=255)                                                                                                 gggaTTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATATGTGCCCTgg                                        <  6:22035/146‑1 (MQ=255)                                                                                                      tCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGtgatg                                                                                                    <  6:28748/84‑1 (MQ=255)                                                                                                      tCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGtgatg                                                                                                    >  5:28748/1‑84 (MQ=255)                                                                                                       ccTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTg                                  >  2:223486/1‑149 (MQ=255)                                                                                                        cTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTgg                                 >  8:193103/1‑149 (MQ=255)                                                                                                         tGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGc                                >  2:279819/1‑149 (MQ=255)                                                                                                           ggTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGt                              <  5:256618/149‑1 (MQ=255)                                                                                                           ggTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACAATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCg                               <  7:15410/148‑1 (MQ=255)                                                                                                           ggTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACAATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCg                               <  7:20155/148‑1 (MQ=255)                                                                                                             ttGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGc                            <  4:31479/149‑1 (MQ=255)                                                                                                                  cccTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCCACAg                       >  6:136814/1‑149 (MQ=255)                                                                                                                   ccTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCCACAgg                      <  7:193103/149‑1 (MQ=255)                                                                                                                   ccTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCCACAgg                      >  5:120210/1‑149 (MQ=255)                                                                                                                      gCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGCCGCAGGTTGAACTCGGACTGGGTAGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCCACAGGTc                    >  1:206975/1‑148 (MQ=255)                                                                                                                                   aTACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCGACAGGTCGACGGTGGCACCgg      >  7:262486/1‑149 (MQ=255)                                                                                                                                    tACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCCACAGGTCGACGGTGGCACCGGc     >  8:266215/1‑149 (MQ=255)                                                                                                                                       ccGGACGGCACGATCGAGCGCTCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCCACAGGTCGACGGTGGCACCGGCCaa  <  7:266215/149‑1 (MQ=255)                                                                                                                                       ccGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGc                                <  8:193910/119‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               CAAAGTTGATCTGCGTGGTCAGGTGGCGGTGCCCCGGTGCCGAGATGAAGAAGTGCACGTGCGCCGGGCGCTGGCCGTGGCGGCCGAGCAGGTCCAGGCATTCCTGGGTTGGGCCCTGCGGGTCGCAGCCATACCCGGACGGCACGATCGAGCGCGCGC‑‑GG‑TAGCGGCCCTCGG‑‑CATCGGTGATGATACGCCGACGCAGGTTGAACTCGGACTGGGTCGAATCGAAGTACGAATAGGTGCCCTGGGTATTGGCGTGCCACAGGTCGACGGTGGCACCGGCCAA  >  NC_002947/4236052‑4236334

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