Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,504,734 T→A 10.4% L174Q (CTG→CAG)  yliI → aldose sugar dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,504,7340TA10.4% 123.9 / 5.7 48L174Q (CTG→CAG) yliIaldose sugar dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (20/23);  new base A (2/3);  total (22/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 = 8.40e-01

TCGGCCGCGGGCAGCTGTCCGAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTTCGTCGGCAAGGACAACGTGCGCCCGGAGATCTGGTCGTACGGCC  >  NC_002947/2504593‑2504880
                                                                                                                                             |                                                                                                                                                   
tCGGCCGCGGGCAGCTGTCCGAAGACCGCGCCCTGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCga‑g                                                                                                                                             <  5:220115/149‑1 (MQ=255)
tCGGCCGCGGGCAGCTGTCCGAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGc                                                                                                                                                 >  6:203921/1‑146 (MQ=255)
     cgcgGGCAGCTGTCCGAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTC‑TCGCCCAGGGCGATGaacag                                                                                                                                        <  6:286083/149‑2 (MQ=255)
          gCAGCTGTCCGAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAAcaaccaac                                                                                                                                   <  5:203921/149‑1 (MQ=255)
                 tCCGAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCTGCTCGCGGCTGGTGTTCGACCTCGAAGGCTACCTGTTCAACGCCCTGGGCGA‑GAACAACCAACGTCCCa                                                                                                                             >  8:276161/1‑148 (MQ=255)
                 tCCGAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTC‑TCGCCCAGGGCGATGaacag                                                                                                                                        <  5:275805/137‑2 (MQ=255)
                 tCCGAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTC‑TCGCCCAGGGCGATGaacag                                                                                                                                        >  6:275805/1‑136 (MQ=255)
                    gAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGc                                                                                                                         >  2:109486/1‑149 (MQ=255)
                          cgcgCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGa                                                                                                                   >  3:97442/1‑149 (MQ=255)
                            cgcCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGAcc                                                                                                                 >  1:139030/1‑149 (MQ=255)
                                cGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGTCCAGGTCCTggg                                                                                                             <  5:12615/149‑2 (MQ=255)
                                cGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGa                                                                                                             <  4:97442/149‑1 (MQ=255)
                                   cTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACaa                                                                                                          <  8:45097/149‑1 (MQ=255)
                                      gagaATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGtt                                                                                                       >  7:85446/1‑149 (MQ=255)
                                         aaTTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTC‑TCGCCCAGGGCGATGaacag                                                                                                                                        >  6:236033/1‑112 (MQ=255)
                                         aaTTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTC‑TCGCCCAGGGCGATGaacag                                                                                                                                        <  5:236033/113‑2 (MQ=255)
                                          aTTTCACCGTGATCTTCCGCCTGCAGCACAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaa                                                                                                   <  7:276161/149‑1 (MQ=255)
                                          aTTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaa                                                                                                   <  8:85446/149‑1 (MQ=255)
                                          aTTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaa                                                                                                   <  2:139030/149‑1 (MQ=255)
                                             tCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaagg                                                                                                 <  8:244048/148‑1 (MQ=255)
                                               aCCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaaggc                                                                                                >  6:127265/1‑147 (MQ=255)
                                                ccGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaaggcaag                                                                                             <  4:157298/149‑1 (MQ=255)
                                                           cGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTa                                                                                   <  8:286255/148‑1 (MQ=255)
                                                            gCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GCACAACCAACGTACCCCTGCCCCGGACCTGGACAAGTTGCAAGGCAAGGTCGGGCGTAtt                                                                                 >  1:79601/1‑149 (MQ=255)
                                                              cagcagCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaaggca                                                                                               <  2:153439/133‑1 (MQ=255)
                                                              cagcagCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaaggca                                                                                               >  1:153439/1‑133 (MQ=255)
                                                                   gcCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCg                                                                           >  5:138220/1‑148 (MQ=255)
                                                                     ccAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTg                                                                                       >  5:121457/1‑134 (MQ=255)
                                                                     ccAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTg                                                                                       <  6:121457/134‑1 (MQ=255)
                                                                      cAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGc                                                                                      >  8:152423/1‑134 (MQ=255)
                                                                      cAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGc                                                                                      <  7:152423/134‑1 (MQ=255)
                                                                           tGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGaa                                                                  >  3:626/1‑149 (MQ=255)
                                                                             tCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaaggcaag                                                                                             >  5:258079/1‑120 (MQ=255)
                                                                             tCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTgcaaggcaag                                                                                             <  6:258079/120‑1 (MQ=255)
                                                                             tCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGt                                                                <  4:626/149‑1 (MQ=255)
                                                                                   ggTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCa                                                           <  5:127265/148‑1 (MQ=255)
                                                                                         cACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGt                                                                                        >  3:210854/1‑113 (MQ=255)
                                                                                         cACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGt                                                                                        <  4:210854/113‑1 (MQ=255)
                                                                                         cACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACaa                                                    >  2:218900/1‑149 (MQ=255)
                                                                                                    gcgGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTtcgtc                                          <  1:109486/148‑1 (MQ=255)
                                                                                                              ttCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTTCGTCGGCAAGGACaa                               >  5:38320/1‑149 (MQ=255)
                                                                                                                  aCCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGcc                                                             <  6:45312/115‑1 (MQ=255)
                                                                                                                  aCCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGcc                                                             >  5:45312/1‑115 (MQ=255)
                                                                                                                               ccTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTTCGTCGGCAAGGACAACGTGCGCCCGGAGATc               <  6:38320/148‑1 (MQ=255)
                                                                                                                               ccTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTTCGTCGGCAAGGACAACGTGCGCCCGGAGATCt              <  1:218900/149‑1 (MQ=255)
                                                                                                                                           ccTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTTCGTCGGCAAGGACAACGTgc                          <  6:322267/125‑1 (MQ=255)
                                                                                                                                           ccTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTTCGTCGGCAAGGACAACGTgc                          >  5:322267/1‑125 (MQ=255)
                                                                                                                                           ccTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTTCGTCGGCAAGGACAACGTGCGCCCGGAGATCTGGTCGTACGGcc  >  4:35902/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TCGGCCGCGGGCAGCTGTCCGAAGACCGCGCCCGGCTGGAGAATTTCACCGTGATCTTCCGCCAGCAGCCCAAGCTGTCGGTGGGTAACCACTTCGGCTCGCGGCTGGTGTTCGACCGCGACGGCTACCTGTTCATCGCCCTGGGCGA‑GAACAACCAACGTCCCACTGCCCAGGACCTGGACAAGTTGCAAGGCAAGGTCGTGCGTATTCTGCCGGACGGCGAAGTGCCCAAAGACAACCCCTTCGTCGGCAAGGACAACGTGCGCCCGGAGATCTGGTCGTACGGCC  >  NC_002947/2504593‑2504880

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