Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,271,523 A→C 100% 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,5230AC100.0% 70.4 / NA 24D45E (GAT→GAGuvrXlesion bypass phage DNA polymerase; phage SPbeta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (9/15);  total (9/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCA  >  BS168_reference_genome/2271393‑2271656
                                                                                                                                  |                                                                                                                                     
aaaGCTTCTTTTGTGGATTTACAATGTGTATCCTTTGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTAt                                                                                                                                <  5:222536/138‑1 (MQ=11)
aaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTAtt                                                                                                                               <  4:91023/139‑1 (MQ=12)
aaaGCTTCATTTGGGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCAAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTAtt                                                                                                                               <  2:114941/139‑1 (MQ=11)
       agagactgGATTTACAATGTGTATCCTTGGATCTTAAGGTATCTCAAATAGTCTCGATCGTGTTTTGATTCAAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCGCTac                                                                                                                         <  6:54597/132‑1 (MQ=11)
                    aCAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCa                                                                                                                  >  6:123075/1‑132 (MQ=255)
                    aCAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCa                                                                                                                  <  5:123075/132‑1 (MQ=255)
                     cAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCTATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGATATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGc                                                                                                           <  3:188781/138‑1 (MQ=11)
                          tgtttccTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTATGTGCTGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTa                                                                                                     <  3:96083/133‑1 (MQ=11)
                            tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa                                                                                                    <  2:25911/138‑1 (MQ=255)
                                    gATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATATTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCACTACAACAGCAAGATAGCATGTTAAAGGGTTAAg                                                                                           >  5:168063/1‑139 (MQ=255)
                                                   cAAATAGTCTCGATCCTGTTTTGATTCAAAAATCTTTTTTAAGTGCATGAGATGCAGATAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCACCATTGTTACAGc                                                                            <  6:168063/139‑1 (MQ=14)
                                                          tctcGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATa                                                                       >  8:198045/1‑137 (MQ=255)
                                                              gATCCTGTTTTGATTCCAAAATGTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAg                                                                 <  5:74518/139‑1 (MQ=255)
                                                                                          tAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACaaa                                      <  7:198045/138‑1 (MQ=255)
                                                                                                     gaTGCAGCTAATAGAACACTTCACTGTCTCTCCGTCTTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACTGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTa                           <  5:11524/138‑1 (MQ=14)
                                                                                                        gCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTg                       >  1:14919/1‑139 (MQ=255)
                                                                                                             tAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGAcaca                                        >  4:66851/1‑117 (MQ=255)
                                                                                                             tAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGAcaca                                        <  3:66851/117‑1 (MQ=255)
                                                                                                                     cactTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGACGCATAAAAGGATTGCATGTCGAGACAAAGTATATTCTTACGTGGAAATAGTGAGtt          >  4:127434/1‑138 (MQ=255)
                                                                                                                        ttCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGc                                                                                          <  7:31318/56‑1 (MQ=255)
                                                                                                                        ttCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGc                                                                                          >  8:31318/1‑56 (MQ=255)
                                                                                                                          cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGAcaca                                        <  7:146788/104‑1 (MQ=255)
                                                                                                                          cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGAcaca                                        >  8:146788/1‑104 (MQ=255)
                                                                                                                          cccTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTaatcaa     >  6:131318/1‑139 (MQ=255)
                                                                                                                               tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAatcaatca  <  2:14919/137‑1 (MQ=255)
                                                                                                                               tctctcCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAatcaatca  <  6:52588/137‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
AAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAAAGGGTTAAGCCCCATTGTTACAGCCGATACAGAAGCATAAAAGGATTTCATATCGACACAAAGTATATTCTTACGTGGAAATTGTGAGTAATCAATCA  >  BS168_reference_genome/2271393‑2271656

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: