Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	BS168_reference_genome	2,271,523	A→C	D45E (GAT→GAG)	uvrX ←	lesion bypass phage DNA polymerase; phage SPbeta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,271,523	0	A	C	98.0%	164.9 / NA	49	D45E (GAT→GAG)	uvrX	lesion bypass phage DNA polymerase; phage SPbeta
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (21/27);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCATT  >  BS168_reference_genome/2271392‑2271658                                                                                                                                    |                                                                                                                                        aaaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                                                                                                                                   <  5:170001/138‑1 (MQ=12) aaaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                                                                                                                                   <  2:80030/138‑1 (MQ=12) aaaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                                                                                                                                   <  7:322936/138‑1 (MQ=12)      cTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTccc                                                                                                                              >  5:244999/1‑138 (MQ=12)              gtgGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCC‑AAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGc                                                                                                                     <  5:239874/138‑1 (MQ=11)                      aCAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAg                                                                                                              <  8:276833/138‑1 (MQ=255)                       cAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGc                                                                                                             <  7:100979/138‑1 (MQ=255)                            tgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGtt                                                                                                        >  7:215982/1‑138 (MQ=255)                            tgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGtt                                                                                                        >  3:320578/1‑138 (MQ=255)                              taaccTTTGATCTTCAGGTATCTCAAAAACTCTCGATCATGTTTTGATTCCAAATTCTTTTTTATGTGCAGGAGATGCAGCTAATAAAACACATCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaaa                                                                                                     <  3:3972/135‑1 (MQ=11)                              tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa                                                                                                      <  7:195950/138‑1 (MQ=255)                              tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTa                                                                                                       >  3:362643/1‑137 (MQ=255)                                  attGGATCTTACGGTATTTCAAATAGTCTCTATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCTGGAGATGCAGCTAATACTACACTTCCCTGTCTCTACGTATTCCCTACAACAGCAAGATAGCATGTTAAAggg                                                                                                  <  6:385119/137‑1 (MQ=11)                                    tGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGtt                                                                                                >  5:222254/1‑138 (MQ=255)                                     ggATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCATGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGtt                                                                                                >  2:366320/1‑137 (MQ=255)                                      gATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAg                                                                                             >  4:135933/1‑139 (MQ=255)                                           tCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCa                                                                                        >  2:24440/1‑139 (MQ=255)                                           tCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCa                                                                                        >  5:335884/1‑139 (MQ=255)                                            cAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGcccc                                                                                         <  6:335884/137‑1 (MQ=255)                                            cAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCa                                                                                        <  2:29284/138‑1 (MQ=255)                                            cAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCAt                                                                                       >  4:339976/1‑139 (MQ=255)                                                  tctcAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAg                                                                                                              >  6:78380/1‑110 (MQ=16)                                                  tctcAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAg                                                                                                              <  5:78380/110‑1 (MQ=16)                                                    tcAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGtt                                                                                                <  4:181790/122‑1 (MQ=255)                                                    tcAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGtt                                                                                                >  3:181790/1‑122 (MQ=255)                                                    tcAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAg                                                                               >  6:6610/1‑139 (MQ=255)                                                      aaaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCCGGAGATGCATCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGc                                                                              <  1:145655/138‑1 (MQ=255)                                                      aaaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACa                                                                                <  6:222254/136‑1 (MQ=255)                                                      aaaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACa                                                                                <  8:215982/136‑1 (MQ=255)                                                      aaaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGc                                                                              <  4:121807/138‑1 (MQ=255)                                                      aaaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGc                                                                              <  1:68915/138‑1 (MQ=255)                                                       aaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGcc                                                                             >  6:86128/1‑138 (MQ=255)                                                                gATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAg                                                                      <  5:315420/136‑1 (MQ=255)                                                                gATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAg                                                                   <  5:86128/139‑1 (MQ=255)                                                                             ttCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGAttt                                                      <  3:206493/139‑1 (MQ=255)                                                                                  aaaTCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCatat                                                 <  1:336552/139‑1 (MQ=255)                                                                                  aaaTCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCatat                                                 <  1:65362/139‑1 (MQ=255)                                                                                                          gCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTg                         >  4:4927/1‑139 (MQ=255)                                                                                                                    ctacactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATtgtg                <  6:244999/138‑1 (MQ=255)                                                                                                                      acactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAg              >  6:375409/1‑138 (MQ=255)                                                                                                                       cactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTa            <  6:172147/139‑1 (MQ=255)                                                                                                                       cactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTa            <  5:146073/139‑1 (MQ=255)                                                                                                                       cactTCACTGTCTAACCGTAAAACCAACAACAGCAAGATATCAAGATAAACCGTAAAGCCCAATTGTAACAGACGAAACAAAAGAAGAAAAGGATAACATATCGACACAAAGTATATTCTTAAGTCGAAATtgaaata             >  8:346894/1‑133 (MQ=255)                                                                                                                        actTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaa           >  4:52860/1‑139 (MQ=255)                                                                                                                        actTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaa           >  4:373088/1‑139 (MQ=255)                                                                                                                          ttCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGAcaca                                          >  8:263642/1‑106 (MQ=255)                                                                                                                          ttCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGAcaca                                          <  7:263642/106‑1 (MQ=255)                                                                                                                            cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaatcaa       >  5:285428/1‑139 (MQ=255)                                                                                                                            cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaatcaa       >  3:385029/1‑139 (MQ=255)                                                                                                                                 tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAatcaatca    <  4:362643/137‑1 (MQ=255)                                                                                                                                 tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCAtt  <  3:135933/139‑1 (MQ=255)                                                                                                                                 tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCAtt  <  3:4927/139‑1 (MQ=255)                                                                                                                                 tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCAtt  <  1:366320/139‑1 (MQ=255)                                                                                                                                 tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCAtt  >  4:77189/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                        AAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCATT  >  BS168_reference_genome/2271392‑2271658

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