Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,328,860 G→A 46.3% W657* (TGG→TGA sasA_18 → Adaptive‑response sensory‑kinase SasA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,328,8600GA46.3% 9.5 / 75.9 54W657* (TGG→TGAsasA_18Adaptive‑response sensory‑kinase SasA
Reads supporting (aligned to +/- strand):  ref base G (17/12);  new base A (13/12);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-02

CCTGGAAGACTACCACTCGCGCCTGACCGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAA  >  ACB143/5328725‑5328991
                                                                                                                                       |                                                                                                                                   
ccTGGAAGACTACCACTCGCGCCTGACCGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTgg                                                                                                                                     <  3:287893/136‑1 (MQ=255)
  tGGAAGACTACCACTCGCGCCTGACCGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGaa                                                                                                                                   >  5:151677/1‑136 (MQ=255)
    gAAGACTACCACTCGCGCCTGACCGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGaa                                                                                                                                   <  4:304523/134‑1 (MQ=255)
     aaGACTACCACTCGCGCCTGACCGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACaa                                                                                                                                <  1:82548/136‑1 (MQ=255)
             cACTCGCGCCTGACCGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATg                                                                                                                         >  7:74964/1‑135 (MQ=255)
             cACTCGCGCCTGACCGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATg                                                                                                                         >  5:198643/1‑135 (MQ=255)
                          ccGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCgg                                                                                                           <  5:160039/136‑1 (MQ=255)
                                 ggCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATcgcc                                                                                                    >  1:248763/1‑136 (MQ=255)
                                 ggCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATcgcc                                                                                                    >  5:159140/1‑136 (MQ=255)
                                 ggCCGCCGAGCTCGATGCCTTGCGCCGCGACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATcgcc                                                                                                    <  4:219778/136‑1 (MQ=255)
                                              gATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCAtgt                                                                                         <  6:159140/134‑1 (MQ=255)
                                              gATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCAtgt                                                                                         <  7:239994/134‑1 (MQ=255)
                                                  ccTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGGTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCg                                                                                    <  8:67407/135‑1 (MQ=255)
                                                   cgtGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGc                                                                                  <  7:244607/134‑1 (MQ=255)
                                                   cTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCggg                                                                                  >  6:178047/1‑135 (MQ=255)
                                                   cTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGc                                                                                  <  7:244604/136‑1 (MQ=255)
                                                   cTTGCGCCGCTACCATCGACAGACCCTGAAGGAGTTGCCCATGGGCGTATGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAATCATGTTGTCGGc                                                                                  <  5:269308/136‑1 (MQ=255)
                                                      gcgcCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATCTGGTCGGCTcg                                                                               >  6:75787/1‑136 (MQ=255)
                                                       cgccgcTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTcg                                                                               <  3:122331/135‑1 (MQ=255)
                                                       cgccgcTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTcg                                                                               <  3:117697/135‑1 (MQ=255)
                                                               ccATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTc                                                                      >  8:93349/1‑136 (MQ=255)
                                                                     ccAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATc                                                                <  6:264058/136‑1 (MQ=255)
                                                                        gACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATcgac                                                             >  8:140693/1‑136 (MQ=255)
                                                                            cTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCGCCATcgacgccc                                                         >  6:36971/1‑133 (MQ=255)
                                                                                aGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTg                                                     <  4:87326/136‑1 (MQ=255)
                                                                                aGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTg                                                     >  1:59509/1‑136 (MQ=255)
                                                                                  gAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGgc                                                   <  8:133831/136‑1 (MQ=255)
                                                                                     ttGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGcgcg                                                <  4:161278/136‑1 (MQ=255)
                                                                                     ttGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGcgcg                                                <  4:248085/136‑1 (MQ=255)
                                                                                     ttGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGcgcg                                                <  4:92875/136‑1 (MQ=255)
                                                                                       gCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGc                                              >  2:211076/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCctg                                           >  7:162970/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCctg                                           >  7:42404/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCctg                                           >  4:288948/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCctg                                           >  8:49578/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCctg                                           >  2:65842/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGATACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCAACGCCCAGCCTTGGCGCGGCctg                                           >  6:8797/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCCCCATCGACGAGCCTTGGCGCGGCctg                                           >  3:55684/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCctg                                           >  3:172457/1‑136 (MQ=255)
                                                                                          cATGGGCGTCTGCTCGCTGACCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGGGACCATCGCCGAGCCTTGGCGCGGccgg                                           >  2:276994/1‑134 (MQ=255)
                                                                                           aTGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCctgc                                          <  3:105620/136‑1 (MQ=255)
                                                                                                  tCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAgg                                   <  8:261686/136‑1 (MQ=255)
                                                                                                     gCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGACTCGAGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCtt                                >  7:304360/1‑136 (MQ=255)
                                                                                                     gCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCtt                                >  3:246360/1‑136 (MQ=255)
                                                                                                       tCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATgg                                                                                                                        <  4:113818/46‑1 (MQ=255)
                                                                                                       tCGCTGGCCAAGGACCAGGAAATCCTGATGTGAAACAAGGCCATgg                                                                                                                        >  3:113818/1‑46 (MQ=255)
                                                                                                         gCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTtcat                             <  2:179100/135‑1 (MQ=255)
                                                                                                                   gACCAGGAAATCCTGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCcg                  <  6:54664/136‑1 (MQ=255)
                                                                                                                          aaaTCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCACCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCa           >  7:205651/1‑136 (MQ=255)
                                                                                                                                tGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGca     >  8:40787/1‑136 (MQ=255)
                                                                                                                                tGATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGAATGTGGTCGGCTCGCTCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCAGca     >  5:59970/1‑136 (MQ=255)
                                                                                                                                 gATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCAGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGcac    >  1:161222/1‑136 (MQ=255)
                                                                                                                                 gATGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCAACGACGCGCCATTGCGCGGCCTGCTGCCAGGCTTAATCAACGTGCCCGCCGCGCACCTGcac    >  5:26341/1‑136 (MQ=255)
                                                                                                                                  aTGTGAAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGcac    <  1:25362/135‑1 (MQ=255)
                                                                                                                                   tgtgGAACAAGGCCATGGAAGAGATCACCGGAATAGACGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGAAAGGCTTCATCAACGTGCCCGCCGCGCCCCTGCACaa  >  2:135439/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                   
CCTGGAAGACTACCACTCGCGCCTGACCGGGCTGGCCGCCGAGCTCGATGCCTTGCGCCGCTACCATCGCCAGACCCTGCAGGAGTTGCCCATGGGCGTCTGCTCGCTGGCCAAGGACCAGGAAATCCTGATGTGGAACAAGGCCATGGAAGAGCTCACCGGCATCGCCGCCAAGCATGTGGTCGGCTCGCGCCTGGTCACCATCGACGAGCCTTGGCGCGGCCTGCTGCAAGGCTTCATCAACGTGCCCGACGAGCACCTGCACAA  >  ACB143/5328725‑5328991

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