Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,368,850 C→A 25.0% E389* (GAA→TAA)  fadA ← acetyl‑CoA C‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,368,8500CA25.0% 71.0 / 2.6 44E389* (GAA→TAA) fadAacetyl‑CoA C‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (15/17);  new base A (2/10);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 8.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATCATTTTATGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCG  >  BS168_reference_genome/3368713‑3368984
                                                                                                                                          |                                                                                                                                        
cAATCATTTTATGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCa                                                                                                                                         <  3:50181/139‑1 (MQ=255)
           tGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑gg                                                                                                                             <  4:150417/138‑1 (MQ=255)
           tGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGc                                                                                                                            <  6:106451/139‑1 (MQ=255)
               cGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCt                                                                                                                        >  7:18920/1‑139 (MQ=255)
                gTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGc                                                                                                                         <  4:320014/137‑1 (MQ=255)
                gTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTc                                                                                                                       >  7:254610/1‑139 (MQ=255)
                       aaaTCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAATTCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTcc                                                                                                                <  4:233056/139‑1 (MQ=255)
                         aTCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTccgc                                                                                                              <  8:18920/139‑1 (MQ=255)
                              ccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTccgccgccg                                                                                                         >  3:119033/1‑139 (MQ=255)
                                    tataCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATacac                                                                                                   <  7:274581/139‑1 (MQ=255)
                                           ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCCCCGATACCAATGGTCa                                                                                            >  2:85288/1‑139 (MQ=255)
                                            tGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTACCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATACACATGGta                                                                                             >  7:290273/1‑136 (MQ=255)
                                              aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCAATTTGTTTTTCCCCC‑TTTTGCAAAATTAAAATA‑CACC‑GGCCGCGCCCCATACCCCCCCGAAAAACATGGTccaga                                                                                         >  4:33588/1‑135 (MQ=255)
                                              aTCGTATGGGACATCCTCGATTAAAAATCCGCCGCATTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATACACATGGTCACg                                                                                          <  4:71388/138‑1 (MQ=255)
                                                            ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGGTCg                                                                           >  4:321112/1‑139 (MQ=255)
                                                            ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCtgg                                                                           >  2:108177/1‑137 (MQ=255)
                                                            ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATACAAATGGTCACGACTACAAAATGCTCg                                                                           >  8:173276/1‑139 (MQ=255)
                                                            ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGGTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCGCCCAGTCCGCCCCCGAAACACCAGGTGCCGGCTGCAAATTGCtgg                                                                           >  1:310844/1‑137 (MQ=255)
                                                              tCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTccgcc                                                                                                             >  8:364317/1‑103 (MQ=255)
                                                              tCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTccgcc                                                                                                             <  7:364317/103‑1 (MQ=255)
                                                              tCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCACCGCCGATACACATGGTCACGACTCCAAATTGCTCGtt                                                                         >  5:151194/1‑139 (MQ=255)
                                                                          ccgccgCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTAAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGGAGGTTTCTCCGtttcct                                                             >  8:288372/1‑137 (MQ=255)
                                                                           cgccgcCTTTTTGTACGTCAGCCGCTTTTTTTCCAATTTGTTTTTCCCCC‑TTTAGAATAATTCAAATA‑CCCC‑GGCCGCCCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGtttcatt                                                            <  1:36394/139‑1 (MQ=255)
                                                                            gccgccTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTccgc                                                                                                              <  5:348036/88‑1 (MQ=255)
                                                                            gccgccTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTccgc                                                                                                              >  6:348036/1‑88 (MQ=255)
                                                                             ccgccTTTTTGTACGTCCGCCGCTTTTTTTCCCATTTTTTTTCCCCCC‑TATCCAATATTTAAAATA‑CACC‑GCCCGCCCCCATCCCCCCGCCCATACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGtttcatttc                                                          <  7:26018/139‑1 (MQ=255)
                                                                                 gtttttGTAGGTCTGCCGCTTTGTTGGCATTTTGTTTTTCCCCA‑TTTAGCATAATTAAAATA‑CCCG‑GGCCGCACCCTTTCCGCCGCCGATACAATGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATgg                                                      <  8:158067/138‑1 (MQ=255)
                                                                                    tttGTCCGTCAGCCGCTTTTTTTGCCATTTGTTTTCCCCCA‑TTTCGCATAATACAAAAC‑CCCC‑GGCCGCCCCCATCCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAt                                                    <  8:243204/138‑1 (MQ=255)
                                                                                    tttGTAGGTCCGCCTTTTTTTTTTCAATTTTTTTTCCCACT‑TTTACAATAATTAAAAAA‑CCCC‑GCCCGCGCCCATCCCGCCGCCGAAACAAAGGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAt                                                    <  8:336707/138‑1 (MQ=255)
                                                                                     ttGTACGTCCGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATACAAAAA‑CCCC‑GCCCGCCCCCATCCCGCCGCCGATACACAGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtg                                                   <  4:124164/138‑1 (MQ=255)
                                                                                        taggtaaGCCGATTTTTTTGCAATTTGTTTTTCCCCT‑TTTAGCATAATAAAAAAA‑CACC‑GGCCGCCCCCATTCCGCCCCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                               <  6:72423/133‑1 (MQ=255)
                                                                                        tACGTCCGCCGCTTTTTTTGCCTTTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CCCC‑GGCCGCCCCCATTCCGCCGCCGATACACATGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                               <  8:192113/139‑1 (MQ=255)
                                                                                         ccgccAGCCGCTTTTTTTGCATTTTGTTTTCCCCCA‑TTTACAATAATACAAAAA‑CCCC‑GCCCGCCCCCATCCCGCCGCCGAAAAACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                               <  7:288372/134‑1 (MQ=255)
                                                                                         acggccgCCGCTTTTTTTCCCATTTTTTTTCCCCCC‑TTTCCCAAAATACAAAAA‑CACC‑GGCCGCCCCCAATCCCCCCCCGATCCACAGGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                              <  6:356062/133‑1 (MQ=255)
                                                                                          cGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATCCACAAGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtgag                                             >  5:356062/1‑139 (MQ=255)
                                                                                          cGTCAGCCGATTTTTTTTCAATTTGTTTTTCCCCC‑TTTAGAATAATTAAAAAC‑CACC‑GCCCGCACCAATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                              <  8:290273/138‑1 (MQ=255)
                                                                                                      tttgtcccATTTGTTTCACCCCA‑TAAAGCATATTAAAAATC‑CTCC‑GCCCGCGCCCAGCCCTCCGCCAATACACTGGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCtt                                  <  6:209237/132‑1 (MQ=255)
                                                                                                                 tGTTTTTCCCCC‑TTAACAAAAAATAAAATA‑CCCC‑GGCCGCTCCAATGCCGCCGCCTAAACACATGGTCACGACTCCAAATTGCTCGGTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                       <  5:311499/138‑1 (MQ=255)
                                                                                                                  ttttttCCCCCATTTAGAATATAAAAAAAA‑CCCG‑GCCCGCGCCCATCCCGCCGCAAATACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                       <  3:234087/137‑1 (MQ=255)
                                                                                                                  tttttaCCCCC‑CTTAACAAAAATAAAAAA‑ACCCGGGCCGCTCCAATCCCGCCGCCGAAACACAGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                       <  1:135808/137‑1 (MQ=255)
                                                                                                                       tCCCCC‑ATTAACAAAATTAAAAAACCCCC‑GGCCGCTCCCATTCCCCCGCCTATAAACATGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCa                  <  3:33588/138‑1 (MQ=255)
                                                                                                                        ccccc‑TTTAGCATAATTCAAAAA‑CCCC‑GGCCGCCCCCATACCCCCGCCGAAACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGc               <  8:39988/139‑1 (MQ=255)
                                                                                                                                   cATAATTCAAATA‑CACC‑GGCCGCTACCAGTCCGCCGCCGATACACAAGGTCACGACTCCAAATTGCTCGTTTCTACGTTACATTTAATGGATGAGTGAGAGTGTAAGCGTTTTGACCGTACAGCCAAGCGCGTGTcc       >  3:79445/1‑137 (MQ=255)
                                                                                                                                      aaTTCAAATA‑CACC‑GGCCGCTCCCATTCCGGCGCCGCTACCCCTGGTCGCGACTCCATATTGGTCGTTTGTGCGTTTCGTTTCATGGTGGAGTGCGTGGGTGAGCTTGGTGGCCGTTCCGCCGAGCGGGGGTccccgcg  >  4:295853/1‑134 (MQ=255)
                                                                                                                                          |                                                                                                                                        
CAATCATTTTATGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCC‑TTTAGCATAATTCAAATA‑CACC‑GGCCGCTCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCG  >  BS168_reference_genome/3368713‑3368984

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