Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,338,494	A→T	15.2%	intergenic (+114/+28)	PP_3812 → / ← PP_5632	lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,338,494	0	A	T	15.2%	104.4 / 8.0	51	intergenic (+114/+28)	PP_3812/PP_5632	lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (23/20);  new base T (4/4);  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

GCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAAC  >  NC_002947/4338348‑4338642                                                                                                                                                   |                                                                                                                                                     gCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAc                                                                                                                                                    >  4:299660/1‑149 (MQ=255)          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAg                                                                                                                                           >  8:257252/1‑149 (MQ=255)            aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAggg                                                                                                                                         <  3:277189/149‑1 (MQ=255)             cTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGc                                                                                                                                        <  3:299660/149‑1 (MQ=255)              ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGcc                                                                                                                                       <  7:257252/149‑1 (MQ=255)                  aGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAGCCGGGAAAACCGCGAAAGGGCCGGc                                                                                                                                    >  2:200833/1‑148 (MQ=255)                         cTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCAGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCCGGTATACCCGCCAAATGGCCGGCAAATCCt                                                                                                                             >  5:10836/1‑148 (MQ=255)                            cAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAAGCTGCCTGCGGTGCGGCCGGTCCGGCCCTGTGGCGGGGAAACACGGGCAAGGGCCGGCAACTCCTGcg                                                                                                                          >  8:227139/1‑147 (MQ=255)                               tGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcg                                                                                                                                                 <  1:358953/122‑1 (MQ=255)                               tGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcg                                                                                                                                                 >  2:358953/1‑122 (MQ=255)                                   gCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCCAATCCTGGTCAGGGcg                                                                                                                   >  3:208666/1‑147 (MQ=255)                                          gCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAGCCGGTAAACCCGCGAAAGGGCCGGCAAATACAGCTCAGGGCTCAAAaccg                                                                                                           >  8:300826/1‑148 (MQ=255)                                            aatcgccccGCCTGGGCGGTTATTTTGCTTCTATTTGCGGAAACACCCCGCAGGCTTCAATGAATCTGCCAGGGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                         <  7:300826/143‑1 (MQ=255)                                                 aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAGCGGGGAAAACCCCGAAAAGGCCGGCAAAAACAGCACAAGGCACAAAcccatacaaa                                                                                                     >  7:50684/1‑141 (MQ=255)                                                   ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAACGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                   >  1:197885/1‑148 (MQ=255)                                                     gCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGAt                                                                                                >  2:208736/1‑149 (MQ=255)                                                             gggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCt                                                                                           <  4:87647/146‑1 (MQ=255)                                                             gggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCCCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCt                                                                                           >  3:87647/1‑146 (MQ=255)                                                             gggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAAACCGCGAAAGGGGCGGCAAATCCTGCTCAGGGCTCAACACCTTACACAATGATCACCTCAg                                                                                        >  2:296062/1‑149 (MQ=255)                                                              ggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATACTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGc                                                                                       >  1:101644/1‑149 (MQ=255)                                                                     tGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCTAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGc                                                                                 <  5:313772/148‑1 (MQ=255)                                                                     tGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGc                                                                                 <  1:208736/148‑1 (MQ=255)                                                                          cTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCGACAATGATCACCCCAGCCTTCGCCAccc                                                                            <  1:200833/148‑1 (MQ=255)                                                                           tATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAGCGGGTAAACCCGCGCAAGGGCCGGCAAATCCTGCTCAGGGCGCAAACCCTTCCACAATGATCACCACAGCCATCGCCAccccg                                                                          >  2:266282/1‑148 (MQ=255)                                                                                    aaaCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                        <  4:248493/110‑1 (MQ=255)                                                                                    aaaCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                        >  3:248493/1‑110 (MQ=255)                                                                                      aCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCt                                                               <  1:259548/149‑1 (MQ=255)                                                                                      aCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCt                                                               <  1:266282/149‑1 (MQ=255)                                                                                        aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTg                                                              <  7:314805/148‑1 (MQ=255)                                                                                        aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGTAAAACCGCGAAAGGGCCCGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAAGATAAACACAGCCATAGCCACCCCGTGGCGATgggag                                                              >  8:314805/1‑145 (MQ=255)                                                                                                      aaTGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCACAGCCTTCGCCACCCCTTGGCGATGGCGGCCTGCCGCGTGGTa                                                >  6:333344/1‑148 (MQ=255)                                                                                                                     ggTGCGGCAGGTCGGGCCTTTTCGCGGTTAAACCCGCTAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTaa                                <  3:351533/149‑1 (MQ=255)                                                                                                                        gCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAg                             <  1:296062/149‑1 (MQ=255)                                                                                                                           gCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                   <  3:140154/76‑1 (MQ=255)                                                                                                                           gCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCCCCCCTTGGCGATGGCTGCATGCCTCTTGGTTCAGCGCGGAGCGGTAACAGGCg                          >  5:174467/1‑149 (MQ=255)                                                                                                                              ggTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACg                       <  5:333344/149‑1 (MQ=255)                                                                                                                                tCCGGCCCTTCCGCGGGCAAACCCGCGAAAGGGCCGGCAAATCCTGCGCAGGGCTCACACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTGGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGGAACAGGCGACGGc                     <  7:227139/149‑1 (MQ=255)                                                                                                                                  cgtcccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCt                   <  6:10836/146‑1 (MQ=255)                                                                                                                                  cGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCt                   >  6:214503/1‑149 (MQ=255)                                                                                                                                    gCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCGGGGCTCAAACCCTTACACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCGCTTGGTACAGCGCGGAGCGGTAACAGGCGCcggccggc                 >  1:254641/1‑149 (MQ=255)                                                                                                                                       cTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg              >  2:151349/1‑149 (MQ=255)                                                                                                                                       cTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAg                             <  3:161386/134‑1 (MQ=255)                                                                                                                                       cTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAg                             >  4:161386/1‑134 (MQ=255)                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGt                                  <  5:316734/128‑1 (MQ=255)                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGt                                  >  6:316734/1‑128 (MQ=255)                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg              >  4:333823/1‑148 (MQ=255)                                                                                                                                          tCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa            <  1:151349/148‑1 (MQ=255)                                                                                                                                           cgcgGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc           <  5:214503/148‑1 (MQ=255)                                                                                                                                                gTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc      >  3:341019/1‑148 (MQ=255)                                                                                                                                                 tAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCg     >  4:5430/1‑148 (MQ=255)                                                                                                                                                 tAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCg     <  3:5430/148‑1 (MQ=255)                                                                                                                                                   aaCCCTCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAAc  <  2:197885/149‑1 (MQ=255)                                                                                                                                                   aaCCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGaa   <  3:333823/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAAC  >  NC_002947/4338348‑4338642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑