Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,413,448 T→G 25.7% T198P (ACC→CCC)  yvrO ← putative ABC transporter (ATP‑binding protein)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,413,4480TG25.7% 49.0 / 11.6 38T198P (ACC→CCC) yvrOputative ABC transporter (ATP‑binding protein)
Reads supporting (aligned to +/- strand):  ref base T (6/22);  new base G (8/2);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 1.86e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTATGAGAGAAAATGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACTTGATTTGGTATCAAGCGCCCCGGTCGGTTCATCAGCTAAAATAAACCTCGGCTGATTGACAAGGGCCCTCG  >  BS168_reference_genome/3413310‑3413582
                                                                                                                                          |                                                                                                                                       
tttATGAGAGAAAATGGAGTTAAAGGCTATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGt                                                                                                                                         <  6:280538/139‑1 (MQ=255)
    tGAGAGAAAATGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAAAACCAACTCTCCGTCTCTTATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTAt                                                                                                                                      <  4:154322/138‑1 (MQ=255)
      agacaAAATGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCa                                                                                                                                    <  5:222123/138‑1 (MQ=255)
          aaaaTGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTAt                                                                                                                                      <  2:93073/132‑1 (MQ=255)
          aaaaTGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCataat                                                                                                                                >  3:155187/1‑138 (MQ=255)
          aaaaTGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGGTTTCctaat                                                                                                                                >  6:238800/1‑138 (MQ=255)
          aaaaTGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGGTAt                                                                                                                                      >  1:93073/1‑132 (MQ=255)
           aaaTGGATTTATAGGAGATGCTAATATGTTCAGGCATGGCTTATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTGCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCAtaataa                                                                                                                              <  4:72960/139‑1 (MQ=255)
           aaaTGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCataata                                                                                                                               <  6:165827/138‑1 (MQ=255)
           aaaTGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCAtaataa                                                                                                                              <  4:425198/139‑1 (MQ=255)
           aaaTGGATTGAAAGGAGATTCGAATATGTACAAGCATGGCTGATAACCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCAtaataa                                                                                                                              <  2:183353/139‑1 (MQ=255)
                  ttGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCtt                                                                                                                       <  6:403012/139‑1 (MQ=255)
                    gAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGGTTTGATTATTATGGTGTTGc                                                                                                                     >  4:166798/1‑139 (MQ=255)
                     aaaGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGc                                                                                                                     <  1:57661/138‑1 (MQ=255)
                     aaaGGAGATTCGAATATGGTCAAGCATGGCTGATAACCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGc                                                                                                                     <  4:375803/138‑1 (MQ=255)
                               cGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATg                                                                                                          <  3:398170/138‑1 (MQ=255)
                               cGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATg                                                                                                          <  1:327120/138‑1 (MQ=255)
                                gattatGTTGACGCTTGGCTGATTTCCCCTCTGACATCTAAAACCACCTCGCGTGGTCTGATAAAATCGGTTCTGGCAGCTTTGTTGGGGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCa                                                                                                        <  4:39890/136‑1 (MQ=255)
                                       ttaaaGCTTGGCTGATATCCCCTCTCTCATATAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGa                                                                                                 <  6:30628/136‑1 (MQ=255)
                                         cAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGaa                                                                                                <  4:166887/138‑1 (MQ=255)
                                                 gCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGGTCTGTCAGCGTTTTTTGCGACGTCCGGGTCGTGGGTTTTAATATTAATTTTATTTCATTACCGGATCCA‑TTCAGAAAACAATga                                                                                        >  4:79874/1‑137 (MQ=255)
                                                       atCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATGATTTTAATTGTCTTGCCTT‑CCCGATGCAGGTCAGAAAACAGTGCCAAAAtt                                                                                 >  8:333717/1‑139 (MQ=255)
                                                            ctctctCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGGCGTGGGGTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTc                                                                            <  6:102000/139‑1 (MQ=255)
                                                            ctctctCATCTAATACAAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTTTCATATTAATTGTCTTTCCGT‑CCCGAGGGAGTTATGAAAACAGTTACCAAATTTGtt                                                                             >  7:89308/1‑138 (MQ=255)
                                                                 tcATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACt                                                                         <  3:168475/137‑1 (MQ=255)
                                                                    tCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGGTTTCATAATAATTGTCTTGGCGT‑CCCGATGCAGTGCAGAAAACACTGCCAAAAATTTTTGACTTGttt                                                                    >  7:127311/1‑139 (MQ=255)
                                                                         tACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGGTTTCATAAAAATTGTGTTTGCTT‑CCCGCTGCAGTTCAGAAAAAAGTGCCAAAAATTTTTTACTTGATTTTGtt                                                               >  5:227270/1‑138 (MQ=255)
                                                                               ctctCCGTCTCCGATAAAAACGGTTCTGTCAGCTTTTTTTGCGACTTCCGGGTCTTGGGTTTTAATTATAATGGTATTGCCAT‑GACGAGTCAGTTCAGCAAAAAGTTCCAAAATTTGTTAAATGGAGTTGGTATCAAgc                                                         >  6:74489/1‑139 (MQ=255)
                                                                                   ccGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGGTTTGATTATAAATGTTTTGGCGT‑CCCGCTGCAGGTCAGAAACAACTGGCACAATTTTTTCAATTGCTTTTGGTTGAAGCGccc                                                     >  5:13341/1‑139 (MQ=255)
                                                                                   ccGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGGTATCATAATTATTGTCTTGCCGT‑GCCGCTGCAGTTCAGAAAACAGTGCCAAAATTTTTTGACTTGATTTGGTATGAAGCGccc                                                     >  1:237883/1‑139 (MQ=255)
                                                                                   ccGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGGTATCATAATAATGGTGTTGCCTT‑GCCGATGCAGTTGAGAAAACAGTGCCAAAATTTGTTCACTTTGTTTGGTTTCAAGCGccc                                                     >  8:91324/1‑139 (MQ=255)
                                                                                       ctctGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACTTGATTTGGTATCAAGCGCCCcggt                                                 <  3:166798/139‑1 (MQ=255)
                                                                                          tGATAAAAACGGTTCTGTCAGCTTTTTTGGCGCCGTCCGGGTCGGGGGTTATCATTATGATGGTCTTGCGTT‑GCCGAGGCAGGTCAGACAACAGTGGCAAAAATTTTTTACGTGATTTGGTATCAAGCGCCCcggtcgg                                              >  8:304577/1‑139 (MQ=255)
                                                                                                                    ttGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACTTGATTTGGTATCAAGCGCCCCGGTCGGTTCATCAGCTAAAATAAACCTCGGc                     <  4:155187/138‑1 (MQ=255)
                                                                                                                       gCGACGTCCGGGTCGTGGGGTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGGTCACTTGATTTGGTATCAAGCGCCCCGGTCGGTTCATCAGCTAAAATAAACCTCGGCTGAt                 <  6:13341/139‑1 (MQ=255)
                                                                                                                        cagcGTCCGGGTCGGGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACTTGATTTGGTATCAAGCGCCCCGGTCGGTTCATCAGCTAAAATAAACCTCGGCTGAtt                <  7:45472/136‑1 (MQ=255)
                                                                                                                             tcagggTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACTTGATTTGGTATCAAGCGCCCCGGTCGGTTCATCAGCTAAAATAAACATCGGCTGATTGACaa           <  5:74489/136‑1 (MQ=255)
                                                                                                                                  gTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACTTGATTTGGTATCAAGCGCCCCGGTCGGTTCATCAGCTAAAATAAACCTCGGCTGATTGACAAGGGcc      <  6:227270/139‑1 (MQ=255)
                                                                                                                                      tGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACTTGATTTGGTATCAAGCGCCCCGGTCGGTTCATCAGCTAAAATAAACCTCGGCTGATTGACAAGGGCCCTCg  <  2:237883/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
TTTATGAGAGAAAATGGATTGAAAGGAGATTCGAATATGTTCAAGCATGGCTGATATCCCCTCTCTCATCTAATACCAACTCTCCGTCTCTGATAAAAACGGTTCTGTCAGCTTTTTTGGCGACGTCCGGGTCGTGGGTTATCATAATAATGGTCTTGCCTT‑CCCGATGCAGTTCAGAAAACAGTGCCAAAATTTGTTCACTTGATTTGGTATCAAGCGCCCCGGTCGGTTCATCAGCTAAAATAAACCTCGGCTGATTGACAAGGGCCCTCG  >  BS168_reference_genome/3413310‑3413582

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