Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	960,972	T→G	10.8%	K5N (AAA→AAC)	katA ←	vegetative catalase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	960,972	0	T	G	10.8%	93.0 / 3.2	35	K5N (AAA→AAC)	katA	vegetative catalase 1
Reads supporting (aligned to +/- strand):  ref base T (11/20);  new base G (2/2);  total (13/22)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGAACACGTTCTCGGTTGAAATGGGCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATA  >  BS168_reference_genome/960838‑961090                                                                                                                                       |                                                                                                                       tCAGGAACACGTTCTCGGTTGAAATGGGCCAATTTTTCGAGTAAAAGTAAATCATGAATTAAAGTTAGTCCGCGAGAACCGGCAGTCATTGAGTTATGATTAACTCCTACCCGAGCCCCCCAGCAAGTAGACAGTTTAt                                                                                                                    >  3:249555/1‑139 (MQ=255)      aaCACGTTCTCGGTTGAAATGGGCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGaa                                                                                                               <  3:92925/139‑1 (MQ=255)      aaCACGTTCTCGGTTGAAATGGGCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGaa                                                                                                               <  3:265453/139‑1 (MQ=255)        ctcGTTCTTGGTTGAAATTGGCCATTTTTTCGAGTAAATGTAGATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACt                                                                                                             <  8:38914/137‑1 (MQ=255)        cacGTTCTCGGTTGAAATGGGCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACt                                                                                                             <  7:77476/139‑1 (MQ=255)           gTTCTCTGTTGAAATTGGCCAATTTATCGAGTAAATGTACATCTTGAATTACAGGTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCAt                                                                                                          <  4:12921/139‑1 (MQ=255)           gTTCTCGGTTGAAATTGGCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCAt                                                                                                          >  4:208259/1‑139 (MQ=255)                       aaTGGGCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGGTTATTTGAACTc                                                                                                            <  8:279290/125‑1 (MQ=255)                       aaTGGGCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGGTTATTTGAACTc                                                                                                            >  7:279290/1‑125 (MQ=255)                           ggCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTgg                                                                                           >  7:283049/1‑138 (MQ=255)                           ggCCAAATTATCGAGTAACTCTACATCTTGAATTAAAGTAGCTCCGCTAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGGTTATTTGAACTCATGTTATCACCTCTTg                                                                                            <  6:153954/137‑1 (MQ=255)                            gCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGaa                                                                                         <  3:83928/139‑1 (MQ=255)                             ccAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTAtt                                                                                                                   <  4:132562/112‑1 (MQ=255)                             ccAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTAtt                                                                                                                   >  3:132562/1‑112 (MQ=255)                                   tttCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTtata                                                                                   <  1:292750/138‑1 (MQ=255)                                    atCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATaa                                                                                  <  5:288650/137‑1 (MQ=255)                                      cGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGaa                                                                               >  2:157116/1‑139 (MQ=255)                                                     cTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCCGCGAGTTGTCAGGTTATTTGAACTCATGTTATGACCTCTTGGAAGTTATAAGAACATAATATAAtctga                                                                >  2:113911/1‑137 (MQ=255)                                                                   ggTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCaa                                                  >  3:94860/1‑139 (MQ=255)                                                                     tCCGCGAGCACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAata                                                <  5:8577/139‑1 (MQ=255)                                                                        gcgAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAAtatta                                             <  2:190343/139‑1 (MQ=255)                                                                              aCCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACACATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTta                                        >  4:276550/1‑138 (MQ=255)                                                                                      tCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAAttatta                               <  4:94860/139‑1 (MQ=255)                                                                                                ttGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGt                      >  4:105544/1‑138 (MQ=255)                                                                                                ttGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGt                      >  2:35864/1‑138 (MQ=255)                                                                                                   aTTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTAt                   <  3:208259/138‑1 (MQ=255)                                                                                                           cAACCGGATCGCACCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTGTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTTTTTTAAAATAGTTTTCTCTATTaa          <  4:249555/139‑1 (MQ=255)                                                                                                           cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTAtt            <  2:169538/137‑1 (MQ=255)                                                                                                           cAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTAtt            <  2:155255/137‑1 (MQ=255)                                                                                                            aaCCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAAt         <  1:157116/139‑1 (MQ=255)                                                                                                                   tcgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGatata  <  3:276550/138‑1 (MQ=255)                                                                                                                    cgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGatata  >  2:195581/1‑138 (MQ=255)                                                                                                                    cgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGatata  >  2:264280/1‑138 (MQ=255)                                                                                                                    cgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGatata  <  1:264280/138‑1 (MQ=255)                                                                                                                    cgcCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGatata  <  1:195581/138‑1 (MQ=255)                                                                                                                                       |                                                                                                                       TCAGGAACACGTTCTCGGTTGAAATGGGCCAATTTTTCGAGTAAATGTACATCTTGAATTAAAGTTGGTCCGCGAGAACCGGCAGTCATTGAGTTTTGATTATCTCCAACCGGAGCGCCCCAGCTAGTTGTCAGTTTATTTGAACTCATGTTATCACCTCTTGGAATTTATAAGAACATAATATCATCTCTAAGTAAAAAGTCAATATTATTTTATAATTATTATAAAATAGTTATCTCTATTAATCGATATA  >  BS168_reference_genome/960838‑961090

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