Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,782,538 G→A 7.9% F125F (TTC→TTT)  atpD ← ATP synthase (subunit beta, component F1)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,782,5380GA 7.9% 153.7 / 3.6 52F125F (TTC→TTT) atpDATP synthase (subunit beta, component F1)
Reads supporting (aligned to +/- strand):  ref base G (27/20);  new base A (2/2);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GGATTAATACGGTTTTACCTACACCGGCACCACCGAACAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAACCGGTACTGAGATTGGCGCT  >  BS168_reference_genome/3782408‑3782673
                                                                                                                                  |                                                                                                                                       
ggATTAATACGGTTTTACCTACACCGGCACCACCGAACAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAgg                                                                                                                                 >  6:70865/1‑139 (MQ=255)
            ttttACCTACACCGGCACCACCGAACAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCctgtctgt                                                                                                                      >  1:211241/1‑138 (MQ=255)
                 ccTACACCGGCACCACCGAACAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCg                                                                                                                >  7:96846/1‑139 (MQ=255)
                 ccTACACCGGCACCACCGAACAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCg                                                                                                                >  2:139391/1‑139 (MQ=255)
                     cacCGGCACCACCGAACAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATc                                                                                                            >  7:194387/1‑139 (MQ=255)
                          gCACCACCGAACAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACATGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCtttt                                                                                                       >  8:186071/1‑139 (MQ=255)
                                acGAACAATCCGATTTTAACGCCCTTAATGTAAGCAGCAACCAAAGAAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCAt                                                                                                 <  6:214510/138‑1 (MQ=255)
                                    aCAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCg                                                                                              <  2:270065/138‑1 (MQ=255)
                                        tCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAgg                                                                                          >  6:85795/1‑138 (MQ=255)
                                        tCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAgg                                                                                          >  5:40768/1‑138 (MQ=255)
                                        tCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCAAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAgg                                                                                          >  5:172416/1‑138 (MQ=255)
                                        tCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCAAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAgg                                                                                          >  5:168408/1‑138 (MQ=255)
                                             ttttACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGACCCgg                                                                                    >  6:202172/1‑139 (MQ=255)
                                               ttACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGc                                                                                   <  1:59118/138‑1 (MQ=255)
                                                   cGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCAAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCAtt                                                                              <  6:168408/139‑1 (MQ=255)
                                                   cGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCAAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCAtt                                                                              <  6:172416/139‑1 (MQ=255)
                                                     cccTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCa                                                                            <  2:86994/139‑1 (MQ=255)
                                                                aGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAg                                                                                                                                  <  7:295270/74‑1 (MQ=255)
                                                                aGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAg                                                                                                                                  >  8:295270/1‑74 (MQ=255)
                                                                aGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATAttt                                                                 >  1:100170/1‑139 (MQ=255)
                                                                aGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATAttt                                                                 >  1:213236/1‑139 (MQ=255)
                                                                aGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTACGCATCCGCAGGAACCGGCTCATACAAATCAATAttt                                                                 >  2:159143/1‑139 (MQ=255)
                                                                     caagcaaATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTcc                                                            >  4:165953/1‑139 (MQ=255)
                                                                         caaATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCgaga                                                        <  1:299932/139‑1 (MQ=255)
                                                                            aTCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACg                                                     >  5:174515/1‑139 (MQ=255)
                                                                                       tAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACg                                          >  6:75105/1‑139 (MQ=255)
                                                                                         aTACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGt                                         <  2:211241/138‑1 (MQ=255)
                                                                                          tACGTGTTTCAATAATTTCAACTTCTGTTGAAAGCTGATCGAATGACGGCGCCTGGCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTcc                                       <  2:100170/139‑1 (MQ=255)
                                                                                          tACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTcc                                       <  8:194387/139‑1 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAAc                                                                                       <  2:49282/86‑1 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAAc                                                                                       >  1:49282/1‑86 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTcc                                       <  3:261869/134‑1 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTcc                                       >  4:261869/1‑134 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAgtg                                  >  8:121543/1‑139 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAgtg                                  >  7:77060/1‑139 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAgtg                                  >  8:112779/1‑139 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAgtg                                  >  7:28196/1‑139 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAgtg                                  >  7:206067/1‑139 (MQ=255)
                                                                                               gTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAgtg                                  >  8:123529/1‑139 (MQ=255)
                                                                                                  acaacaaTTACAACTTCTGTGTAAACCTGAACCACTGACGGCGCCTGGCTGTGAATCGGACCCTTTTTCCAAGCCGCAGGAACCGGCTCATTCAAATCAATATTTTCACCGAGAACGTTAAATACACGTCCAAGTGtt                                <  6:40768/137‑1 (MQ=255)
                                                                                                         tttCCACTTCTGTTTAAAGCTGATCGCCTGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGACAACGTTAAATACACGTCCAAGTGTTACATCAcc                        <  6:104745/139‑1 (MQ=255)
                                                                                                         tttCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCAcc                        <  7:200939/139‑1 (MQ=255)
                                                                                                             aaCTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAAc                     <  7:112779/138‑1 (MQ=255)
                                                                                                             aaCTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAAc                     <  3:145311/138‑1 (MQ=255)
                                                                                                             aaCTTCTGTTGAAAGCTGAGCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAAc                     <  8:96846/138‑1 (MQ=255)
                                                                                                             aaCTTCTGTTGAAAGCGGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAAc                     <  6:107683/138‑1 (MQ=255)
                                                                                                              aCTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAACCg                   >  3:312720/1‑139 (MQ=255)
                                                                                                                             tGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAACCGGTACTGAGATTGgcg    >  5:84637/1‑139 (MQ=255)
                                                                                                                               aTCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAACCGGTACTGAGATTGgcg    <  7:123529/137‑1 (MQ=255)
                                                                                                                               aTCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAACCGGTACTGAGATTGgcg    <  8:206067/137‑1 (MQ=255)
                                                                                                                               aTCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAACCGGTACTGAGATTGGCGCt  >  4:211822/1‑139 (MQ=255)
                                                                                                                               aTCGAATGAAGGCGCCTGGCTGTGAATCGTATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAACCGGTACTGAGATTGgcgc   <  7:90696/138‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                       
GGATTAATACGGTTTTACCTACACCGGCACCACCGAACAATCCGATTTTACCGCCCTTAATGTAAGGAGCAAGCAAATCAACAACTTTAATACCTGTTTCAAGAATTTCAACTTCTGTTGAAAGCTGATCGAATGAAGGCGCCTGTCTGTGAATCGGATCCTTTTTCGCATCCGCAGGAACCGGCTCATTCAAATCAATATTTTCTCCGAGAACGTTAAATACACGTCCAAGTGTTACATCACCAACCGGTACTGAGATTGGCGCT  >  BS168_reference_genome/3782408‑3782673

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