Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 962,064 A→C 18.2% N224T (AAC→ACC)  ssuB → aliphatic sulfonate ABC transporter (ATP‑binding protein)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome962,0640AC18.2% 100.7 / 3.1 50N224T (AAC→ACC) ssuBaliphatic sulfonate ABC transporter (ATP‑binding protein)
Reads supporting (aligned to +/- strand):  ref base A (26/14);  new base C (2/8);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 1.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAGAAAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACTCATTGTA  >  BS168_reference_genome/961926‑962199
                                                                                                                                           |                                                                                                                                       
ggaAGAAAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAATCCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGaaa                                                                                                                                         >  1:210284/1‑139 (MQ=255)
ggaAGAAAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCTTAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGaaa                                                                                                                                         >  7:155197/1‑139 (MQ=255)
ggaAGAAAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGaaa                                                                                                                                         >  7:237119/1‑139 (MQ=255)
ggaAGAAAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGaaa                                                                                                                                         >  3:86223/1‑139 (MQ=255)
ggaAGAAAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGaaa                                                                                                                                         >  3:56777/1‑139 (MQ=255)
 gaagaaAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAAc                                                                                                                                        >  8:176060/1‑139 (MQ=255)
 gaagaaAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAAc                                                                                                                                        >  7:71887/1‑139 (MQ=255)
 gaagaaAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAAc                                                                                                                                        >  5:269566/1‑139 (MQ=255)
               cGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCAGCGAAACCGGACGACTCCGGa                                                                                                                          >  1:144248/1‑139 (MQ=255)
                         cgctaCTCATCATATACCTTAATCTTTTTATTTTTCAACTGCGCTGGATATCCTTACGCCCAACCCGGGAAAAACCATTAC‑GCTGCTCCCCATCCCCGTGGCGTAGCCCCGAAACCGCACGACTCCGGATTTCCAAGCg                                                                                                                <  4:27369/135‑1 (MQ=255)
                          gTGACTCATTATATATATGAATGTTTTTATTTAGTAAATGAGCTTGCGATCATTAAGGCCAGACCGGGTAAAATCCATAA‑GCTGATGCCCATTCACCTGGGGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCg                                                                                                                <  8:94015/138‑1 (MQ=255)
                           tGACTCATGATCTAGCTGAATCTTTTTATTTCGGAAATGCGCTTGCTATCCTTAAGGCCAACCCGGGAAAAATCCATAC‑GGTGATGCCCATTCACCGGGCTTACCCGCGAAACCGCACGACTCCGGATTTCCAAGCGAt                                                                                                              <  2:5658/139‑1 (MQ=255)
                           agacataTGTTATATATTAATTTTTTAATTTATCAAATTCGCTGGCTATCATTACGCCCAACCCGGGCAAAATCAAAAA‑GCGTAGCCACATTCACTTGGCGTAGCCGCCAACCCGCACGACTCCGGATTTCCAAGCGAt                                                                                                              <  6:359583/133‑1 (MQ=255)
                              cccATGATATACTTGATTCTGTTTATTTAGCAAATGAGCTGGCTAACCTACCCGCCACCCCGGGTAAAACCCATAC‑GCTGCGGCCCATCCCCCTGGCGTACCCCCCACCCCGCACGACTCCGGATTTCCAAGCGATcag                                                                                                           <  4:56777/137‑1 (MQ=255)
                                       tacaTGATTTTGTTTATTTAGTAATTGATCTTGCGAGCCTTACGCCCAACCCGTGGAAAATACTTAC‑GCTGAGGCCCATTCACGTGGCGTATCCCCCAACCGGGACGACTCCGGATTTCCAAGCGATCAGACagagag                                                                                                   <  8:155197/135‑1 (MQ=255)
                                        aGATGAATCTTTTAATTTAGCAAATGAGCTTGCGATCCTAACGCCCAACCCGGGCAAAACCCATAA‑GCTGAGCCCCATTCACCTGGCGTACCCCCCAACCCGCACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTg                                                                                                 <  2:108992/139‑1 (MQ=255)
                                            gAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACagagagctgtct                                                                                             >  5:359583/1‑133 (MQ=255)
                                             aaTCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAACCCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTc                                                                                            >  5:333100/1‑139 (MQ=255)
                                             aaTCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTc                                                                                            >  1:286710/1‑139 (MQ=255)
                                               tCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGAt                                                                                                              >  6:142411/1‑119 (MQ=255)
                                               tCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGAt                                                                                                              <  5:142411/119‑1 (MQ=255)
                                               tCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCgg                                                                                          >  4:359082/1‑139 (MQ=255)
                                               tCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCg                                                                                           >  4:319992/1‑138 (MQ=255)
                                               tCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCg                                                                                           >  1:333568/1‑138 (MQ=255)
                                               tCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCg                                                                                           >  3:277442/1‑138 (MQ=255)
                                               tCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCg                                                                                           >  6:61896/1‑138 (MQ=255)
                                                                  aGCTTGCGATCCTTAAGGACAAACCGGGGAAAATCCATAA‑CCTGATGCCCATTCACCTGGCGTATCCGCGAACCCGGACGACTCCGTATTTCCACGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACtgaagc                                                                       >  6:213686/1‑138 (MQ=255)
                                                                   gCTTGCGATCCTTACGCCAAACCCGCGGAAAATCCATAA‑GCTGATGCCCATGCACCTGGCTTACCCGCGAACCCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGAt                                                                      <  6:81985/139‑1 (MQ=255)
                                                                       gctaTCATTAAGGACACACCGGGTAAAATCCATAC‑GCTGATGCCCATTCCCCTGGCGTAACCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTg                                                                   <  6:96906/135‑1 (MQ=255)
                                                                        cgctcattACGGCCAAACCGGGCAAAATCCAAAA‑GCTGATCCCCATCCCCCTGGCGTAGCCGCTAACCCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGa                                                                 <  2:18942/133‑1 (MQ=255)
                                                                              ttAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACg                                                            >  6:157793/1‑138 (MQ=255)
                                                                                acgcccAACCCGGGCAAAATCCATAA‑GCGGAGCCCCATCCCCCTGGCGTACCCGCTAACCCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAg                                                         <  4:86223/135‑1 (MQ=255)
                                                                                 aGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAg                                                         >  3:90393/1‑138 (MQ=255)
                                                                                       aaCCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGAt                                                   >  2:212107/1‑138 (MQ=255)
                                                                                           ggggTAAAATCATAATGCTGAGCCCCTTTCCCTTGGCGTAGCCGCGACCCCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAgg                                               <  2:286710/139‑1 (MQ=255)
                                                                                                               aTGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGgcggaagc                           >  4:86043/1‑139 (MQ=255)
                                                                                                               aTGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGgcggaagc                           >  7:226405/1‑139 (MQ=255)
                                                                                                               aTGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGAGTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGgcggaagc                           >  4:1484/1‑139 (MQ=255)
                                                                                                                 gCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGgcggaagcg                          >  8:302906/1‑138 (MQ=255)
                                                                                                                           gtGGCGGCTCCGCCACACCGCCCGACGCCGGATTTCCAAGCGCTCAGACCGCGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGaaaaa               <  8:87918/138‑1 (MQ=255)
                                                                                                                             ggCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCa                                                 <  5:356190/103‑1 (MQ=255)
                                                                                                                             ggCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCa                                                 >  6:356190/1‑103 (MQ=255)
                                                                                                                             ggCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAg             >  3:250516/1‑139 (MQ=255)
                                                                                                                             ggCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAg             >  4:90845/1‑139 (MQ=255)
                                                                                                                              gcttATCGGCTAAACGGGACGACTCCGGATTTCCAAGCGATCAGACAGTGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAgg            <  2:144248/136‑1 (MQ=255)
                                                                                                                                 tATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACt         <  5:61896/139‑1 (MQ=255)
                                                                                                                                 tATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACt         <  5:157793/139‑1 (MQ=255)
                                                                                                                                  aaccGCTAAACCGCACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACt         <  3:320124/136‑1 (MQ=255)
                                                                                                                                  aTCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACt         <  3:319992/138‑1 (MQ=255)
                                                                                                                                   tCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTGCGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACt         >  6:349552/1‑137 (MQ=255)
                                                                                                                                   tCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTGCGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACt         <  5:349552/137‑1 (MQ=255)
                                                                                                                                   tCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACTCa       >  3:257700/1‑139 (MQ=255)
                                                                                                                                   cccccgACACCGCACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACt         <  6:269566/133‑1 (MQ=255)
                                                                                                                                        gAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACTCATTGTa  >  7:78452/1‑139 (MQ=255)
                                                                                                                                        gAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACTCATTGTa  >  5:1131/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GGAAGAAAAAAACAACGATGATTCTCGTGACTCATGATATAGATGAATCTGTTTATTTAGGAAATGAGCTTGCGATCCTTAAGGCCAAACCGGGGAAAATCCATAA‑GCTGATGCCCATTCACCTGGCGTATCCGCGAAACCGGACGACTCCGGATTTCCAAGCGATCAGACAGAGAGTGCTGTCGGAATTCGAGAAAACAGAAGATTTGGAGTACGCAGAAGGATCAGGGATATAGAAAGGGCGGAAGCGGAAATGAAAAAAGGACTCATTGTA  >  BS168_reference_genome/961926‑962199

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