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_genome2,271,5230AC98.1% 179.1 / NA 53D45E (GAT→GAG) uvrXlesion bypass phage DNA polymerase; phage SPbeta
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (28/24);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGAAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCATT  >  BS168_reference_genome/2271388‑2271658
                                                                                                                                       |                                                                                                                                       
gatgaAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTcc                                                                                                                                      >  7:60205/1‑139 (MQ=12)
    aaaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                                                                                                                                   <  4:106443/138‑1 (MQ=12)
           tCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCAAAAATCTTTTTTAATTGCAGGAGATGAAGCTAATACTACACTTCCCTGACACTCCGTATACCCTac                                                                                                                           >  8:205455/1‑139 (MQ=11)
             aTTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTacaa                                                                                                                         <  3:167126/139‑1 (MQ=12)
                   ggATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCa                                                                                                                    <  5:294496/138‑1 (MQ=255)
                    gATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGc                                                                                                                     >  7:329818/1‑136 (MQ=255)
                        tACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATa                                                                                                               <  4:245799/138‑1 (MQ=255)
                         aCAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAg                                                                                                              >  4:10887/1‑138 (MQ=255)
                          cAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGc                                                                                                             <  8:329818/138‑1 (MQ=255)
                               tgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGtt                                                                                                        >  1:81812/1‑138 (MQ=255)
                                 tATCCTTGTATCTTCAGGTATCTCAAATATTCTCGATCCTTTTTTGATTCCAAAAACTTTTTTAAGTGCAGGATATGAATAAAATACTACACTTCCATGTCTCTCATTATTACCTACAACAGAAATATAGCATTTTa                                                                                                       >  8:169788/1‑137 (MQ=2)
                                 tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaaa                                                                                                     <  8:185197/139‑1 (MQ=255)
                                 tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa                                                                                                      <  3:269857/138‑1 (MQ=255)
                                 tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTa                                                                                                       >  6:144527/1‑137 (MQ=255)
                                       tGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTa                                                                                               <  8:103964/139‑1 (MQ=255)
                                        ggATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa                                                                                                      >  2:252588/1‑131 (MQ=21)
                                        ggATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa                                                                                                      <  1:252588/131‑1 (MQ=21)
                                        ggATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTa                                                                                               <  8:21253/138‑1 (MQ=255)
                                         gATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAg                                                                                             >  6:281918/1‑139 (MQ=255)
                                               cAGTTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATCCTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCa                                                                                        <  1:230833/138‑1 (MQ=255)
                                               cAGGTATCTCATATTGTCTCGATCCTGTTTTGTTTCCAAAATCTTTTTTAAGTGCAGGCGATGCTGCTAATACTAGACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCa                                                                                        <  7:205455/138‑1 (MQ=11)
                                               cAGGTATCTCAAATAGTTTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCAt                                                                                       >  6:100677/1‑139 (MQ=255)
                                               cAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGcccc                                                                                         <  4:115813/137‑1 (MQ=255)
                                               cAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCa                                                                                        <  5:144527/138‑1 (MQ=255)
                                                     tctcAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGtt                                                                                   <  4:160275/137‑1 (MQ=255)
                                                     tctcAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGtt                                                                                   >  3:160275/1‑137 (MQ=255)
                                                     tctcAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTa                                                                                  >  1:31833/1‑138 (MQ=255)
                                                         aaaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACa                                                                                <  8:60205/136‑1 (MQ=255)
                                                         aaaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGc                                                                              <  5:281918/138‑1 (MQ=255)
                                                          aaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGc                                                                              <  6:290980/137‑1 (MQ=255)
                                                          aaTAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGc                                                                              >  5:290980/1‑137 (MQ=255)
                                                               tctcGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa                                                                                                      <  2:154291/108‑1 (MQ=255)
                                                               tctcGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa                                                                                                      >  1:154291/1‑108 (MQ=255)
                                                                         gTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGGTAAGCACAATTGTTACAGCCGATACAGAAGCATaa                                                              >  3:195011/1‑138 (MQ=14)
                                                                          ttttGATTCCAACATCTTTTTTAAGTGCTTGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATaa                                                              <  5:277113/137‑1 (MQ=14)
                                                                                             tttATGTGCTGGAGATGAAGCTAATAATACACTTCCCTTTCTATCCGTATTCACTACAAAAGAAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACATCAGATACAGAAGCATAAAAGGATTTCATATCGAcaca                                          <  7:169788/138‑1 (MQ=14)
                                                                                                           aTGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTAct                           >  6:9105/1‑138 (MQ=255)
                                                                                                             gCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTg                         >  8:63209/1‑139 (MQ=255)
                                                                                                             gCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTg                         >  3:332114/1‑139 (MQ=255)
                                                                                                                       ctacactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGT‑GAAATTGTGAg              <  4:335081/139‑1 (MQ=255)
                                                                                                                          gacaacccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTa            <  5:66806/134‑1 (MQ=255)
                                                                                                                          cactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATATCATGTTAAAGTGTTAAGACACATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGt             >  2:238039/1‑138 (MQ=14)
                                                                                                                          cactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGt             >  3:303083/1‑138 (MQ=255)
                                                                                                                          cactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTa            <  2:81812/139‑1 (MQ=255)
                                                                                                                          cactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCACATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGt             >  8:189704/1‑138 (MQ=255)
                                                                                                                           actTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaa           >  8:77580/1‑139 (MQ=255)
                                                                                                                           actTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaa           >  6:117927/1‑139 (MQ=255)
                                                                                                                            ctTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaat          <  5:117927/139‑1 (MQ=255)
                                                                                                                               gccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaatca        >  2:73288/2‑138 (MQ=255)
                                                                                                                               cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaatcaa       >  4:182516/1‑139 (MQ=255)
                                                                                                                               cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaatcaa       >  4:373124/1‑139 (MQ=255)
                                                                                                                               cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaatcaa       >  6:111867/1‑139 (MQ=255)
                                                                                                                               cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaatca        >  7:133585/1‑138 (MQ=255)
                                                                                                                                    tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAatcaatca    <  3:10887/137‑1 (MQ=255)
                                                                                                                                    tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCAtt  <  8:133585/139‑1 (MQ=255)
                                                                                                                                    tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCAtt  >  3:186268/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                       
GATGAAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCATT  >  BS168_reference_genome/2271388‑2271658

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