Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,987,550 T→C 7.1% intergenic (‑3/+181) pfkA ← / ← accA 6‑phosphofructokinase/acetyl‑CoA carboxylase (carboxyltransferase alpha subunit)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,987,5500TC 7.1% 170.2 / 3.8 56intergenic (‑3/+181)pfkA/accA6‑phosphofructokinase/acetyl‑CoA carboxylase (carboxyltransferase alpha subunit)
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base C (2/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TCAATCCCGCGTATCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCTTACCTATTATACAATTAAATAAATTGTTTTT  >  BS168_reference_genome/2987415‑2987672
                                                                                                                                       |                                                                                                                          
tCAATCCCGCGTATCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCa                                                                                                                         <  5:91880/139‑1 (MQ=255)
   aTCCCGCGTATCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCAtt                                                                                                                       >  3:134523/1‑138 (MQ=255)
      ccGCGTATCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGTACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCAc                                                                                                                    <  7:124678/138‑1 (MQ=255)
         cgTATCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCa                                                                                                                >  1:108792/1‑139 (MQ=255)
          gTATCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAg                                                                                                               <  8:99658/139‑1 (MQ=255)
             tCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCa                                                                                                             <  5:115732/138‑1 (MQ=255)
             tCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCa                                                                                                             <  3:67969/138‑1 (MQ=255)
                  tGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACata                                                                                                        >  2:236424/1‑138 (MQ=255)
                    tAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTGATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACatata                                                                                                      <  4:206623/138‑1 (MQ=255)
                      aaTACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGa                                                                                                   <  5:203046/139‑1 (MQ=255)
                           cGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTGCATTCACCTAAGCAACATAAATTATTaa                                                                                               >  4:35062/1‑138 (MQ=255)
                           cGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCAtt                                                                                                                       <  8:72292/114‑1 (MQ=255)
                           cGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCAtt                                                                                                                       >  7:72292/1‑114 (MQ=255)
                            gTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAaac                                                                                             <  3:286930/139‑1 (MQ=255)
                              aaaCTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAaaca                                                                                            >  4:291256/1‑138 (MQ=255)
                                  ttCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATAAATGATTAaacataaca                                                                                       >  7:78478/1‑139 (MQ=255)
                                    ctaCGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCTTTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCTCCTGAGCAACATATATGATTAaacataaca                                                                                       <  4:152537/135‑1 (MQ=255)
                                     aaCGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGa                                                                                    >  7:68594/1‑139 (MQ=255)
                                            tGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCCCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGc                                                                             >  1:233947/1‑139 (MQ=255)
                                            tGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCCCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGc                                                                             >  1:20350/1‑139 (MQ=255)
                                                  aTCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGa                                                                                                   <  3:46058/111‑1 (MQ=255)
                                                  aTCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGa                                                                                                   >  4:46058/1‑111 (MQ=255)
                                                  aTCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTaa                                                                        >  7:33035/1‑138 (MQ=255)
                                                  aTCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTaa                                                                        >  1:79909/1‑138 (MQ=255)
                                                      cTTTTCTGACTACTTCGCGAAATGCTGCGTTCATTCCCGGGCATTCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCTTTCACCTCAGCAACATATATGATTAAAAATAACATGATCATTGCCCTAAAATCa                                                                   <  8:78478/139‑1 (MQ=255)
                                                                tactGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTaaa                                                                                              >  7:146189/1‑102 (MQ=255)
                                                                tactGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTaaa                                                                                              <  8:146189/102‑1 (MQ=255)
                                                                  ctGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGAt                                                       <  3:67734/139‑1 (MQ=255)
                                                                  ctGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGAt                                                       <  1:236424/139‑1 (MQ=255)
                                                                   tGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACata                                                                                                        >  4:119761/1‑89 (MQ=255)
                                                                   tGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACata                                                                                                        <  3:119761/89‑1 (MQ=255)
                                                                   tGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTaaa                                                                                              >  1:102278/1‑99 (MQ=255)
                                                                   tGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTaaa                                                                                              <  2:102278/99‑1 (MQ=255)
                                                                    gcgcGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTATTACCCCGATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGAtct                                                     >  7:135582/1‑139 (MQ=255)
                                                                     cgcgAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCAtt                                                                                                                       >  8:227211/1‑72 (MQ=255)
                                                                     cgcgAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCAtt                                                                                                                       <  7:227211/72‑1 (MQ=255)
                                                                     cgcgAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGaaaa                                                            >  3:196784/1‑131 (MQ=255)
                                                                     cgcgAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGaaaa                                                            <  4:196784/131‑1 (MQ=255)
                                                                     cgcgAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCCCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGAtc                                                      <  2:233947/137‑1 (MQ=255)
                                                                     cgcgAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCCCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGAtc                                                      <  2:20350/137‑1 (MQ=255)
                                                                      gagAACTGCAGAGTTCATTCCCGGTGAATCCCCGCCGCTCGTTTATACCCCTATTCTTTTCATTCTACAATCACCTCAGCAACTTATATGATTAAACATAACATGATCATTGCCCTAAAACCAATGTAAAGGGATCTcc                                                   <  6:249308/137‑1 (MQ=255)
                                                                         aaCTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCAtc                                                <  2:108792/139‑1 (MQ=255)
                                                                           ctgctgCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCAtcct                                              <  3:291256/139‑1 (MQ=255)
                                                                                      aTTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAAc                                   >  1:17325/1‑139 (MQ=255)
                                                                                      aTTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAAc                                   >  1:158017/1‑139 (MQ=255)
                                                                                      aTTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGGTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAAc                                   >  1:197846/1‑139 (MQ=255)
                                                                                         cccGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCt                                 <  8:68594/138‑1 (MQ=255)
                                                                                           cGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCATCATATAAGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCACTATGAAACCTTAc                              <  8:139378/139‑1 (MQ=255)
                                                                                            ggggAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCTTAcc                             <  2:197846/139‑1 (MQ=255)
                                                                                                 aTCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCTTACCtatta                        <  3:137176/139‑1 (MQ=255)
                                                                                                  tCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCTTACCtatta                        >  8:185062/1‑138 (MQ=255)
                                                                                                                        cTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGAtc                                                      >  7:201295/1‑86 (MQ=255)
                                                                                                                        cTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGAtc                                                      <  8:201295/86‑1 (MQ=255)
                                                                                                                        cTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCTTACCTATTATACAATTAAATAAATTGttttt  >  6:101724/1‑138 (MQ=255)
                                                                                                                                ttCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCTTACCTATTATACAAt                  <  3:245393/114‑1 (MQ=255)
                                                                                                                                ttCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCTTACCTATTATACAAt                  >  4:245393/1‑114 (MQ=255)
                                                                                                                                       |                                                                                                                          
TCAATCCCGCGTATCCGTTGTAAATACCGTAAACTTCAACGTCATGATAGATCGCTTTTCTGACTACTGCGCGAACTGCTGCGTTCATTCCCGGGGAATCCCCGCCGCTCGTTAATACCCCTATTCGTTTCATTCTCCATTCACCTCAGCAACATATATGATTAAACATAACATGATCATTGCCCTAAAAACAATGAAAAGGGATCTCCATCCTCCCTATGAAACCTTACCTATTATACAATTAAATAAATTGTTTTT  >  BS168_reference_genome/2987415‑2987672

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