Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,582,559 A→T 8.0% L290* (TTG→TAG) ‡ treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,582,5590AT 8.0% 143.4 / 3.3 50L290* (TTG→TAG) ‡treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base T (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.43e-01

CACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTC  >  NC_002947/4582415‑4582701
                                                                                                                                                 |                                                                                                                                              
cacaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGc                                                                                                                                            >  6:245802/1‑149 (MQ=255)
 acaGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGc                                                                                                                                            >  4:91592/1‑148 (MQ=255)
    ggTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCTTGGCCATg                                                                                                                                        <  1:281524/149‑1 (MQ=255)
     gTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATg                                                                                                                                        >  5:72279/1‑148 (MQ=255)
      tAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATGTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTa                                                                                                                                      >  7:176776/1‑149 (MQ=255)
       aGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCTTGGCCATg                                                                                                                                        <  8:4021/146‑1 (MQ=255)
       aGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCTTGGCCATg                                                                                                                                        >  7:4021/1‑146 (MQ=255)
           gcgcTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATgc                                                                                                                                 >  5:21834/1‑149 (MQ=255)
                 gcgGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGc                                                                                                                            <  5:245802/148‑1 (MQ=255)
                           aCCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGaaa                                                                                                                 >  5:172357/1‑149 (MQ=255)
                             cATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGCGAGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAg                                                                                                                >  1:254507/1‑149 (MQ=255)
                             cATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGt                                                                                                               <  4:188325/149‑1 (MQ=255)
                                     gggTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGcc                                                                                                       <  6:72279/149‑1 (MQ=255)
                                         cAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCatat                                                                                                   >  4:129346/1‑149 (MQ=255)
                                          aGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGGCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATg                                                                                                  <  1:330109/149‑1 (MQ=255)
                                           gTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATgg                                                                                                 >  6:226941/1‑149 (MQ=255)
                                                tcatGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATtc                                                                                            <  6:21834/149‑1 (MQ=255)
                                                  gtGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATtcgt                                                                                          >  6:138013/2‑149 (MQ=255)
                                                   tGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATtcgtc                                                                                         <  5:138013/149‑1 (MQ=255)
                                                    ggTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATtcgtcg                                                                                        >  1:123574/1‑149 (MQ=255)
                                                          gcAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgccgtcg                                                                                  <  5:226941/149‑1 (MQ=255)
                                                          gcAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgcc                                                                                      <  2:314551/145‑1 (MQ=255)
                                                          gcAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTcgtcgcc                                                                                      >  1:314551/1‑145 (MQ=255)
                                                            aaaaaGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCgtcgccgtcgcc                                                                                >  6:150765/1‑149 (MQ=255)
                                                            aaaaaGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCgtcgccgtcgcc                                                                                >  4:238304/1‑149 (MQ=255)
                                                             aaaaGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGTCATCAACGGAAAGTGGAAAGCCATATTGCATTCGTCGCCGTCGCCt                                                                               >  8:65439/1‑149 (MQ=255)
                                                              aaaGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCGATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGACGTCGCCtt                                                                              >  7:258991/1‑149 (MQ=255)
                                                                     gCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGaa                                                                       <  6:172357/149‑1 (MQ=255)
                                                                      cccACTGGC‑AGTTTGCCGGGATTTCCGGGGTATGTCTCAGGATGTCGTTTATAGGGAAGAGGACTTCCATGTCCAAGGCCATGTACATTCGCGTCATCAAAGGAAACTGGAAAGCCATATGGCATTCGTCGCCGTCTCCTTCGCCGAAg                                                                      >  4:218499/1‑149 (MQ=255)
                                                                      cccACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAg                                                                      <  6:8256/149‑1 (MQ=255)
                                                                             gC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCAAAGTACggg                                                                <  3:238304/148‑1 (MQ=255)
                                                                                            tttCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCCGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGc                                                 <  1:109733/149‑1 (MQ=255)
                                                                                                cGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACt                                             >  3:164276/1‑149 (MQ=255)
                                                                                                 ggggTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTg                                            <  8:114633/149‑1 (MQ=255)
                                                                                                    gTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGtt                                         <  2:254507/149‑1 (MQ=255)
                                                                                                    gTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGtt                                         <  2:177556/149‑1 (MQ=255)
                                                                                                    gTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGtt                                         <  4:46798/149‑1 (MQ=255)
                                                                                                            cAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCgg                                 <  5:150765/149‑1 (MQ=255)
                                                                                                               gATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGc                                >  5:120420/1‑147 (MQ=255)
                                                                                                               gATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCca                              <  2:100901/149‑1 (MQ=255)
                                                                                                                        gATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGc                     >  7:233065/1‑149 (MQ=255)
                                                                                                                           aGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGACCAGCAGCATGCGAt                  >  1:233456/1‑149 (MQ=255)
                                                                                                                                aGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCgggg             <  3:129346/149‑1 (MQ=255)
                                                                                                                                 gCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGt            <  2:123574/149‑1 (MQ=255)
                                                                                                                                     tCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAATTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGtt        <  8:118082/149‑1 (MQ=255)
                                                                                                                                     atgtacATGGCCTTGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGtt        >  2:82085/6‑149 (MQ=255)
                                                                                                                                         ccATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTg       <  6:120420/146‑1 (MQ=255)
                                                                                                                                           aTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTc  >  5:167110/1‑149 (MQ=255)
                                                                                                                                           aTGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTc  >  6:66057/1‑149 (MQ=255)
                                                                                                                                            tGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTc  >  4:142968/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CACAGGTAGTCGCGCTCGCGGTCGGTGACCATTTCCAGGGTCAGTTCATCATGGTTGCGCAAAAAGATCGCCCACTGGC‑AGTTGGCCGGGATTTCCGGGGTCTGGCGCAGGATGTCGGTGATAGGGAAGCGGTCTTCCATGGCCAAGGCCATGTACATGCGCGGCATCAACGGAAAGTGGAAAGCCATATGGCATTCGTCGCCGTCGCCTTCGCCGAAGTACGGGCGCGTGTCTTCAGGCCACTGGTTGGCTTCGGCCAGCAGCATGCGATCGGGGTAGTTGGCGTC  >  NC_002947/4582415‑4582701

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