Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,338,493	A→T	13.8%	intergenic (+113/+29)	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,493	0	A	T	13.8%	93.8 / 8.7	51	intergenic (+113/+29)	PP_3812/PP_5632	lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (24/20);  new base T (3/4);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

GCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCG  >  NC_002947/4338348‑4338639                                                                                                                                                  |                                                                                                                                                   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                                                                                                     >  3:248493/1‑110 (MQ=255)                                                                                    aaaCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                     <  4:248493/110‑1 (MQ=255)                                                                                      aCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCt                                                            <  1:266282/149‑1 (MQ=255)                                                                                      aCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCt                                                            <  1:259548/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                          >  4:161386/1‑134 (MQ=255)                                                                                                                                       cTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAg                          <  3:161386/134‑1 (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  <  3:5430/148‑1 (MQ=255)                                                                                                                                                 tAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCg  >  4:5430/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCG  >  NC_002947/4338348‑4338639

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