Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,644,663 C→T 7.0% D212N (GAC→AAC)  degU ← two‑component response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,644,6630CT 7.0% 182.0 / 2.9 57D212N (GAC→AAC) degUtwo‑component response regulator
Reads supporting (aligned to +/- strand):  ref base C (24/29);  new base T (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.67e-01

ACTCTATTGTAAGGGAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTGCAGCACTTCACATTCCCGCCTAGTTA  >  BS168_reference_genome/3644539‑3644800
                                                                                                                            |                                                                                                                                             
aCTCTATTGTAAGGGAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTAcattcatt                                                                                                                                 <  2:30583/138‑1 (MQ=255)
    tATTGTAAGGGAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTCAATAAttttttt                                                                                                                            >  2:67912/1‑138 (MQ=255)
      ttGTAAGGGAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTAGATTGAtttttttt                                                                                                                           >  2:43310/1‑136 (MQ=255)
           aGGGAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCTTTTCTACCCAGCCATTTTTAATGGCGCCCACAACGGCTTGCGTCCGGTCGTGGAGAGTGATTTTTTGGAAtata                                                                                                                     >  2:266892/1‑139 (MQ=255)
           aGGGAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTCATTAATTTTTTTTTAAat                                                                                                                      >  5:58674/1‑138 (MQ=255)
           aGGGAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTCAATGATTTTTTTTTAAata                                                                                                                     >  2:232510/1‑139 (MQ=255)
              gAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTTGTTTACATTCATTTTTTGTAAAATATTg                                                                                                                  >  6:29405/1‑139 (MQ=255)
              gAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTAAAATCATTTTTTTTTAAatatat                                                                                                                  >  6:308577/1‑137 (MQ=255)
              gAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGGTCATTGATTTTTTTTTAAATAttt                                                                                                                  >  4:219400/1‑138 (MQ=255)
              gAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGGCCGGTTGTTTACATTCATTTTTTGTAAAATATTg                                                                                                                  >  6:33634/1‑139 (MQ=255)
              gAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTCCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTAATTAATTTTTTTTTAAATAtttg                                                                                                                 >  8:218296/1‑137 (MQ=255)
                aaCGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTTCt                                                                                                                >  6:66804/1‑139 (MQ=255)
                aaCGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTCAATAATTTTTTTTTAAATAtttca                                                                                                                >  6:83654/1‑136 (MQ=255)
                   ggCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAATCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGAc                                                                                                            <  5:135462/139‑1 (MQ=255)
                        ccTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTACCTTGTTTTTTTTTTAAATATTTATGTCATGGt                                                                                                     >  1:26081/1‑139 (MQ=255)
                               aaaGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAAc                                                                                              <  1:232510/139‑1 (MQ=255)
                                       acagaGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTtct                                                                                      <  1:208109/135‑1 (MQ=255)
                                        ccTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTtct                                                                                      <  1:148303/138‑1 (MQ=255)
                                        ccTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTtct                                                                                      <  1:67912/138‑1 (MQ=255)
                                        ccTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTtct                                                                                      <  1:43310/138‑1 (MQ=255)
                                         cTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTTGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCa                                                                                    <  5:33634/139‑1 (MQ=255)
                                         cTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTTGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCa                                                                                    <  5:29405/139‑1 (MQ=255)
                                         cTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCa                                                                                    <  5:308577/139‑1 (MQ=255)
                                            gTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTTCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTg                                                                                 <  2:182212/139‑1 (MQ=255)
                                               aaaCGCAAGTCTACTATTATACTAAAGCCTTTCTACCCAGCCATTTTTAAAGGCGACCACAACGGCGTTCGGCCGGGCGTTTAAAATCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATa                                                                              <  1:50783/139‑1 (MQ=255)
                                                 aGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATaaa                                                                            <  7:271066/139‑1 (MQ=255)
                                                      aGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTAGATTCATTTTTTTTAAAAAATTTGTGGCAAGGTTTTTTACGGTTTTTTGACGGATAAACaggg                                                                       >  2:208109/1‑136 (MQ=255)
                                                       gTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTTCATTAATTTTTTTTTAAATATTTGTTACATGGTTTTTTAAGGGTTTGTTACAGATaaaaaaata                                                                      >  6:251283/1‑136 (MQ=255)
                                                       gTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGa                                                                      <  6:97118/139‑1 (MQ=255)
                                                       gTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGa                                                                      <  4:39903/139‑1 (MQ=255)
                                                       gTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTACATTAAATTTTTTTAAAATATTTGGGTCAATGTTTTTTAACGGTTTGTCACAGATAAAAAATGa                                                                      >  2:243488/1‑139 (MQ=255)
                                                            cTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCAc                                                                 <  5:83654/139‑1 (MQ=255)
                                                            cTATTATACTATCTAATTTCTAGCCAGCGATTTTTTATGGCGACCACAACGGCTTGGGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTAAc                                                                 <  7:105579/139‑1 (MQ=255)
                                                                  taCTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGGTGATTCATTTTTTTTTAAATATTGGGGACCAGGGTTTTTAAGGTTTTGTGCCTGTTTAACCATGTTTCAcccaaaa                                                           >  4:64700/1‑134 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTTCATTAATTTTTTGTATAATATTTGTGACATGGTTTTTTAAGGTTTTGTCACTGATAAAAAATTTTTCACCAATTCCGCgg                                                      >  4:209704/1‑139 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACcgcg                                                       <  6:84272/138‑1 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTCATTAATTTTTTTTTTAATATTTGTGACATTGGTTTTTAAGGTTTTATAACTTATAAAAAAAAATTCACCAATACCGCgg                                                      >  5:65093/1‑139 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTCATTAATTTTTTTTTAAATATTTGTGACATTGTTTTTAAAGGGTTTGTCACGGATAAAAAATGATTGACCAAAAcaccgg                                                      >  8:102342/1‑134 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTCATTAATTTTTTTTTAAATATTTATGACAATGTTTTTTAAGGTTTTATAAATGATAAACAATTATTAACCAAAACCGCgg                                                      >  8:230816/1‑139 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTTCAATCATTTTTTGTAAAATATTTCTGGCATTGTTTTTTACGGTTTTATCACTGATAAAAAATTATTAACCAATAAcgccg                                                      >  7:20568/1‑137 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGTACATTAAATTTTTTTAAAAAATTGCTGACATTGGTTTTAACGGGTTTATGACTGATAAAAAAAGATTACCCAATTCCCCgg                                                      >  2:10863/1‑139 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTGGTCATTTAATTTTTTTTAAATATTGGTGGCATGGGTTTTTAGGGTTTTGTGACTGTGAAACAAAGTTTGCCAAATACCgggg                                                      >  5:271206/1‑139 (MQ=255)
                                                                       tctcATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGGCCGGTCGTGTACAGTCGTTTTTTGTAAAATATTGGGGCCGTGGTGTGTGAGGGTTTTGGCACTGGTAAAGGGTGGTTGAGCAAGACCGCgg                                                      >  4:244164/1‑139 (MQ=255)
                                                                                 aCCCAGCCATTTTTAATGGCGACCACACCGGCTTGCGTCCGGTCGTGGAGCTTGAATTTTTTTTAAAAATTTGGTCCATGGGTTTTAACGGGGTTGTACCGGTTTAAAAAAGATTGACCAATAACCCGGGTGCTTTTTc                                            >  2:194372/1‑139 (MQ=255)
                                                                                         aTTTTTTATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCa                                    <  1:267439/139‑1 (MQ=255)
                                                                                         aTTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCa                                    <  1:266892/139‑1 (MQ=255)
                                                                                         aTTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCa                                    <  4:309566/139‑1 (MQ=255)
                                                                                         aTTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGACAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCa                                    <  3:13478/139‑1 (MQ=255)
                                                                                             ttAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGc                                 <  2:280262/138‑1 (MQ=255)
                                                                                             ttAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGc                                 <  5:1279/138‑1 (MQ=255)
                                                                                                      aCCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTgcagca                       <  5:251283/139‑1 (MQ=255)
                                                                                                           aaCTGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTGCAGCACTTca                  <  3:36/139‑1 (MQ=255)
                                                                                                                     gTCCGGTCGTTTACATTGATTTTTTGTATAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTGCAGCACTTCACATTCCCGcc        <  1:157386/139‑1 (MQ=255)
                                                                                                                     gTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTGCAGCACTTCACATTCCCGcc        <  6:65093/139‑1 (MQ=255)
                                                                                                                           tCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTGCAGCACTTCACATTCCCGCCTAGTTa  <  5:66804/139‑1 (MQ=255)
                                                                                                                           tCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTGCAGCACTTCACATTCCCGCCTAGTTa  <  6:68296/139‑1 (MQ=255)
                                                                                                                           tCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTGCAGCACTTCACATTCCCGCCTAGTTa  <  1:243488/139‑1 (MQ=255)
                                                                                                                            |                                                                                                                                             
ACTCTATTGTAAGGGAAACGGCAGCCTAAAAAAAGACCTGCCTAGTAAAAGGCAAGTCTCCTATTATACTATCTCATTTCTACCCAGCCATTTTTAATGGCGACCACAACGGCTTGCGTCCGGTCGTTTACATTCATTTTTTGTAAAATATTGCTGACATGGTTTTTAACGGTTTTCTCACTGATAAACAATGATTCACCAATACCGCGGTTGCTTTTTCCGTCTGCAAGCATCTGCAGCACTTCACATTCCCGCCTAGTTA  >  BS168_reference_genome/3644539‑3644800

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: