Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,271,424 T→C 100% R78R (AGA→AGG 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,4240TC100.0% 98.7 / ‑6.6 33R78R (AGA→AGGuvrXlesion bypass phage DNA polymerase; phage SPbeta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (12/20);  minor base A (1/0);  total (13/20)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAA  >  BS168_reference_genome/2271292‑2271558
                                                                                                                                    |                                                                                                                                      
tCCTGCATCTATAAAAGATTCATCAATTGAATACGTATGGACACTTTTCTCAGGAAGAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGTATTTACAATGTGTATCCTTgg                                                                                                                                   <  1:50041/138‑1 (MQ=14)
tCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGa                                                                                                                                  >  6:200380/1‑139 (MQ=255)
                        aaTTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGtct                                                                                                           <  4:294672/138‑1 (MQ=255)
                            gAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGAt                                                                                                       <  7:123814/138‑1 (MQ=255)
                                       aaaacaTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGATTCATTTGTGGATTTACAATTTCTATACTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGAtt                                                                                           >  4:114749/4‑139 (MQ=11)
                                            aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  8:99955/1‑138 (MQ=37)
                                            aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  1:63679/1‑138 (MQ=37)
                                            aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  6:209233/1‑138 (MQ=37)
                                             tttCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaaaa                                                                                     <  8:8978/139‑1 (MQ=37)
                                              ttCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAg                                                                                                                           <  4:297110/100‑1 (MQ=255)
                                              ttCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAg                                                                                                                           >  3:297110/1‑100 (MQ=255)
                                                   aGGAACAAATCTGTAAAACAGCTATGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTATGATTCCAACATCtttt                                                                               >  2:154133/1‑139 (MQ=255)
                                                          aaTCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTg                                                                        >  8:52338/1‑139 (MQ=21)
                                                          aaTCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTg                                                                        >  2:169607/1‑139 (MQ=21)
                                                           aTCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTg                                                                        >  2:113296/1‑138 (MQ=21)
                                                                  aaaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGAt                                                                <  5:34520/139‑1 (MQ=21)
                                                                  aaaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCGTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGAGCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGAt                                                                <  6:139473/139‑1 (MQ=255)
                                                                  aaaCAGCGTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTCTCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGAt                                                                <  2:248442/139‑1 (MQ=255)
                                                                   aaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATg                                                               <  5:237011/139‑1 (MQ=21)
                                                                       gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  3:105938/139‑1 (MQ=21)
                                                                        cTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTg                                                                                                                                    >  8:231002/1‑65 (MQ=255)
                                                                        cTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTg                                                                                                                                    <  7:231002/65‑1 (MQ=255)
                                                                        cTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           >  6:65095/1‑138 (MQ=21)
                                                                                                    aaaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                               <  3:200763/138‑1 (MQ=12)
                                                                                                             aTTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTacaa                     <  5:209233/139‑1 (MQ=12)
                                                                                                                 gtgGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAg                  <  5:65095/138‑1 (MQ=255)
                                                                                                                    gATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                               <  2:267178/122‑1 (MQ=12)
                                                                                                                    gATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                               >  1:267178/1‑122 (MQ=12)
                                                                                                                          cAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGc         <  1:288890/138‑1 (MQ=255)
                                                                                                                          cAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGc         <  5:200380/138‑1 (MQ=255)
                                                                                                                          cAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGc         <  7:52338/138‑1 (MQ=255)
                                                                                                                          cAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGc         <  1:169607/138‑1 (MQ=255)
                                                                                                                                 tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa  <  1:113296/138‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                      
TCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAA  >  BS168_reference_genome/2271292‑2271558

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: