Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,580,596 C→G 13.8% G959A (GGC→GCC)  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,5960CG13.8% 135.7 / 11.3 49G959A (GGC→GCC) treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base C (25/17);  new base G (3/4);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGCAGCAGCATAGTCGAAGGATCGCAGCACGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGA  >  NC_002947/4580453‑4580739
                                                                                                                                               |                                                                                                                                               
tcGCAGCAGCATAGTCGAAGGATCGCAGCACGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGt                                                                                                                                           >  1:549673/1‑150 (MQ=255)
tcGCAGCAGCATAGTCGAAGGATCGCAGCACGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGt                                                                                                                                           >  1:510442/1‑150 (MQ=255)
    agcagcaTAGTCGAAGGATCGCAGCACGCCGCTGACCTATTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGGCCAGGTGctc                                                                                                                                       <  1:684862/150‑3 (MQ=255)
       agcaTAGTCGAAGGATCGCAGGACGCCGCTGACATCGTTGTAAGGGCTGTGTTTGGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTc                                                                                                                                    <  2:403796/150‑1 (MQ=255)
         caTAGTCGAAGGATCGCAGCACGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGGCCAGGTGctcgtggt                                                                                                                                  <  2:456513/150‑8 (MQ=255)
         caTAGTCGAAGGATCGCAGCACGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGgcgc                                                                                                                                  >  1:359425/1‑150 (MQ=255)
          aTAGTCGAAGGATCGCAGCACGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGcc                                                                                                                                 <  2:359425/150‑1 (MQ=255)
               cGAAGGATCGCAGCACGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTa                                                                                                                            <  2:510442/150‑1 (MQ=255)
                        gcagcaCGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGGCCAGGTGct                                                                                                                                        <  2:90325/129‑2 (MQ=255)
                        gcagcaCGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGGCCAGGTGct                                                                                                                                        >  1:90325/1‑128 (MQ=255)
                                ccgcTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCa                                                                                                           <  1:229698/150‑1 (MQ=255)
                                    tGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCAGCTTCAAAATCTATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCgt                                                                                                                               <  1:254361/126‑1 (MQ=255)
                                    tGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCgt                                                                                                                               >  2:254361/1‑126 (MQ=255)
                                        aTCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGTCGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGc                                                                                                   <  1:649634/150‑1 (MQ=255)
                                        aTCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGAACGGATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAgg                                                                                                   >  1:484115/1‑149 (MQ=255)
                                          ccTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTc                                                                                                                                    >  2:20799/1‑115 (MQ=255)
                                          ccTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTc                                                                                                                                    <  1:20799/115‑1 (MQ=255)
                                            ttGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGc                                                                                                                 >  1:158486/1‑132 (MQ=255)
                                            ttGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGc                                                                                                                 <  2:158486/132‑1 (MQ=255)
                                                               cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgt                                                                            >  1:485710/1‑150 (MQ=255)
                                                                ccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtgtg                                                                           >  1:768881/1‑150 (MQ=255)
                                                                 cGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGc                                                                          >  1:529353/1‑150 (MQ=255)
                                                                         ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAg                                                                  >  1:661449/1‑150 (MQ=255)
                                                                         ttCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAg                                                                  >  1:459832/1‑150 (MQ=255)
                                                                                    gggcgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAg                                                                                                               >  1:41277/1‑94 (MQ=255)
                                                                                    gggcgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAg                                                                                                               <  2:41277/94‑1 (MQ=255)
                                                                                        gggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCt                                                                         >  2:252490/1‑128 (MQ=255)
                                                                                        gggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCt                                                                         <  1:252490/128‑1 (MQ=255)
                                                                                                cTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgt                                           >  2:344806/1‑150 (MQ=255)
                                                                                                 ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAGGgtggcct                                                                               <  2:338279/113‑5 (MQ=38)
                                                                                                 ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAGGgtggcct                                                                               >  1:338279/1‑109 (MQ=38)
                                                                                                           aTGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAg                                <  2:661449/150‑1 (MQ=255)
                                                                                                                     cATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTAt                      <  1:344806/150‑1 (MQ=255)
                                                                                                                      aTCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATc                     >  2:840917/1‑150 (MQ=255)
                                                                                                                           ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAt                                                                              <  1:645707/88‑1 (MQ=255)
                                                                                                                           ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAt                                                                              >  2:645707/1‑88 (MQ=255)
                                                                                                                           ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                 <  1:806397/133‑1 (MQ=255)
                                                                                                                           ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                 >  2:806397/1‑133 (MQ=255)
                                                                                                                           ccTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGt                >  2:842142/1‑150 (MQ=255)
                                                                                                                                accacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt           <  1:842142/150‑1 (MQ=255)
                                                                                                                                  cacGAGCACCTGGCCGAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTTGCTATCGCTGTCCAgggt         >  2:408831/1‑150 (MQ=255)
                                                                                                                                  cacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgtat         >  2:642011/1‑148 (MQ=255)
                                                                                                                                     ctgCACCTGGGCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGa                                                                                                     <  2:428303/53‑1 (MQ=38)
                                                                                                                                     ctgCACCTGGGCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGa                                                                                                     >  1:428303/3‑55 (MQ=255)
                                                                                                                                        cACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTcgcg   >  2:477655/1‑150 (MQ=255)
                                                                                                                                         aCCTGGGCCAGGTGCAGGTCGCCGTGTACACGCATTAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGa  >  2:333553/1‑150 (MQ=255)
                                                                                                                                         aCCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGa  >  1:580158/1‑150 (MQ=255)
                                                                                                                                           cTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTg                                                                        >  2:827122/1‑78 (MQ=255)
                                                                                                                                           cTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTg                                                                        <  1:827122/78‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TCGCAGCAGCATAGTCGAAGGATCGCAGCACGCCGCTGACATCCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGA  >  NC_002947/4580453‑4580739

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