Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,580,600 T→A 18.0% Q943L (CAG→CTG) ‡ treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,580,6000TA18.0% 103.2 / 17.5 50Q943L (CAG→CTG) ‡treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base T (23/18);  new base A (2/7);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.27e-01

TTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGAT  >  NC_002947/4580453‑4580739
                                                                                                                                                   |                                                                                                                                            
ttGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtg                                                                                                                                             >  4:107229/1‑149 (MQ=255)
ttGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtg                                                                                                                                             >  8:46431/1‑149 (MQ=255)
  gTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgc                                                                                                                                            <  5:223425/148‑1 (MQ=255)
   tAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTTCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgctt                                                                                                                                          <  7:204598/149‑1 (MQ=255)
    aaGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtgcttg                                                                                                                                         >  1:55946/1‑149 (MQ=255)
         cTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATAGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAg                                                                                                                                    >  6:194086/1‑149 (MQ=255)
           gtgtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAggg                                                                                                                                  >  5:287646/1‑149 (MQ=255)
           gtgtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAggg                                                                                                                                  <  2:243041/149‑1 (MQ=255)
            tgtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAgg                                                                                                                                   <  2:75836/147‑1 (MQ=255)
            tgtTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGACACCCTGAGCAAGCAGGCTCAgg                                                                                                                                   <  7:244820/147‑1 (MQ=255)
              ttttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGgtgt                                                                                                                               >  3:138639/1‑149 (MQ=255)
              ttttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAGGgtgg                                                                                                                               <  4:58585/149‑2 (MQ=255)
                  gCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGa                                                                                                                            <  3:107229/148‑1 (MQ=255)
                   cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGACACCCTGAGCAAGCAGGCTCAGGgtggcctgc                                                                                                                          <  2:215183/149‑7 (MQ=255)
                   cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGAATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAt                                                                                                                           <  4:249214/148‑1 (MQ=255)
                    ccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgt                                                                                                                         >  5:25919/1‑149 (MQ=255)
                     cGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgtg                                                                                                                        >  1:123436/1‑149 (MQ=255)
                        cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGTTAGGCATAACCCTGCACCACGAGCAACTGGACCAGCTGCAGGTCGCAGTGTACAAGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTg                                                                                                                     >  5:83953/1‑149 (MQ=255)
                               cAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTGGATGTGCTggg                                                                                                                   <  2:233675/144‑1 (MQ=255)
                               cAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAgg                                                                                                              >  6:296285/1‑149 (MQ=255)
                               cAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAgg                                                                                                                                   <  1:131109/128‑1 (MQ=255)
                               cAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAgg                                                                                                                                   >  2:131109/1‑128 (MQ=255)
                                 aGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGcc                                                                                                            <  4:138639/149‑1 (MQ=255)
                                               ttCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAgg                                                                                                                                   >  3:192858/1‑112 (MQ=255)
                                               ttCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAgg                                                                                                                                   <  4:192858/112‑1 (MQ=255)
                                                   ccTTCAAAATCGATGATGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGTGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAg                                                                                         <  5:296285/149‑1 (MQ=255)
                                                     ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGgt                                                                                       >  8:105721/1‑149 (MQ=255)
                                                                 gAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAg                                                                            <  2:55946/148‑1 (MQ=255)
                                                                          aTCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGa                                                                                 >  8:281061/1‑134 (MQ=255)
                                                                          aTCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGa                                                                                 <  7:281061/134‑1 (MQ=255)
                                                                                    accacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAg                                                        >  1:143930/1‑149 (MQ=255)
                                                                                    accacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAg                                                        >  1:203407/1‑149 (MQ=255)
                                                                                     ccacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                       <  5:194086/149‑1 (MQ=255)
                                                                                      cacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                   <  8:216566/73‑1 (MQ=255)
                                                                                      cacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                   >  7:216566/1‑73 (MQ=255)
                                                                                           gCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACt                                                                        <  3:108317/126‑1 (MQ=255)
                                                                                           gCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACt                                                                        >  4:108317/1‑126 (MQ=255)
                                                                                                   cccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTg                                         <  7:105721/149‑1 (MQ=255)
                                                                                                             ggTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCcagc                                  >  5:139667/1‑146 (MQ=255)
                                                                                                             ggTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCcagc                                  <  6:139667/146‑1 (MQ=255)
                                                                                                               tCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTc                             <  5:162732/149‑1 (MQ=255)
                                                                                                                  ccGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAg                          >  6:313059/1‑149 (MQ=255)
                                                                                                                  ccGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGGCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAg                          <  6:287646/149‑1 (MQ=255)
                                                                                                                       tACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGt                     >  2:98244/1‑149 (MQ=255)
                                                                                                                         cacGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑GGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTgg                   <  5:313059/149‑1 (MQ=255)
                                                                                                                             cATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtca               <  6:25919/149‑1 (MQ=255)
                                                                                                                             cATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGt                 >  4:34031/1‑147 (MQ=255)
                                                                                                                               tAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagt             <  8:279999/149‑1 (MQ=255)
                                                                                                                                  gCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagtt            >  5:226585/1‑147 (MQ=255)
                                                                                                                                  gCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagtt            >  2:128751/1‑147 (MQ=255)
                                                                                                                                  gCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagtt            >  2:133443/1‑147 (MQ=255)
                                                                                                                                       ccACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGc                                      >  1:278589/1‑116 (MQ=255)
                                                                                                                                       ccACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGc                                      <  2:278589/116‑1 (MQ=255)
                                                                                                                                        cACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCt     >  4:143402/1‑148 (MQ=255)
                                                                                                                                           ccTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTAGGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGAt  >  1:234624/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGT‑TGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGAT  >  NC_002947/4580453‑4580739

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