Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,338,414	T→A	7.5%	intergenic (+34/+108)	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,414	0	T	A	7.5%	152.7 / 3.0	53	intergenic (+34/+108)	PP_3812/PP_5632	lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base T (27/22);  new base A (2/2);  total (29/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 = 5.47e-01

TACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCG  >  NC_002947/4338270‑4338562                                                                                                                                                 |                                                                                                                                                      tACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGc                                                                                                                                                   >  2:214433/1‑149 (MQ=255)  aCAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTg                                                                                                                                                    <  8:349485/147‑1 (MQ=255)  aCAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTg                                                                                                                                                    >  7:349485/1‑147 (MQ=255)     aCGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTc                                                                                                                                                >  5:210516/1‑148 (MQ=255)     aCGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATACACCCGCCCAGGCGGGTCTATTGCtcat                                                                                                                                               <  1:11158/149‑4 (MQ=255)        gAATGATTCATGACTATGGGCTGCTC‑ACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAttt                                                                                                                                           >  2:344111/1‑149 (MQ=255)         aaTGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAttt                                                                                                                                           >  1:199904/1‑149 (MQ=255)                cATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTGGCTTCTATTTGCGGaaa                                                                                                                                    >  6:307932/1‑149 (MQ=255)                     cTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGcc                                                                                                                               <  2:97402/149‑1 (MQ=255)                      tATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGccc                                                                                                                              >  4:29101/1‑149 (MQ=255)                      tATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGccc                                                                                                                              <  7:335298/149‑1 (MQ=255)                               ctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCtt                                                                                                                      <  6:194960/148‑1 (MQ=255)                               ctCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCtt                                                                                                                      >  1:269440/1‑148 (MQ=255)                                tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAGCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTc                                                                                                                     >  5:348096/1‑148 (MQ=255)                                tCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTc                                                                                                                     >  6:69974/1‑148 (MQ=255)                                  aaCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAAt                                                                                                                  >  5:153362/1‑149 (MQ=255)                                  aaCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAAt                                                                                                                  <  1:214433/149‑1 (MQ=255)                                  aaCCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTAAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAAt                                                                                                                  <  7:141588/149‑1 (MQ=255)                                     ccGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGaa                                                                                                               <  8:123380/149‑1 (MQ=255)                                               aaGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCg                                                                                                     >  4:295261/1‑149 (MQ=255)                                                aGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCgg                                                                                                    >  1:192805/1‑149 (MQ=255)                                                  tGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTg                                                                                                  <  2:192805/149‑1 (MQ=255)                                                    ccGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTc                                                                                                                                                <  5:178770/101‑1 (MQ=255)                                                    ccGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTc                                                                                                                                                >  6:178770/1‑101 (MQ=255)                                                    ccGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGc                                                                                                 <  6:238884/148‑1 (MQ=255)                                                      ggCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑Gc                                                                                             <  6:210516/149‑1 (MQ=255)                                                         aGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGc                                                                                                 >  7:309526/1‑143 (MQ=255)                                                                   gcgcTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCtt                                                                                >  2:306585/1‑149 (MQ=255)                                                                       tGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                                         <  7:339035/89‑1 (MQ=255)                                                                       tGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGc                                                                                                                                         >  8:339035/1‑89 (MQ=255)                                                                              ggCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCATTAGCGGGTaaa                                                                     >  3:100129/1‑149 (MQ=255)                                                                                cGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAAcc                                                                   <  2:199904/149‑1 (MQ=255)                                                                                        gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAg                                                           <  6:153362/149‑1 (MQ=255)                                                                                        gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                           <  3:295261/149‑1 (MQ=255)                                                                                          aCTTAAAGTGCATCTGCAGTGAGGCCATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGTGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAgg                                                          <  1:67460/149‑1 (MQ=255)                                                                                           cTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCATTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGc                                                        <  1:306585/149‑1 (MQ=255)                                                                                             tAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCATTAGCCGGGAAAACCCCGAAAAGGCCg                                                      >  5:1833/1‑149 (MQ=255)                                                                                              aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCgg                                                     >  8:259726/1‑149 (MQ=255)                                                                                               aaGTGCATCTGCAGTGAGGCGATCAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATCAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCgg                                                     <  3:246799/148‑1 (MQ=255)                                                                                                            gTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                       <  5:69974/149‑1 (MQ=255)                                                                                                            gTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                       <  2:269440/149‑1 (MQ=255)                                                                                                                  atagaaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                          >  4:46121/6‑124 (MQ=38)                                                                                                                  atagaaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                          <  3:46121/119‑1 (MQ=38)                                                                                                                        gCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCt                           >  2:95748/1‑149 (MQ=255)                                                                                                                         cAAAAGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCAATAGCGGGTAAACCCCCGAAAGGGCCCGGAAATACTGCACAGGGCTCAAACCCta                          >  4:349535/1‑148 (MQ=255)                                                                                                                         cAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt                          >  3:121189/1‑149 (MQ=255)                                                                                                                          aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                         >  5:368304/1‑149 (MQ=255)                                                                                                                          aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                         <  1:95748/149‑1 (MQ=255)                                                                                                                          aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                         >  2:50111/1‑149 (MQ=255)                                                                                                                           aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGTTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                        <  1:344111/149‑1 (MQ=255)                                                                                                                               aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTAGCGGGTTAACCCGCCAAACGGCCCGCAAAACCGGCTCAGGGCTCAAACCCATCcaca                     >  2:327237/1‑148 (MQ=255)                                                                                                                                 ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                   >  1:243463/1‑148 (MQ=255)                                                                                                                                 ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCATTCGCGGGGGAAACCGCGAAACGGCCGGCACATACAGGTCAGGGCGGAAACCCGGCcaca                     >  1:157397/1‑146 (MQ=255)                                                                                                                                                ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATACTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTc   >  3:302985/1‑149 (MQ=255)                                                                                                                                                 tttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTAGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCg  >  5:84817/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      TACAACGGAATGATTCATGACTATGGGCTGCTCAACCCGTTGAGCCAAGTGCCGGCAGTAAAGGCGGCGCTGCGGCAGGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCG‑GCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCG  >  NC_002947/4338270‑4338562

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