Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,992,352 C→T 100% L755L (TTG→TTA dnaEC ← DNA polymerase III (alpha subunit), DnaE3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,992,3520CT100.0% 130.4 / NA 38L755L (TTG→TTAdnaECDNA polymerase III (alpha subunit), DnaE3
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/16);  total (22/16)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACAGATCCGTTTTCAACTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGCAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTAAAGCCATAGTTTGCGAATTTGACGATTAAATCGTAGAC  >  BS168_reference_genome/2992221‑2992480
                                                                                                                                   |                                                                                                                                
aCAGATCCGTTTTCAACTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCg                                                                                                                            <  1:40015/138‑1 (MQ=255)
aCAGATCCGTTTTCAACTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCg                                                                                                                            >  6:3215/1‑138 (MQ=255)
aCAGATCCGTTTTCAACTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCg                                                                                                                            >  5:121721/1‑138 (MQ=255)
aCAGATCCGTTTTAAACTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCCCTCCCCTTTGCTTCGTAAACTATTTGAGAGATTTTATCCTCATTGACAATAACGCATGGAAGTAATCCGc                                                                                                                           >  2:89979/1‑139 (MQ=255)
      ccGTTATCAACTTTCAACGAAAAAGTGTTTTTTATCAATGCCGGAGGAAGTATACGAATCCCGCTCCCCTTTCCTTCGTAAATATATTGAGAGTTTTTTTCCTCTTTGCCAATCACGATTTTCATTAACCCGCACTTa                                                                                                                      <  2:89281/138‑1 (MQ=255)
            tCAACTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAatat                                                                                                                <  8:209428/138‑1 (MQ=255)
              aaCTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAAc                                                                                                             <  6:123353/139‑1 (MQ=255)
              aaCTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAAc                                                                                                             >  4:103780/1‑139 (MQ=255)
              aaCTGTAAAAGGAAAACTGCTTTTATTCACTGACTGATGTAGGATACGAATCCCGCTCCCCTTTGCTTCTTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATaa                                                                                                              <  3:23243/138‑1 (MQ=255)
                   tAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGAt                                                                                                         >  4:103358/1‑138 (MQ=255)
                     aaCGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATaa                                                                                                       <  2:45492/138‑1 (MQ=255)
                          aaaaCTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCt                                                                                                 >  4:172943/1‑139 (MQ=255)
                                 cTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCaa                                                                                            <  1:89979/137‑1 (MQ=255)
                                    ttATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACg                                                                                       >  3:210377/1‑139 (MQ=255)
                                    ttATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACg                                                                                       >  3:104252/1‑139 (MQ=255)
                                            tGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTgg                                                                               >  5:74407/1‑139 (MQ=255)
                                                            aaaaTCCCTCTCCCATTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTCTCAGTATTCCGCACATAAATAATAACTCATAATGAGATTTCAAATACGCAAGCTGGAACCCGATCCTGCTGt                                                               >  8:110335/3‑139 (MQ=255)
                                                                tCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGc                                                           <  3:103358/139‑1 (MQ=255)
                                                                  ccGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCa                                                         >  7:4774/1‑139 (MQ=255)
                                                                     cTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCg                                                      >  5:217894/1‑139 (MQ=255)
                                                                             tGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCAcc                                                                           >  6:80865/1‑110 (MQ=255)
                                                                             tGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCAcc                                                                           <  5:80865/110‑1 (MQ=255)
                                                                              gCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCtt                                              <  3:172943/138‑1 (MQ=255)
                                                                                  cGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGc                                                           <  3:73761/121‑1 (MQ=255)
                                                                                  cGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGc                                                           >  4:73761/1‑121 (MQ=255)
                                                                                  cGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTcctgt                                         >  7:117059/1‑135 (MQ=255)
                                                                                  cGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGt                                          >  1:38875/1‑138 (MQ=255)
                                                                                  cGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGt                                          >  6:188632/1‑138 (MQ=255)
                                                                                  cGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGt                                          >  5:56358/1‑138 (MQ=255)
                                                                                   gTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTa                                        <  6:74407/139‑1 (MQ=255)
                                                                                    tAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATa                                                                                                                      >  7:139292/1‑60 (MQ=255)
                                                                                    tAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATa                                                                                                                      <  8:139292/60‑1 (MQ=255)
                                                                                       aGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTAAAGc                                    <  4:210377/139‑1 (MQ=255)
                                                                                                   aTTTTATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTAAAGCCATAGTTTGCg                         >  4:801/1‑138 (MQ=255)
                                                                                                     tttATCCTCATTGCCAATGACGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTAAAGCCATAGTTTGCGAAt                      >  5:22187/1‑139 (MQ=255)
                                                                                                             cATTGCCAATGACGCTTGTCAATAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTAAAGCCATAGTTTGCGAATTTGACGAt              >  1:68455/1‑139 (MQ=255)
                                                                                                                         cGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTAAAGCCATAGTTTGCGAATTTGACGATTAAATCGTAGAc  <  6:22187/139‑1 (MQ=255)
                                                                                                                         cGCTTGTCAGTAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTAAAGCCATAGTTTGCGAATTTGACGATTAAATCGTAGAc  <  3:57637/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                
ACAGATCCGTTTTCAACTGTAAACGGAAAACTGCTTTTATTCACTGACGGAGGAAGGATACGAATCCCGCTCCCCTTTGCTTCGTAAAGATATTGAGAGATTTTATCCTCATTGCCAATGACGCTTGTCAGCAATCCGCACATAAAATATAACGGATAATGAGCTTTCAAATACGCAAGCTGGCACCCGATCATGCTGTATGCCACCGCATGGCTTCTGTTAAAGCCATAGTTTGCGAATTTGACGATTAAATCGTAGAC  >  BS168_reference_genome/2992221‑2992480

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: