Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,581,699 C→T 8.5% L591L (CTG→CTA)  treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,581,6990CT 8.5% 128.4 / 2.7 47L591L (CTG→CTA) treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base T (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

ACACCGGCGTGGTGGCGGTGCCGAAACGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAGGTGGCAAGGTCAGCAGAAACGGCA  >  NC_002947/4581554‑4581847
                                                                                                                                                 |                                                                                                                                                    
acacCGGCGTGGTGGCGGTGCCGAAACGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTgcagca                                                                                                                                                  >  1:7867/1‑150 (MQ=255)
    cGGCGTGGTGGCGGTGCCGAAACGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGttc                                                                                                                                              >  1:245557/1‑150 (MQ=255)
          ggtggCGGTGCCGAATCTCCCCCCGTAGCACAGGCCCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCAt                                                                                                                                        <  2:373452/150‑1 (MQ=255)
           gtggCGGTGCCGAAACGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGTAGCAGTTCTTCCATg                                                                                                                                       >  2:53684/1‑150 (MQ=255)
           gtggCGGTGCCGAAACGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGTAGCAGTTCTTCCATg                                                                                                                                       >  2:53526/1‑150 (MQ=255)
              gCGGTGCCAGACCGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATgcgc                                                                                                                                    <  2:613080/150‑1 (MQ=255)
                     cGAAACGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCa                                                                                                                             >  2:463873/1‑150 (MQ=255)
                      gAAACGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAg                                                                                                                            >  2:154283/1‑150 (MQ=255)
                            cacgcCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCaa                                                                                                                      >  1:399618/1‑150 (MQ=255)
                                         aGGCGCGCATCAGCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCg                                                                                                         <  1:46856/150‑1 (MQ=255)
                                          ggCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGTAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCgg                                                                                                        <  1:53526/150‑1 (MQ=255)
                                          ggCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGTAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCgg                                                                                                        <  1:53684/150‑1 (MQ=255)
                                                    atcgatcgGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTc                                                                                              >  1:1049374/1‑150 (MQ=255)
                                                             cccTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTg                                                                                     >  1:696804/1‑150 (MQ=255)
                                                               cttccttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGa                                                                                   <  1:1037439/150‑1 (MQ=255)
                                                               cttccttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGa                                                                                   <  1:691367/150‑1 (MQ=255)
                                                                ttccttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGttcttc                                                                                                                                           <  2:940992/93‑1 (MQ=255)
                                                                ttccttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGttcttc                                                                                                                                           >  1:940992/1‑93 (MQ=255)
                                                                ttccttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGAt                                                                                  <  1:85499/150‑1 (MQ=255)
                                                                ttccttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGAt                                                                                  >  1:797486/1‑150 (MQ=255)
                                                                    ttGCTGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGc                                                                              <  2:499114/150‑1 (MQ=255)
                                                                    ttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGc                                                                              >  1:132198/1‑150 (MQ=255)
                                                                    ttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGc                                                                              >  1:179143/1‑150 (MQ=255)
                                                                    ttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGc                                                                              >  1:460813/1‑150 (MQ=255)
                                                                    ttGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCACAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGc                                                                              <  1:205353/150‑1 (MQ=255)
                                                                         ctctaaCCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCt                                                                         >  1:588252/5‑150 (MQ=255)
                                                                               ccAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCAt                                                                   >  1:704172/1‑150 (MQ=255)
                                                                                 aGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCNTGc                                                                 <  2:696804/150‑1 (MQ=255)
                                                                                            ggCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGcct                                                      <  2:994383/150‑3 (MQ=255)
                                                                                            ggCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTg                                                      <  2:132198/150‑1 (MQ=255)
                                                                                             gCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCACTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGAGGCTgg                                                     >  1:211633/1‑150 (MQ=255)
                                                                                              cAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGc                                                    <  2:714128/150‑1 (MQ=255)
                                                                                                 gTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGACAGATGGGGATGCGATCGCGGCTGGCcag                                                 >  2:362106/1‑150 (MQ=255)
                                                                                                           aGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAg                                       <  1:939274/150‑1 (MQ=255)
                                                                                                              aTGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCCGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAg                                    <  1:304408/150‑1 (MQ=255)
                                                                                                                     tttGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCAt                             >  2:75608/1‑150 (MQ=255)
                                                                                                                       tGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGCCTgggc                                                                                  <  2:792174/95‑3 (MQ=255)
                                                                                                                       tGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGCCTgggc                                                                                  >  1:792174/1‑93 (MQ=255)
                                                                                                                         cAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAgg                         <  2:797486/150‑1 (MQ=255)
                                                                                                                              tATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAGGTGGCa                    >  1:699438/1‑150 (MQ=255)
                                                                                                                                tCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTc                                                                                              <  2:943447/74‑1 (MQ=255)
                                                                                                                                tCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTc                                                                                              >  1:943447/1‑74 (MQ=255)
                                                                                                                                   cTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAGGTGGCAAGGTc               <  1:1127799/150‑1 (MQ=255)
                                                                                                                                        gggggCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAAc                                          >  2:643278/1‑118 (MQ=255)
                                                                                                                                        gggggCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAAc                                          <  1:643278/118‑1 (MQ=255)
                                                                                                                                          gggCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAGGTGGCAAGGTCAGCAGaa        >  2:321614/1‑150 (MQ=255)
                                                                                                                                          gggCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAGGTGGCAAGGTCAGCAGaa        >  1:292962/1‑150 (MQ=255)
                                                                                                                                              ttGCAGCAGTTCTTCCGTGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAGGTGGCAAGGTCAGCAGAAACgg    >  2:711351/1‑150 (MQ=255)
                                                                                                                                                gcagcaGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAGGTGGCAAGGTCAGCAGAAACGGCa  >  2:646395/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ACACCGGCGTGGTGGCGGTGCCGAAACGCACGCCGTAGCACAGGCGCACGGCATCGATCGGCCCTTCCTTGCCGGCGAACCAGCGGCGCTTGGGCAGGTACTGCGGCAGGATGGTGCTTTGCAGCGTATCGCTTGCGGGGGCTTGCAGCAGTTCTTCCATGCGCTTGCGCAGTACCAATGTGGTCAATTCGGGTAACCCCTCGGTGGGCTGGATGTGCCAGCTGGGCATGCGATCGCGGCTGGCCAGCAGGAACCAGTAGAACGCATAAGGTGGCAAGGTCAGCAGAAACGGCA  >  NC_002947/4581554‑4581847

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