Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,346,335 T→G 7.8% L578R (CTG→CGG) ‡ acsA‑II → acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,346,3350TG 7.8% 151.4 / 5.0 51L578R (CTG→CGG) ‡acsA‑IIacetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base T (22/25);  new base G (2/2);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

GAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTGGGCGATAC  >  NC_002947/5346187‑5346476
                                                                                                                                                    |                                                                                                                                             
gAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCt                                                                                                                                               >  5:31579/1‑149 (MQ=255)
             tgctgcATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGACCAGGGCCAGCagctcg                                                                                                                                  <  8:3332/149‑6 (MQ=255)
                tgcATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATc                                                                                                                               >  3:39835/1‑149 (MQ=255)
                  cATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATc                                                                                                                               <  1:140020/147‑1 (MQ=255)
                  cATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATc                                                                                                                               <  5:121088/147‑1 (MQ=255)
                    tGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCgg                                                                                                                             >  2:99896/1‑147 (MQ=255)
                          cGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCg                                                                                                                     >  7:118912/1‑149 (MQ=255)
                                 gAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGcc                                                                                                              >  3:183010/1‑149 (MQ=255)
                                 gAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGcc                                                                                                              >  1:210210/1‑149 (MQ=255)
                                      cgcgGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGc                                                                                                         <  4:153628/149‑1 (MQ=255)
                                              tCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGACCAGGGCCAGCag                                                                                                                                      >  6:225571/1‑111 (MQ=255)
                                              tCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGACCAGGGCCAGCag                                                                                                                                      <  5:225571/112‑2 (MQ=255)
                                                ggCTACCCCCACGTCCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGt                                                                                                <  3:147031/148‑1 (MQ=255)
                                                ggCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTg                                                                                               <  1:99896/149‑1 (MQ=255)
                                                ggCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTg                                                                                               <  8:176960/149‑1 (MQ=255)
                                                 gCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTgg                                                                                              <  3:212847/149‑1 (MQ=255)
                                                 gCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTgg                                                                                              >  4:212847/1‑149 (MQ=255)
                                                     cccccACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGgcgc                                                                                          <  6:31579/149‑1 (MQ=255)
                                                      ccccACGACCTCAACGGCCAAGTCGTCTACGCCTTCGTCACGACAAAGAATCGCGTGACACAAGATGACACCCTCAAGGCCGAGATGCTGGCCCTGGTCAGCAAGGAAATCGCCAGCTTCGCAAACCCTGAGATGATTCAGACGgcgc                                                                                          >  4:182401/1‑148 (MQ=255)
                                                               cTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTAcc                                                                                >  5:183567/1‑149 (MQ=255)
                                                                 cAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTAccc                                                                               >  1:121706/1‑148 (MQ=255)
                                                                   aGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAg                                                                                                          <  4:88292/119‑1 (MQ=255)
                                                                   aGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAg                                                                                                          >  3:88292/1‑119 (MQ=255)
                                                                        cAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGAcgcg                                                                       >  8:230488/1‑149 (MQ=255)
                                                                         aaGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACgcgc                                                                      <  2:121706/149‑1 (MQ=255)
                                                                            gCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCa                                                                   <  7:230488/149‑1 (MQ=255)
                                                                                 tACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGc                                                                >  7:144910/1‑147 (MQ=255)
                                                                                      cTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCa                                                         >  2:139473/1‑149 (MQ=255)
                                                                                        tCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATg                                                       >  3:239646/1‑149 (MQ=255)
                                                                                           tCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATgcgg                                                    <  7:123691/149‑1 (MQ=255)
                                                                                             aCGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATgcggcg                                                  >  7:77275/1‑149 (MQ=255)
                                                                                                 ccATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGtata                                              <  3:40253/149‑1 (MQ=255)
                                                                                                   aTGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACt                                            <  6:183567/149‑1 (MQ=255)
                                                                                                    tGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTg                                           <  4:39835/149‑1 (MQ=255)
                                                                                                      aaTGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTgcg                                         <  2:210210/149‑1 (MQ=255)
                                                                                                      aaTGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTgc                                          <  4:239646/148‑1 (MQ=255)
                                                                                                        tGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGa                                                            >  8:206175/1‑128 (MQ=255)
                                                                                                        tGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGa                                                            <  7:206175/128‑1 (MQ=255)
                                                                                                               acacCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGc                                <  8:144910/149‑1 (MQ=255)
                                                                                                                cacCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGcc                               <  3:177499/149‑1 (MQ=255)
                                                                                                                 acCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCt                              >  8:82883/1‑149 (MQ=255)
                                                                                                                   cAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTg                             >  8:84266/1‑148 (MQ=255)
                                                                                                                     gATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTgc                                          <  1:9642/133‑1 (MQ=255)
                                                                                                                     gATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTgc                                          >  2:9642/1‑133 (MQ=255)
                                                                                                                                cAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGaa                >  5:179062/1‑148 (MQ=255)
                                                                                                                                 aaGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGaaaa              <  7:84266/149‑1 (MQ=255)
                                                                                                                                 aaGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGaaaa              >  3:74047/1‑149 (MQ=255)
                                                                                                                                      cGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTgg         <  2:21593/149‑1 (MQ=255)
                                                                                                                                       tccttgctgACCAGGGCCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTggg        >  7:212714/5‑149 (MQ=255)
                                                                                                                                       gAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTggg        <  8:77275/149‑1 (MQ=255)
                                                                                                                                        aGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTGGGc       <  1:8167/149‑1 (MQ=255)
                                                                                                                                             ctgGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTGGGCGATAc  >  5:193454/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GAAAGTGCACTGGTGCTGCATGACAGCGTGGCCGAAGCCGCGGTGGTCGGCTACCCCCACGACCTCAAGGGCCAAGGCGTCTACGCCTTCGTCACGACCATGAATGGCGTGACACCAGATGACACGCTCAAGGCCGAGCTGCTGGCCCTGGTCAGCAAGGAAATCGGCAGCTTCGCCAAGCCTGAGCTGATTCAGTGGGCGCCAGCGCTACCCAAGACGCGCTCAGGCAAGATCATGCGGCGTATACTGCGCAAAATCGCCTGCAATGAGCTGGAAAACCTGGGCGATAC  >  NC_002947/5346187‑5346476

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