Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,839,474 C→A 29.2% K133N (AAG→AAT yrbE ← putative inositol‑related oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,839,4740CA29.2% 54.1 / 3.2 38K133N (AAG→AATyrbEputative inositol‑related oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (8/15);  new base A (7/2);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 4.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGATCGACATATCCATAAACAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCAGAGAGAAACTAACGGGCTTCTCACA  >  BS168_reference_genome/2839340‑2839590
                                                                                                                                       |                                                                                                                    
tGGATCGACATATCCATAAACAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTa                                                                                                                   >  5:236190/1‑138 (MQ=255)
   aTCGACCTATCCATAAACAATCCGCCTGAAGTACGGACATAATCTATATTCGGTGGTTCTGGGTCCCTGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAAt                                                                                                               <  5:129407/139‑1 (MQ=255)
       aCATATCCATAAACAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGGTTTTTTATAAATTAGTGGTGTCCCGTTCACCATTTTGAAAAAACGTTTTTATTTTTTTTATATGAg                                                                                                            >  1:149056/1‑138 (MQ=255)
       aCATATCCATAAACAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGGTTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTGTTTAAATTAgg                                                                                                           >  5:104768/1‑139 (MQ=255)
       aCATATCCATAAACAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGGGTTTTTATAAATTAGGTGTTCCGATTTGACCAATTTGAACAATAGGTTTTTTTTTTTTTTAAAtatg                                                                                                            >  1:277636/1‑136 (MQ=255)
                tAAACAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGCGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAgcg                                                                                                  <  4:102912/139‑1 (MQ=255)
                tAAACAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAgcg                                                                                                  <  5:83455/139‑1 (MQ=255)
                  aaCAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATTAGGTGTACCCATTTCACCAATTTAAACAATCGTTTTGTTTTTCTTAAAATGAGGAGCAAAccgcg                                                                                                >  7:58070/1‑137 (MQ=255)
                     aaTCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGt                                                                                             <  6:236190/139‑1 (MQ=255)
                       tCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGTTTTTTAATAAATGAGGTGTACCGGTTTCACCATTTTTAACAATCGTTTTGTTTTTTTTTAAATGTTGATAAAAGcggcgggt                                                                                            >  5:255999/1‑136 (MQ=255)
                       tCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGTTTTTTAATAAAAGAGCTGTACCGAGTTGACCAATTAGAACAATCCTGTTGAGTTTGATAAAAAGAGGATCAAAGCCGCCGtt                                                                                            >  8:172574/1‑138 (MQ=255)
                       tCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTTATAAATGTGGTGTACCGATTTCACCATTTTCAACAAACGTTTTGATTTTATTAAAATGAGGATCAAAGcggcgggt                                                                                            >  3:290722/1‑136 (MQ=255)
                       tCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTAACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATAAAAGCGGCGGtt                                                                                            >  8:199387/1‑138 (MQ=255)
                       tCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGAATTCACCAGTTTAAACAATATTTTTGATTTTATTTAAATTAGGATCAAAGagggggat                                                                                            >  4:217680/1‑136 (MQ=255)
                        ccGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTGAACAAACGTTTTGTTTTTATTTAAATTAGGATCAAAGCGGCGGGTa                                                                                           >  2:217997/1‑138 (MQ=255)
                          gccgcaCGTCCGGACAAAATCTATTTTCGGGGGTGCGGGGTCCGGGGACGGTATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTAAAAATGAGGATCAACGCGGCGGTTAAAt                                                                                        <  2:277636/134‑1 (MQ=255)
                                  tACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTg                                                                                <  1:217997/139‑1 (MQ=255)
                                   agagACAGAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTg                                                                                <  3:35923/135‑1 (MQ=255)
                                                      gggggTTCTGGGTCCCGGGACGTGGATTTTAATAAATTAGGTGTTCCGAATTCACCAGAATCAACAAAAGTTTTTGTTTTAGTAAAAAGTAGATAAAACGGGCGCTTAAATACGACTTGAAGAGTTTCCCCCTGTTTcc                                                            >  6:167022/1‑139 (MQ=255)
                                                         ggTTCTGGGTCCCGGGACGTGTTTTTTAATAAATGAGGTGTTAAGATTTAACCATTTTAAACAAACCTTTTTATTTTATTTAAATTTGGATTAAAGCGGCGGTTATATACGCCATGAATTGGTACACCGTTTTTCCggg                                                         >  6:247111/1‑138 (MQ=255)
                                                         ggTTCTGGGGCCCGGGACGTGATTTTTAATAAAGGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGa                                                         <  6:255999/139‑1 (MQ=255)
                                                           ttCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGAATTCACCATTTTCAACAATCGTTTTGATTTTATTAAAATGATGATCAAAGCGGCGGGTAAATACGACATGAAGTGTTAACCCGTGTTTACGGaaa                                                       >  5:136735/1‑137 (MQ=255)
                                                           ttCTGGGTCCCGGGACGTGATTTTTAAAAAATTAGGTGTACCGAATTAACCAATTTAATCAATCGTATTTATTTTATTAAAATTTTGATCAAAGCCGCGGTTTAATCCGAAATGAAGTGTTACCCCGGGTTTCCGGAcc                                                       >  2:189780/1‑138 (MQ=255)
                                                                 gTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCa                                                 <  3:94439/139‑1 (MQ=255)
                                                                  tCCCGGGACGTTATTTTTAATAAATGAGGTGTACCTATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAATTTCGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCaa                                                <  5:167022/139‑1 (MQ=255)
                                                                    ccGGTACGTTATTTTTAATAATTGAGGTGTACCGATTTCACCATTTTAAAAAATCGTTTTTATTTTATTAAAATGATGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAg                                               <  2:154648/138‑1 (MQ=255)
                                                                       ggACGTGTTTTTTAATAAATGAGGTGTTCCGGTTTCACCATTTTGAACAATCGTTTTGATTTTGTTAAAATGAGGATCAAAGCGGCGGGTTAATACGACTTGAAGTGTTACCCCGTGTTTGCGGACCGACGCCAgcgcg                                           >  4:94439/1‑138 (MQ=255)
                                                                           gTGATTTTTAATAACTGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCAc                                       <  3:217680/139‑1 (MQ=255)
                                                                            tGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTAAACAATCGTTTTTATTTTCTTAAAATGAGGATCAAATCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCAGGACCGCCGCAACCCCTTCACt                                      >  6:295258/1‑139 (MQ=255)
                                                                                          tGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCagag                        <  4:174969/139‑1 (MQ=255)
                                                                                             ggttcaCCGATTTCACCATTTTAACAAATCTTTTTTATTTTATTAAAATGAGTATCAACGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCAgagaga                     <  5:35825/134‑1 (MQ=255)
                                                                                                        ttCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCAGAGAGAAACTAACGGGc          >  3:157538/1‑139 (MQ=255)
                                                                                                           aCCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCAGAGAGAAACTAACGGGCTtc       >  5:109419/1‑139 (MQ=255)
                                                                                                               ttttCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCAGAGAGAAACTAACGGGCTTCTcac   <  5:117822/139‑1 (MQ=255)
                                                                                                                tttCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCAGAGAGAAACTAACGGGCTTCTcac   >  5:172466/1‑138 (MQ=255)
                                                                                                                 ttCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTATCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCAGAGAGAAACTAACGGGCTTCTcaca  <  1:189780/138‑1 (MQ=255)
                                                                                                                         tCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCg                                                                   >  3:112912/1‑65 (MQ=255)
                                                                                                                         tCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCg                                                                   <  4:112912/65‑1 (MQ=255)
                                                                                                                                       |                                                                                                                    
TGGATCGACATATCCATAAACAATCCGCCTGACGTACGGACATAATCTATATTCGGGGGTTCTGGGTCCCGGGACGTGATTTTTAATAAATGAGGTGTACCGATTTCACCATTTTCAACAATCGTTTTGATTTTCTTAAAATGAGGATCAAAGCGGCGGTTAAATCCGACTTGAAGTGTTACCCCGTGTTTCCGGACAGCCGCAAGCGCTTCACTCGTTTCGTCCAGAGAGAAACTAACGGGCTTCTCACA  >  BS168_reference_genome/2839340‑2839590

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: