Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,080,047 C→T 16.6% A209V (GCC→GTC)  ecsC → putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,080,0470CT16.6% 114.0 / 12.8 50A209V (GCC→GTC) ecsCputative hydrolase
Reads supporting (aligned to +/- strand):  ref base C (20/20);  new base T (2/6);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TCATCGAAAATTGGGAGACGGAGAAAAAGAGCATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAGCTGATGATAGGTGATCGCAGCGCCGGCAAG  >  BS168_reference_genome/1079915‑1080185
                                                                                                                                    |                                                                                                                                          
tCATCGAAAATTGGGAGACGGAGAAAAAGAGCATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGTCAATta                                                                                                                                      >  4:380092/1‑139 (MQ=255)
tCATCGAAAATTGGGAGACGGAGAAAAAGAGCATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATta                                                                                                                                      >  3:313213/1‑139 (MQ=255)
 cATCGAAAATTGGGAGACGGAGAAAAAGAGCATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGACGCTTTTTGAGCTTTTGCCGGGCGTAGGAGCGGCAGGCGGCGGAATTGCCAATta                                                                                                                                      >  5:262256/1‑138 (MQ=255)
             ggatccGGATAAAAAGAGCATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCtgcttgct                                                                                                                         <  4:50840/134‑1 (MQ=255)
             ggAGACGGTGAAAATGAGCATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCtgcttgct                                                                                                                         <  6:15630/139‑1 (MQ=255)
                             agCATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACgg                                                                                                         >  4:106567/1‑139 (MQ=255)
                               cATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGct                                                                                                       >  1:409011/1‑138 (MQ=255)
                                   gACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGAc                                                                                                   <  4:91674/139‑1 (MQ=255)
                                    aCTGGAGAGTGTTCCAGCAGTATTTTCGGGATTATTTCGACGCCGTGAAGCTTTTTCAGGTTTTGCCGGGCGTAGGTGCGGCAGTGGGCGGAATTGTCAATAATAAGCTGCATGCTCAGCTTGGATTATCTGCAAGACa                                                                                                  <  8:127126/139‑1 (MQ=255)
                                    aCTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGATTTTTCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGAc                                                                                                   >  3:322325/1‑138 (MQ=255)
                                      tGGTGAGTGTTCCAGCAGGATTATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTGTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATg                                                                                                <  1:104498/139‑1 (MQ=255)
                                             tgtTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCa                                                                                         >  1:91254/1‑139 (MQ=255)
                                               ttCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGGCAATTATAAGCTGCATGGTCAGCTTGGAGGAACGGCAAGCCATGTTTTTCATc                                                                                       >  1:424180/1‑139 (MQ=255)
                                               ttCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCtgcttgct                                                                                                                         >  7:374036/1‑105 (MQ=255)
                                               ttCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCtgcttgct                                                                                                                         <  8:374036/105‑1 (MQ=255)
                                               ttCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCAt                                                                                        <  2:166944/138‑1 (MQ=255)
                                               ttCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCAt                                                                                        <  8:34504/138‑1 (MQ=255)
                                               ttCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATc                                                                                       >  2:378181/1‑139 (MQ=255)
                                                 ccAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGTCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCt                                                                                      <  5:141663/138‑1 (MQ=255)
                                                 ccAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGTCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTg                                                                                     <  4:302271/139‑1 (MQ=255)
                                                  cagcagGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCtgcttgct                                                                                                                         >  4:54102/1‑102 (MQ=255)
                                                  cagcagGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCtgcttgct                                                                                                                         <  3:54102/102‑1 (MQ=255)
                                                   agcagGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGTCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAg                                                                                   <  3:380092/139‑1 (MQ=255)
                                                   agcagGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAg                                                                                   <  4:313213/139‑1 (MQ=255)
                                                     cagGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTGAGCTTTTTCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACAAGTTTTTCATCTGAGGa                                                                                 >  3:86948/1‑139 (MQ=255)
                                                       ggAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGata                                                                               >  2:230264/1‑139 (MQ=255)
                                                       ggAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGGAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGACGata                                                                               >  3:383107/1‑139 (MQ=255)
                                                       ggAATATCGCGATTATATCGACGTAGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGata                                                                               <  4:86948/139‑1 (MQ=255)
                                                                        tagACGTCGTGAAGCGGTTTAAGCTTTTTCCGGGGGTAGGAGCGGCAGTCGGCGGATATGCCAATTTAAATCTGCTTTCTAACAGTGGAGAAACCGCAAGCCAAGTTTTTAATATGAGGATTTTTAATTAAAAAGACgg                                                              >  8:205951/3‑139 (MQ=255)
                                                                           aCGTCGTGAAGCTTTTTCAGCGTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGTGCAACGGCAAGACATGTTTTTCATGTGAGGATATTAAAGGAAACAGCCggcg                                                            >  3:21596/1‑136 (MQ=255)
                                                                             gtcgtGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGaa                                                          >  8:116113/1‑138 (MQ=255)
                                                                                             aGCTTTTGTCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGGTTGCTCAGATTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGt                                          <  4:322325/138‑1 (MQ=255)
                                                                                                ttttGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGc                                                                                                                                            <  1:149138/37‑1 (MQ=255)
                                                                                                ttttGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGc                                                                                                                                            >  2:149138/1‑37 (MQ=255)
                                                                                                      cGGGCGTCGGAGCGGCAGTCGGCGGAATTGTCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAg                                <  1:377204/139‑1 (MQ=255)
                                                                                                      cGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAg                                <  4:311119/139‑1 (MQ=255)
                                                                                                      cGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAg                                <  1:378181/139‑1 (MQ=255)
                                                                                                         gCGTAGGAGCGGCAGTCGGCGGATTTGTCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCTAGGCATGTTTTTCATCGGAGGATTTTTAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTGTGTGTATAGCt                              >  5:236349/1‑138 (MQ=255)
                                                                                                                  cGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAACCACTCCGGCTGTTTTTTTATAGCTGATGATAGGt                    >  5:187779/1‑139 (MQ=255)
                                                                                                                  cGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAGCTGATGATAGGt                    <  5:277480/139‑1 (MQ=255)
                                                                                                                    gCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGata                                                                               <  2:348699/78‑1 (MQ=255)
                                                                                                                    gCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAACCACTCCGGCTGTTTTTTTATAGCTGATGATAGGTGa                  >  7:417520/1‑139 (MQ=255)
                                                                                                                    gCAGTCGGCGGAATTGACAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGata                                                                               >  1:348699/1‑78 (MQ=255)
                                                                                                                       gTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCt                                                                                      <  1:268212/68‑1 (MQ=255)
                                                                                                                       gTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTCGAGAAACGGCAAGACATGTTTTTCATCt                                                                                      >  2:268212/1‑68 (MQ=255)
                                                                                                                         cggcggAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAGCt                              >  6:112621/1‑122 (MQ=255)
                                                                                                                         cggcggAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAGCt                              <  5:112621/122‑1 (MQ=255)
                                                                                                                            cggAATTGTCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGTAGAATAAACACTCCGGCTGTTTTTTTATAGCTGATGATAGGTGATCGCAgcg          <  4:327990/139‑1 (MQ=255)
                                                                                                                               aaTTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACa                                                                                                  >  2:375531/1‑48 (MQ=255)
                                                                                                                               aaTTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACa                                                                                                  <  1:375531/48‑1 (MQ=255)
                                                                                                                                   gCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAGCTGATGATAGGTGATCGCAGCGCCGGCa    <  2:91254/138‑1 (MQ=255)
                                                                                                                                   gCCAATTATAACCTGCTGGATCAGCTTGGAGAACCGGAAAGAAATGTTTTTCATCTGAGTATATTAAAGCAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAGCTGATGATAGGTGATCGCAGCGCCGGCa    <  4:21596/138‑1 (MQ=255)
                                                                                                                                    ccAATTATCAGCTGCTTGCTCTGGTTGGAGAAACGGAAAGAAATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAGCTGATGATAGGTGATCGCAGCGCCGGCAAg  <  4:360823/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                          
TCATCGAAAATTGGGAGACGGAGAAAAAGAGCATCGACTGGAGAGTGTTCCAGCAGGAATATCGCGATTATATCGACGTCGTGAAGCTTTTTCAGCTTTTGCCGGGCGTAGGAGCGGCAGTCGGCGGAATTGCCAATTATAAGCTGCTTGCTCAGCTTGGAGAAACGGCAAGACATGTTTTTCATCTGAGGATATTAAAGGAAACAGCCGGAGAATAAACACTCCGGCTGTTTTTTTATAGCTGATGATAGGTGATCGCAGCGCCGGCAAG  >  BS168_reference_genome/1079915‑1080185

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