Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,758,999	A→G	9.7%	S360P (TCT→CCT)	sacC ←	levanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,758,999	0	A	G	9.7%	112.7 / 4.2	41	S360P (TCT→CCT)	sacC	levanase
Reads supporting (aligned to +/- strand):  ref base A (23/14);  new base G (2/2);  total (25/16)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

AAAACCAAATTCAGCAGCTGAACCAGGGCTTACTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTTCTCCATGGGGATGTCGGAACATCATTCG  >  BS168_reference_genome/2758871‑2759130                                                                                                                                 |                                                                                                                                    aaaaCCAAATTCAGCAGCTGAACCAGGGCTTACTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTAAGAGGTTCcgcg                                                                                                                            >  4:207054/1‑138 (MQ=255) aaaaCCAAATTCAGCAGCTGAACCAGGGCTTACTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGAGATGGTCAGATTCTTCCACTTCTTATAGGTTCcgcg                                                                                                                            >  1:114626/1‑138 (MQ=255)               cagcTGAACCAGGGCTTACTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGctc                                                                                                             <  2:69692/139‑1 (MQ=255)               cagcTGAACCAGGGCTTACTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGGCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGctc                                                                                                             <  7:295063/139‑1 (MQ=255)                      aCCAGGGCTTACTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCAcc                                                                                                      <  3:207054/139‑1 (MQ=255)                               tACTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTgg                                                                                             >  2:38745/1‑139 (MQ=255)                                aCTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCCCCGGGGTTTggc                                                                                            >  6:215984/1‑138 (MQ=255)                                 cTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGa                                                                                            <  4:33197/138‑1 (MQ=255)                                 cTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGa                                                                                            <  2:267090/138‑1 (MQ=255)                                 cTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGa                                                                                            <  2:93280/138‑1 (MQ=255)                                 cTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATAGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCATATTCTTCCACTTCTTAGAGGTTCCGCGATTGGTTTCCAGCTCTTTCACCGGTGTTTGGAc                                                                                           <  3:29971/139‑1 (MQ=255)                                     aaaTTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCAct                                                                                       >  7:120535/1‑139 (MQ=255)                                      aaTTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCActc                                                                                      >  1:165105/1‑139 (MQ=255)                                            gCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAAccc                                                                                >  5:82815/1‑139 (MQ=255)                                            gCATATATTTCATAGGCATCTCCAGAATGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCAGCGAATGGTTTCCAGCTCTTTAACCCGTGGTTGGACCACTCAAAccc                                                                                >  5:199920/1‑139 (MQ=255)                                                       aTAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGGGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAAc                                                                     >  8:47641/1‑139 (MQ=255)                                                       aTAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGGGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAAc                                                                     >  6:89819/1‑139 (MQ=255)                                                            cttctcCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGGCTGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCttt                                                                <  1:55477/137‑1 (MQ=255)                                                              tctcCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTACACGCttta                                                               >  2:277342/1‑137 (MQ=255)                                                                     aTTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGTTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTCACCCCTTCGGTAAACGCTTTCAAtttt                                                         >  1:133232/1‑137 (MQ=255)                                                                     aTTGCCCAGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACTCTTTCAATTTTaa                                                       <  8:247929/139‑1 (MQ=255)                                                                       tGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAAct                                                     >  6:295795/1‑139 (MQ=255)                                                                         ccccGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCCCTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAActct                                                   >  2:78244/1‑139 (MQ=255)                                                                         ccccGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGAc                                                                                           <  4:37889/99‑1 (MQ=255)                                                                         ccccGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGAc                                                                                           >  3:37889/1‑99 (MQ=255)                                                                         ccccGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAActct                                                   >  4:302590/1‑139 (MQ=255)                                                                         ccccGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAActct                                                   >  3:6136/1‑139 (MQ=255)                                                                          cccGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAActctc                                                  <  3:133161/139‑1 (MQ=255)                                                                                      ttATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACg                                      <  4:260784/139‑1 (MQ=255)                                                                                       tATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCCGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGt                                     >  4:29853/1‑139 (MQ=255)                                                                                            cTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCAc                                <  3:302590/139‑1 (MQ=255)                                                                                               gCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGCCCCCTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTtc                             >  6:127958/1‑139 (MQ=255)                                                                                               gCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGCCCACTCTAACCCCTTCGGTACACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTtc                             >  8:178844/1‑139 (MQ=255)                                                                                               gCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACtt                              <  1:38745/138‑1 (MQ=255)                                                                                               gCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTtc                             >  6:280778/1‑139 (MQ=255)                                                                                                                     tcCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTTCTCCATGGGGATGTCGGAAcat        >  1:62119/1‑138 (MQ=255)                                                                                                                     tcCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTTCTCCATGGGGATGTCGGAAcat        >  8:30119/1‑138 (MQ=255)                                                                                                                      ccACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTTCTCCATGGGGATGTCGGAAcatca      >  5:212036/1‑139 (MQ=255)                                                                                                                        aCTTCTTAGGGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTTCTCCATGGGGATGTCGGAACATCAtt    <  5:89819/139‑1 (MQ=255)                                                                                                                        aCTTCTTAGGGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTTCTCCATGGGGATGTCGGAACATCAtt    <  7:47641/139‑1 (MQ=255)                                                                                                                          ttcttAGAGGTTCCGCGAATGGTTTCCATCGCTTTCACAGGTGTTTGGACCACTGTAACCCCTTCGGTAAAAGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCCATTCTCCATGGGGATGTGGGATCATCAttag  >  4:14473/1‑137 (MQ=255)                                                                                                                                 |                                                                                                                                    AAAACCAAATTCAGCAGCTGAACCAGGGCTTACTTTAAATTCTGCATTTATTTCATAGGCATCTCCAGATTGCCCCGCTAACACATTATGACTTGCAGGCGATATGGTCAGATTCTTCCACTTCTTAGAGGTTCCGCGAATGGTTTCCAGCTCTTTCACCGGTGTTTGGACCACTCTAACCCCTTCGGTAAACGCTTTCAATTTTAACTCTCTTGGAATGGACGTTGCACTTCTCCATGGGGATGTCGGAACATCATTCG  >  BS168_reference_genome/2758871‑2759130

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