Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,572,399 T→A 8.1% V740D (GTC→GAC)  treY → maltooligosyl trehalose synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,572,3990TA 8.1% 184.6 / 5.2 59V740D (GTC→GAC) treYmaltooligosyl trehalose synthase
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base A (3/2);  total (33/26)
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.00e+00

CTGGAGCGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGCACGCTGGACGATGCCACGCCGGCC  >  NC_002947/4572253‑4572547
                                                                                                                                                  |                                                                                                                                                    
cTGGAGCGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAAGGGCCTGGACCa                                                                                                                                                   >  1:1069685/1‑150 (MQ=255)
         gcCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGtt                                                                                                                                          >  2:953009/1‑150 (MQ=255)
         gcCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGtt                                                                                                                                          <  1:1192568/150‑1 (MQ=255)
         gcCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGtt                                                                                                                                          <  2:548910/150‑1 (MQ=255)
           cGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAAGGGCCTGGACCAGGccattgaga                                                                                                                                        <  2:1069685/150‑8 (MQ=255)
                 aGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGTGGACCATCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGa                                                                                                                                  <  1:317765/150‑1 (MQ=255)
                 aGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGTGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGa                                                                                                                                  <  2:791089/150‑1 (MQ=255)
                             gcgcCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGcc                                                                                                                                               <  2:422391/125‑1 (MQ=255)
                             gcgcCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGcc                                                                                                                                               >  1:422391/1‑125 (MQ=255)
                                  tGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGccc                                                                                                                 <  1:906257/150‑1 (MQ=255)
                                       tCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCt                                                                                                            >  2:1125844/1‑150 (MQ=255)
                                          ttACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTa                                                                                                         >  1:620462/1‑150 (MQ=255)
                                             ccTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGACCAGGccattga                                                                                                                                          >  2:516781/1‑109 (MQ=255)
                                             ccTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGACCAGGccattga                                                                                                                                          <  1:516781/114‑6 (MQ=255)
                                              cTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGgt                                                                                                                       >  1:496029/1‑132 (MQ=255)
                                              cTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGgt                                                                                                                       <  2:496029/132‑1 (MQ=255)
                                              cTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAg                                                                                                     <  1:558404/150‑1 (MQ=255)
                                              cTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAg                                                                                                     <  2:109437/150‑1 (MQ=255)
                                                                gACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGAACGGGCGCAACGTGCTGTGGACGGCACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGAc                                                                                   >  1:444903/1‑150 (MQ=255)
                                                                    gCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCa                                                                               <  2:1066420/150‑1 (MQ=255)
                                                                     cGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAg                                                                              >  2:443234/1‑150 (MQ=255)
                                                                      gagaACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGc                                                                             <  2:311099/150‑1 (MQ=255)
                                                                      gagaACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGc                                                                             >  1:392423/1‑150 (MQ=255)
                                                                            cAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTc                                                                       >  2:598593/1‑150 (MQ=255)
                                                                            cAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTc                                                                       <  1:953009/150‑1 (MQ=255)
                                                                             aGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCg                                                                      >  2:1015801/1‑150 (MQ=255)
                                                                              gCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGa                                                                     >  2:1190171/1‑150 (MQ=255)
                                                                                aaCTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGAcc                                                                   >  2:933867/1‑150 (MQ=255)
                                                                                 aCTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGAccc                                                                  <  2:988742/150‑1 (MQ=255)
                                                                                      gcAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGccc                                                                                                                 <  1:809558/98‑1 (MQ=255)
                                                                                      gcAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGccc                                                                                                                 >  2:809558/1‑98 (MQ=255)
                                                                                         aGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAAcc                                                          <  1:1125844/150‑1 (MQ=255)
                                                                                             acTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAAGGAATACTGGGACTTCAGCCTGGTCGACCCGGACAAccgccg                                                      >  2:825715/2‑150 (MQ=255)
                                                                                              cTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACcgccgc                                                     >  2:1078485/1‑150 (MQ=255)
                                                                                               tGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACcgccgcg                                                    >  2:534477/1‑150 (MQ=255)
                                                                                               tGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACcgccgcg                                                    >  1:588286/1‑150 (MQ=255)
                                                                                                ggCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACcgccgcgc                                                   >  1:879647/1‑150 (MQ=255)
                                                                                                 gCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACcgccgcgcc                                                  <  1:692236/150‑1 (MQ=255)
                                                                                                   cGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGt                                                <  2:750997/150‑1 (MQ=255)
                                                                                                   cGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGt                                                >  2:72442/1‑150 (MQ=255)
                                                                                                   cGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGt                                                >  1:1197110/1‑150 (MQ=255)
                                                                                                   cGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGt                                                <  1:598593/150‑1 (MQ=255)
                                                                                                            gcAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGc                                       >  1:37416/1‑150 (MQ=255)
                                                                                                                  gctgGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGgtgt                                                                                                                     <  1:777916/66‑1 (MQ=255)
                                                                                                                  gctgGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGgtgt                                                                                                                     >  2:777916/1‑66 (MQ=255)
                                                                                                                      gCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGt                             >  2:557956/1‑150 (MQ=255)
                                                                                                                         tGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGc                          <  1:443234/150‑1 (MQ=255)
                                                                                                                              cGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTg                                                                                       <  2:802680/84‑1 (MQ=255)
                                                                                                                              cGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTg                                                                                       >  1:802680/1‑84 (MQ=255)
                                                                                                                                  gcTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTgg                                                                                      >  1:250775/1‑81 (MQ=255)
                                                                                                                                  gcTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTgg                                                                                      <  2:250775/81‑1 (MQ=255)
                                                                                                                                    tctcAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTAcc                                                                                                       <  2:807191/62‑1 (MQ=255)
                                                                                                                                    tctcAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTAcc                                                                                                       >  1:807191/1‑62 (MQ=255)
                                                                                                                                     ctcAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGCACGCTGGACGAt              >  1:662510/1‑150 (MQ=255)
                                                                                                                                     ctcAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGCACGCTGGACGAt              >  1:670992/1‑150 (MQ=255)
                                                                                                                                        aaTGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGCACGCTGGACGATGcc           >  2:727205/1‑150 (MQ=255)
                                                                                                                                             ccTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGCACGCTGGACGATGCCACgcc      >  1:197273/1‑150 (MQ=255)
                                                                                                                                             ccTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGCACGCTGGACGATGCCACgcc      >  1:1191987/1‑150 (MQ=255)
                                                                                                                                                 gTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGCACGCTGGACGATGCCACgccggcc  <  1:1078485/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTGGAGCGCGCCGAACGAGGCCTACGAAGGCGCCTGCGCTCATTACCTCGACGGCCTGCTGCTGGACAGCGAGAACCAGCAACTGCGCAAGTCGCTGGCCGACGCGGCGCAGCTGCTGGCCTGCCCCGGCGCTCTCAATGGCCTGGTCCAGGCCCTGTTGCGCATGACCACACCCGGTGTGCCCGACCTGTACCAGGGTAACGAATACTGGGACTTCAGCCTGGTCGACCCGGACAACCGCCGCGCCGTGGACTACGCTTCGAGGCGTCGCACGCTGGACGATGCCACGCCGGCC  >  NC_002947/4572253‑4572547

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