Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,521,291	C→T	8.5%	W165* (TGG→TGA)	lip ←	lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,521,291	0	C	T	8.5%	134.6 / 2.0	47	W165* (TGG→TGA)	lip	lipase
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base T (2/2);  total (23/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 = 9.94e-01

TGACCTCATAGGCCCCTTCCCCACCCCAGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACGTGACCGATGCCACCCGCCCCG  >  NC_002947/5521143‑5521439                                                                                                                                                      |                                                                                                                                                     tGACCTCATAGGCCCCTTCCCCACCCCAGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGc                                                                                                                                                      >  8:313644/1‑149 (MQ=255)  gACCTCATAGGCCCCTTCCCCACCCCAGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGc                                                                                                                                                      <  3:48610/148‑1 (MQ=255)      tCATAGGCCCCTTCCCCACCCCAGGTTTGCGGCATGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCacac                                                                                                                                                 >  1:236694/1‑149 (MQ=255)          aGGCCCCTTCCCCACCCCAGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACcag                                                                                                                                              >  6:253715/1‑148 (MQ=255)          aGGCCCCTTCCCCACCCCAGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACcag                                                                                                                                              >  2:209684/1‑148 (MQ=255)              cccTTCCCCACCCCAGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAgcagca                                                                                                                                         >  3:302806/1‑149 (MQ=255)               ccTTCCCCACCCCAGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCAc                                                                                                                                        <  6:356535/149‑1 (MQ=255)                    cccaccccaGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCa                                                                                                                                    >  3:220792/1‑148 (MQ=255)                                   cGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGcaca                                                                                                                     <  1:144324/148‑1 (MQ=255)                                         gCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGa                                                                                                              <  4:220792/149‑1 (MQ=255)                                           ccTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACgg                                                                                                                                      >  8:121147/1‑123 (MQ=255)                                           ccTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAg                                                                                                             >  1:80781/1‑148 (MQ=255)                                               tgGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGAt                                                                                                          <  2:236694/147‑1 (MQ=255)                                                      gCCTGGTTGAACAGCGCCACCCCGGCGCTTGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCATGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCt                                                                                                  <  7:76445/148‑1 (MQ=255)                                                      gCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGa                                                                                                              <  6:234141/136‑1 (MQ=255)                                                      gCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGa                                                                                                              >  5:234141/1‑136 (MQ=255)                                                       ccTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTgg                                                                                                <  1:209684/149‑1 (MQ=255)                                                       ccTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTgg                                                                                                <  4:302806/149‑1 (MQ=255)                                                        cTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTggg                                                                                               <  2:80781/149‑1 (MQ=255)                                                        cTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTggg                                                                                               <  3:367113/149‑1 (MQ=255)                                                             tGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGaa                                                                                          >  8:14691/1‑149 (MQ=255)                                                              gAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACATCCTTGAGGATGCGCCCCTGGGGTGaa                                                                                          <  4:263907/148‑1 (MQ=255)                                                                       cACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTg                                                                                >  5:82190/1‑149 (MQ=255)                                                                             ggCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGc                                                                          <  5:253715/149‑1 (MQ=255)                                                                                gcTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGtt                                                                       >  7:296633/1‑149 (MQ=255)                                                                                    gTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCa                                                                    <  7:253139/148‑1 (MQ=255)                                                                                         cGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAggtggt                                                              >  8:202810/1‑149 (MQ=255)                                                                                            cTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCgg                                                           <  7:202810/149‑1 (MQ=255)                                                                                               gTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGTCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGc                                                          >  2:94118/1‑147 (MQ=255)                                                                                                  cGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCTCAGCTTTGTGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCatc                                                      <  5:259970/148‑3 (MQ=255)                                                                                                   gAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGTCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTc                                                    <  1:94118/149‑1 (MQ=255)                                                                                                          gCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCa                                             >  2:130263/1‑149 (MQ=255)                                                                                                           cataTCTACCGGCAGCGGGTCACGACGC‑CAGGCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCAt                                            <  3:339977/146‑1 (MQ=255)                                                                                                            acaTCTACCGGCAGCGGGTCACGACGCGCAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCAt                                            >  1:204292/1‑149 (MQ=255)                                                                                                              aTCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGTCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGt                                         >  5:29018/1‑149 (MQ=255)                                                                                                                   ccGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGtt                                        <  7:313644/145‑1 (MQ=255)                                                                                                                   ccGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGc                                    >  5:104738/1‑149 (MQ=255)                                                                                                                   ccGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGc                                    <  1:130263/149‑1 (MQ=255)                                                                                                                      gCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAgg                                                                  >  6:57030/1‑116 (MQ=255)                                                                                                                      gCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAgg                                                                  <  5:57030/116‑1 (MQ=255)                                                                                                                      gCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGcccc                                 <  6:82190/149‑1 (MQ=255)                                                                                                                        aGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGc                               >  5:221575/1‑149 (MQ=255)                                                                                                                             gTCACGACGC‑CAGCCGGTTTCCATTCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGTGCCCCGCCACCg                          <  6:29018/149‑1 (MQ=255)                                                                                                                                acgacgC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACg                       <  6:104738/149‑1 (MQ=255)                                                                                                                                acgacgC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACg                       >  1:201510/1‑149 (MQ=255)                                                                                                                                acgacgC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACg                       >  2:289066/1‑149 (MQ=255)                                                                                                                                  gacgC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACGTg                     >  8:186259/1‑149 (MQ=255)                                                                                                                                    cgC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACGTGAc                   <  1:289066/149‑1 (MQ=255)                                                                                                                                                    aGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACGTGACCGATGCCACCCGccc    >  1:343698/1‑149 (MQ=255)                                                                                                                                                      ccACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACGTGACCGATGCCACCCGCCCCg  <  7:14691/149‑1 (MQ=255)                                                                                                                                                      ccACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACGTGACCGATGCCACCCGCCCCg  <  3:326880/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TGACCTCATAGGCCCCTTCCCCACCCCAGGTTTGCGGCAGGCCCTGTGGATAAGCCTGGTTGAACAGCGCCACCCCGGCGCTGGTCAGCGACTGGTGCGAGGCGTGCACATCTACCGGCAGCGGGTCACGACGC‑CAGCCGGTTTCCAGCCACACCAGCAGCACGGCAAACCCATGCAGCACAGCCTTGAGGATGCGCCCCTGGGGTGAATTGCCCGGTGCCGTGCGTTCCAGGTGGTCGGCCAGCTCAGAACCATGGTTGGGCCCCGCCACCGACGTGACCGATGCCACCCGCCCCG  >  NC_002947/5521143‑5521439

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