Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,511,930 G→T 29.9% S88* (TCA→TAA)  ykoA ← putative membrane protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,511,9300GT29.9% 70.9 / 13.7 54S88* (TCA→TAA) ykoAputative membrane protein of unknown function
Reads supporting (aligned to +/- strand):  ref base G (9/28);  new base T (15/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAACCTGTTTGTTCTCTTTTCCGCGCTTTTCTAATAGCCGTTTA  >  BS168_reference_genome/1511794‑1512052
                                                                                                                                         |                                                                                                                           
aGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTAATTTGa                                                                                                                            >  8:11401/1‑138 (MQ=255)
aGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TGTTATTTTGa                                                                                                                            >  7:171034/1‑138 (MQ=255)
 gCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑GTTTATTTTTaa                                                                                                                           >  7:95520/1‑138 (MQ=255)
  ccGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGG‑GTTTATTTTTaaa                                                                                                                          >  7:31721/1‑138 (MQ=255)
  ccGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGGTTAACGAAATGCGACAAACAAAATAAAAACGACTTGCTGGCCTAGGACAGCACGCGA‑TTTTAAATTGaaa                                                                                                                          >  8:120617/1‑138 (MQ=255)
      aaCGGGCTGGGACTGATCGCGGATTATCTAATTGTGGGAACTTCGATGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAaagaa                                                                                                                     <  7:349605/139‑1 (MQ=255)
      aaCGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAaagaa                                                                                                                     <  6:177287/139‑1 (MQ=255)
           gCTGGTACTTAGCGCGGAAGATCGAATTGTCGCCAAAGCTAATTTTGTCTTTTTCCCGTTTAACTAAATGCGTAAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGAc                                                                                                               <  7:252026/139‑1 (MQ=255)
            cTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑GTTTTTTTTGTAAAAGAAG‑CGGcc                                                                                                              >  5:96359/1‑139 (MQ=255)
            cTGGGACTGATCGCGGAAGATCGAACTGTCGGCACATCGAAGTTTGTCTTTTTACCGTTTAACGAAATTCGTCAAACAAAATAAAAAAGCCTTTCTGGCCTAGGACCGCAGGCGT‑TTTTTTTTTTTAAaataatcagc                                                                                                                >  4:251204/1‑133 (MQ=255)
                gACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGTGTTTTTTTTTTTAAAAAAAG‑‑GACcaaat                                                                                                          >  5:252277/1‑137 (MQ=255)
                gACTGAGCGCGGACTATCTAATTGGCGGCACATCGAAGTTTGTCTTTTTCCCGTTAAACCAAATGCGCAAAAAAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACcaatg                                                                                                          <  4:108741/139‑1 (MQ=255)
                 aCTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑GTTTATTTTGTAAAAGAAT‑CGCCcaatg                                                                                                          >  4:384821/1‑138 (MQ=255)
                  cttaTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACcaatgca                                                                                                        <  3:101916/136‑1 (MQ=255)
                  cTGATCGCGTAAGATCGAATTGTCGGCACATCGAAGTTTTTATTTTACCCTTTTAACTAAATGCGAAAAACAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACcaatgca                                                                                                        <  3:308512/139‑1 (MQ=255)
                  cTGATCGCGGAAGATCGAATTGTCGGCAATTCGAAGTTTGTATTTTTCCCGTTTAACCAAACGCGTCAAACAAAATAAACACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACcaatgca                                                                                                        <  7:333144/139‑1 (MQ=255)
                  cTGATCGCGGAAGACGGATTTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGACAAACAAAAAAAACACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACcaatgca                                                                                                        <  5:309738/139‑1 (MQ=255)
                  aggATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACTAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACcaatgca                                                                                                        <  1:342549/137‑1 (MQ=255)
                   tGATCGCGGAATTTCTAATTGGCGGAACATCGATGTTTGTTTTTTGCCCGTTTACCGAAATGCGAAAACAAAAATAAACCCGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACcaatgcaa                                                                                                       <  7:269195/139‑1 (MQ=255)
                     aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTTTTTTGTAAAATAAG‑CGACcaattcaatg                                                                                                     >  3:417603/1‑139 (MQ=255)
                     aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTTATTTTAAAAATAAA‑CGAAcaaatcaatg                                                                                                     >  2:30421/1‑139 (MQ=255)
                     aTCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑GTTTATTTTGAAAAAGAAG‑CCACcaatgccatg                                                                                                     >  7:293159/1‑139 (MQ=255)
                       cgctGAAGATCCAATTTGCGGCAAATCGATGTTTGTATTTTACCCGTTTAACTAAATGCTCAAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAatgcaatgca                                                                                                   <  3:94806/139‑1 (MQ=255)
                               aTCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTTTTTTTAAAAAAAAG‑CGACCAAAGCAAAGCATaaaaa                                                                                             >  2:342549/1‑136 (MQ=255)
                               aTCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TGTTTTTTTTTAAAAAAAG‑GCACCCAAGCAAAGGATAAAAACg                                                                                           >  1:207488/1‑139 (MQ=255)
                                   aaTTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGtgat                                                                                       <  7:110309/139‑1 (MQ=255)
                                        tgggCACATCTATGTTTTTCTTTTTCCCGTTTAACGAAATGCGTCAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAaca                                                                                  <  8:67849/137‑1 (MQ=255)
                                          ggCACATCGAATTTTGTCTTTTTCCCGTTTAACGAAATTCGTAAAAAAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTTAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACt                                                                                <  7:331254/139‑1 (MQ=255)
                                              cacccaaGTTTGTCTTTTACCCTTAAAACAAAAGCCTTCAAAAAAAAAAAAAACGCCTTGCGGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAACAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAAc                                                                            <  7:164534/134‑1 (MQ=255)
                                                   aaGTTTGTCTTTTAACCGTTAAACGAAATGCTTAAAACAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTTAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGata                                                                       <  8:95520/139‑1 (MQ=255)
                                                    tttttGTTTTTTACCGTTTTAACGAAATGCGTAAAAAAAAAAAAAAACGCCTTGGTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGata                                                                       <  8:388060/136‑1 (MQ=255)
                                                    attttGTTTTTTACCCGTTTCACGAATTCGTTCAAAAAAAAAAAAAACTCCTTGCTGGCCTATGACATCAGGCTT‑TTTTATTTTGAAAAAGAAG‑AGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGata                                                                       <  8:31721/136‑1 (MQ=255)
                                                      tttGTCTTTTTCCCTTTTAACTAAATGGTTAAAAAAAATTAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGataat                                                                     <  7:11401/138‑1 (MQ=255)
                                                      tttGTCTTTTTCCCGTTTAACTAAAGCCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCTT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGataat                                                                     <  8:21548/138‑1 (MQ=255)
                                                         gtgttttaCCCGTTAACCCAAATGCGCCAAACAAAAAAAACCCGCCTTGCGGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGc                                                                  <  3:425836/135‑1 (MQ=255)
                                                         gTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTTTTTTTTAAAAGAAGCCGCCCAATCCATTG‑AAAAAACCGAGCTGAACACAGATACGAGAATATcc                                                                 >  4:75346/1‑139 (MQ=255)
                                                         gTCTTTTACCCTTTAAACGAAATGCGTCAAAAAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGAAGAACACTGAACCGATAATAGc                                                                  <  1:91965/138‑1 (MQ=255)
                                                              tttCCCTTTAAACTAAATCCTTCAAAAAAAAAAAACACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGca                                                             <  3:384821/138‑1 (MQ=255)
                                                                 cccGTTTAACTAATTGCGAAAAAAAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCacac                                                          <  1:30421/138‑1 (MQ=255)
                                                                 cccGTTAAACGAAATGCTTCAAAAAAAAAAAAAACGCCTTGCTTGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGAAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCacac                                                          <  7:358277/138‑1 (MQ=255)
                                                                             aaTGCGAAAAAAAAAAAAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTg                                             <  5:388787/139‑1 (MQ=255)
                                                                               tGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTTTTTTTAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACAAGGAACCCATAATAGCCCGCCCACAGAAGAAAAGGTGaa                                           >  3:188859/1‑139 (MQ=255)
                                                                               tGCGTCAAAAAAAAAAAAAACGCATTGCTGGCCTAGGAAAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGaa                                           <  4:373833/139‑1 (MQ=255)
                                                                                   tCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTTTTTTTTAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACAATGAAGAAAAGGTGAACCtg                                       >  3:373833/1‑139 (MQ=255)
                                                                                           aaaaaacaCGCCTTGCTGTCATAGGACAGCAGGCTT‑TTTTATTTTGAAAAATAAG‑CGACAAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAACCTGTTTGTtct                               <  5:107616/136‑1 (MQ=255)
                                                                                                                  ggACAGCAGGCGT‑GTTTATTTTTAAAAAGCAG‑CGCACAATGCAATTCATAAAACCGTGCTGAACCCGGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAACCTGTTTGTTGGCGTTTCCGCGGGTTTCTAATAGc        >  1:132032/1‑139 (MQ=255)
                                                                                                                        cagGCGT‑TTTTATTTTTAAAAAGAAG‑CGCCCAAAGCAATGCATAAAAACCTGATGAACACTGCACCGATAAAATCCAGCACACGGAATAAAAGGTGAACCTGTTTGTTGTCTTTTACGCGCTTTTGTAATAGCCGTTTa  >  2:227521/1‑139 (MQ=255)
                                                                                                                        cagGCGT‑TTTTATTTTTAAAAAGAAG‑CGACCCAAGCCATGCATAAAACCCTGATGAACACAGAACCGATAATAGCCAGCACACAGAAGAAAACGTGTCCCTGTTTGTTGTCTATTACGCGCTTTTCTAATAGCCGGTTa  >  7:310063/1‑139 (MQ=255)
                                                                                                                        cagGCGT‑TTTTATTTTTAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAAcc                                         >  6:70175/1‑100 (MQ=255)
                                                                                                                        cagGCGT‑TTTTATTTTTAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAACCTGTTTGTTCTCTTTTCcg                       >  3:186516/1‑118 (MQ=255)
                                                                                                                        cagGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACCCTGAACCGATAATAGCCAGCCCACAGAAGAAAAGGTGAACCGGTTTGTTCTCTTTTGCGCGCTTTTCTAATAGCCGTTTa  >  6:365672/1‑139 (MQ=255)
                                                                                                                        cagGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAATGTGAACCTGTTTGTTCTCTTTTCCGCGCTTTTCTAATAGCCGTTTa  >  5:204794/1‑139 (MQ=255)
                                                                                                                        cagGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAAcc                                         <  5:70175/100‑1 (MQ=255)
                                                                                                                        cagGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAACCTGTTTGTTCTCTTTTCcg                       <  4:186516/118‑1 (MQ=255)
                                                                                                                                         |                                                                                                                           
AGCCGAAACGGGCTGGGACTGATCGCGGAAGATCGAATTGTCGGCACATCGAAGTTTGTCTTTTTCCCGTTTAACGAAATGCGTCAAACAAAATAAAAACGCCTTGCTGGCCTAGGACAGCAGGCGT‑TTTTATTTTGAAAAAGAAG‑CGACCAATGCAATGCATAAAACCGTGATGAACACTGAACCGATAATAGCCAGCACACTGAAGAAAAGGTGAACCTGTTTGTTCTCTTTTCCGCGCTTTTCTAATAGCCGTTTA  >  BS168_reference_genome/1511794‑1512052

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