Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,337,477	A→G	13.6%	S25G (AGC→GGC)	PP_3812 →	lipase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,337,477	0	A	G	13.6%	67.2 / 3.9	30	S25G (AGC→GGC)	PP_3812	lipase
Reads supporting (aligned to +/- strand):  ref base A (13/12);  new base G (2/2);  total (15/15)
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.82e-01

AGTCGCAGCGATGGCCCCAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGG  >  NC_002947/4337334‑4337619                                                                                                                                                |                                                                                                                                               aGTCGCAGCGATGGCCCCAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTg                                                                                                                                          <  2:348327/150‑1 (MQ=255)     gCAGCGATGGCCCCAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGgtgt                                                                                                                                      >  1:7296/1‑150 (MQ=255)     gCAGCGATGGCCCCAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGgtgt                                                                                                                                      >  1:56584/1‑150 (MQ=255)       aGCGATGGCCCCAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCg                                                                                                                                    >  2:346504/1‑150 (MQ=255)                cccAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGaa                                                                                                                                  >  1:170576/1‑143 (MQ=255)                cccAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGaa                                                                                                                                  <  2:170576/143‑1 (MQ=255)                 ccAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAacacac                                                                                                                          <  2:84694/150‑1 (MQ=255)                                       aCTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTg                                                                                                    >  2:357879/1‑150 (MQ=255)                                                    gCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTg                                                                                       >  1:371141/1‑150 (MQ=255)                                                      acaGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCt                                                                                                              <  1:583985/125‑1 (MQ=255)                                                      acaGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCt                                                                                                              >  2:583985/1‑125 (MQ=255)                                                        aGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCaa                                                                                   >  1:469806/1‑150 (MQ=255)                                                               gATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGGAGCCCTTGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTg                                                                            <  1:235404/150‑1 (MQ=255)                                                                      ccATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGc                                                                     <  1:357879/150‑1 (MQ=255)                                                                         tGAACATTGTCAACAAAGCCCACACCCCCTCCCTCTTTCCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGACAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTca                                                                  <  1:346504/150‑1 (MQ=255)                                                                                               caccacCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGt                                            <  2:7296/150‑1 (MQ=255)                                                                                                accaccTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCTCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCacac                                                               <  1:576803/130‑1 (MQ=255)                                                                                                accaccTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCTCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCacac                                                               >  2:576803/1‑130 (MQ=255)                                                                                                     ctccctcTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGCCCGCGCGTCAATCTTTCAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGAc                                      <  2:16502/150‑1 (MQ=255)                                                                                                               gcgcTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGACCGGTGCCCAg                            >  1:445031/4‑150 (MQ=255)                                                                                                               gCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAacac                                                                                                                            <  1:453991/54‑1 (MQ=255)                                                                                                               gCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAacac                                                                                                                            >  2:453991/1‑54 (MQ=255)                                                                                                                             cAGTGCCTTCGCCGCAGGCGGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCt                                                                             >  1:44312/1‑87 (MQ=255)                                                                                                                             cAGTGCCTTCGCCGCAGGCGGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCt                                                                             <  2:44312/87‑1 (MQ=255)                                                                                                                             cAGTGCCTTCGCCGCAGGCGGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCt                                                                             <  2:44302/87‑1 (MQ=255)                                                                                                                             cAGTGCCTTCGCCGCAGGCGGCCCTGGTGTCGAACAACACACCCAGGCCTTCCACGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCt                                                                             >  1:44302/1‑87 (MQ=255)                                                                                                                                tGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGa           <  2:56584/150‑1 (MQ=255)                                                                                                                                   cTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCt        <  2:371141/150‑1 (MQ=255)                                                                                                                                    ttCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGGGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGACCGGTGCCCGGGCTTCGGTAAAGGGCGGCCTa       >  1:336177/1‑150 (MQ=255)                                                                                                                                         cgcAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCATCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCggg  >  1:238022/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                               AGTCGCAGCGATGGCCCCAAAGGCCAGGCAAGACGCCCACTACCAACACTTGCACAGAAAAGGATTCGCCCATGAACATTGTCCACAAAGCCCTCACCACCTCCCTCTTGGCCCTGTCCGTTTCCAGTGCCTTCGCCGCAGGCAGCCCTGGTGTCGAACAACACACCCAGGCCTTCCTCGAAGCCCTGGAGAAAGGGGGTGGCAAGCCGCTGGAACAGCTCACACCAAACGACGCCCGCGCGGTACTGACCGGTGCCCAGGCTTCGGTAAAGGTCGACCTATCGGG  >  NC_002947/4337334‑4337619

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑