Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,570,431 T→C 11.1% L84P (CTG→CCG)  treY → maltooligosyl trehalose synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,570,4310TC11.1% 89.6 / 2.4 36L84P (CTG→CCG) treYmaltooligosyl trehalose synthase
Reads supporting (aligned to +/- strand):  ref base T (16/16);  new base C (2/2);  total (18/18)
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.95e-01

ACCTGTATGCCTCGCCCATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTATGCAGAGTTCTTCGACATCCAGTGGCACTCCAGCGAC  >  NC_002947/4570284‑4570561
                                                                                                                                                   |                                                                                                                                   
aCCTGTATGCCTCGCCCATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTAGCTTCGTCAACCCGGGTCTTGGCGGCGAGGCCGCGCTGGAGCTCCTGGTTGCCTCCCTGCGCTAGCATGGCATGGGGCTGa                                                                                                                                   >  2:360478/1‑150 (MQ=255)
aCCTGTATGCCTCGCCCATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCGAGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGa                                                                                                                                   >  2:501211/1‑150 (MQ=255)
          cTCGCCCATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGacac                                                                                                                        >  1:33548/1‑150 (MQ=255)
               ccATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTc                                                                                                                             <  1:258458/140‑1 (MQ=255)
               ccATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTc                                                                                                                             <  1:258486/140‑1 (MQ=255)
               ccATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTc                                                                                                                             >  2:258486/1‑140 (MQ=255)
               ccATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTc                                                                                                                             >  2:258458/1‑140 (MQ=255)
                        aGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACcaca                                                                                                          <  2:33548/150‑1 (MQ=255)
                          gCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAAc                                                                                                              >  2:8209/1‑144 (MQ=255)
                          gCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAAc                                                                                                              <  1:8209/144‑1 (MQ=255)
                                        gcgcCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGc                                                                                          <  2:170319/150‑1 (MQ=255)
                                          gcCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCg                                                                                        >  1:25146/1‑150 (MQ=255)
                                             aTGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACa                                                                                     >  2:382174/1‑150 (MQ=255)
                                                      acgTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCtggt                                                                            <  1:502465/150‑1 (MQ=255)
                                                           gTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACg                                                                                                                       <  2:352381/102‑1 (MQ=255)
                                                           gTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACg                                                                                                                       >  1:352381/1‑102 (MQ=255)
                                                                tCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGc                                                                                          >  2:337697/1‑126 (MQ=255)
                                                                tCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGc                                                                                          <  1:337697/126‑1 (MQ=255)
                                                                      cTGCGTCAACCCGGAGCTGGGCGGCGAGCCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCANATGGCCNTAGGCGGTGCCGNCAACCCCTGGTGGCAAAGCCTGCTGGc                                                            <  1:458705/150‑1 (MQ=255)
                                                                                ccGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACgg                                                  >  1:324338/1‑150 (MQ=255)
                                                                                ccGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCCGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACgg                                                  >  1:246707/1‑150 (MQ=255)
                                                                                    aGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCa                                              <  1:200656/150‑1 (MQ=255)
                                                                                       tGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTc                                           >  1:256415/1‑150 (MQ=255)
                                                                                            ggcgCGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGCACTCACCACACGCACATAATCCAACCACATGGCCCTAGGCGGTGCCGACACCCCCTGGTGGGAAAGCCGGCTGGCCTGGGGACGGCGCAGTCCGTAt                                      <  1:369392/150‑1 (MQ=255)
                                                                                            ggcgAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGTCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTAt                                      <  1:360478/150‑1 (MQ=255)
                                                                                                           gAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCCGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTa                                       <  2:361518/134‑1 (MQ=255)
                                                                                                           gAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCCGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTa                                       >  1:361518/1‑134 (MQ=255)
                                                                                                              cgcCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCCGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTATGCAGAGTTCTTCGACATc                    <  2:246707/150‑1 (MQ=255)
                                                                                                                            ccTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTATGCAGAGttcttc                          >  2:161475/1‑130 (MQ=255)
                                                                                                                            ccTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTATGCAGAGttcttc                          <  1:161475/130‑1 (MQ=255)
                                                                                                                             cTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGgcg                                                <  2:230027/107‑1 (MQ=255)
                                                                                                                             cTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGgcg                                                >  1:230027/1‑107 (MQ=255)
                                                                                                                                          gcatggGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCa                                              >  1:367334/1‑96 (MQ=255)
                                                                                                                                          gcatggGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCa                                              <  2:367334/96‑1 (MQ=255)
                                                                                                                                                ggCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTATGCAGAGTTCTTCGACATCCAGTGGCACTCCAGCGAc  >  2:205722/1‑134 (MQ=255)
                                                                                                                                                ggCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTATGCAGAGTTCTTCGACATCCAGTGGCACTCCAGCGAc  <  1:205722/134‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                   
ACCTGTATGCCTCGCCCATCCTCAAGGCCCGCGCGGGCTCGCGCCATGGCTACGACGTGGTCGATCCTACCTGCGTCAACCCGGAGCTGGGCGGCGAGGCCGCGCTGGAGCGCCTGGTTGCCGCCCTGCGCCAGCATGGCATGGGGCTGATCCTCGACACGGTGTCCAACCACATGGCCGTAGGCGGTGCCGACAACCCCTGGTGGCAAAGCCTGCTGGCCTGGGGACGGCGCAGTCCGTATGCAGAGTTCTTCGACATCCAGTGGCACTCCAGCGAC  >  NC_002947/4570284‑4570561

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