Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,839 C→G 13.6% intergenic (+76/+677) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8390CG13.6% 118.1 / 10.9 51intergenic (+76/+677)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (22/21);  new base G (3/4);  total (26/25)
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.09e-01

TCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGC  >  NC_002947/4740693‑4740986
                                                                                                                                                  |                                                                                                                                                     
tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                                                                    <  2:349122/149‑1 (MQ=255)
  ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                                                                   <  1:268181/148‑1 (MQ=255)
     cGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcg                                                                                                                                               >  3:90952/1‑149 (MQ=255)
         cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                                                                           <  3:180302/149‑1 (MQ=255)
         cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                                                                            <  1:91866/148‑1 (MQ=255)
           cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  5:155687/149‑1 (MQ=255)
            gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                                                                        <  3:18973/149‑1 (MQ=38)
               aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGcc                                                                                                                                     <  4:90952/149‑1 (MQ=255)
                       acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                                                        <  6:202456/138‑1 (MQ=255)
                       acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                                                                        >  5:202456/1‑138 (MQ=255)
                       acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCACGCGAAGAGGCCGGCACAgc                                                                                                                             >  8:80706/1‑148 (MQ=255)
                                 aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         >  8:153787/1‑127 (MQ=255)
                                 aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                         <  7:153787/127‑1 (MQ=255)
                                      cgcgACAGCGTTGCCCTTAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGACGGCTGCCCTCGGGCAGCCTTTTTAATTTGACCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                                                                             <  3:226276/149‑1 (MQ=255)
                                                      tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                         <  3:163833/106‑1 (MQ=25)
                                                      tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                         >  4:163833/1‑106 (MQ=25)
                                                        aaaGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTATACGCCCGCt                                                                                            >  8:301340/1‑148 (MQ=255)
                                                           caCCGAAAGTAAGCGGCAGTAGCCTAACGAAAGAAGGCTGCCCTCGGGCAGCTTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACg                                                                                        <  8:116639/148‑1 (MQ=255)
                                                                   gTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGACCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCtt                                                                                >  2:54581/1‑149 (MQ=255)
                                                                    tAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCttt                                                                               >  7:198557/1‑149 (MQ=255)
                                                                        cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGCGGCCGGCACAGACACCCCAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                           >  5:352618/1‑149 (MQ=255)
                                                                        cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                                                           <  8:198557/149‑1 (MQ=255)
                                                                           cAGTAGCCTACCGCATCAAGGCGGCCCTCGGGCAGCTTTTTTTATTTGTCGTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTa                                                                        <  4:350610/149‑1 (MQ=255)
                                                                                         aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGCCCCCACAAAACCTCAGGTTTTCTCCGCCCGCTCCCGCAGCCCTTTCAGGGGATTGCACGGCGCATc                                                         >  7:56001/1‑149 (MQ=255)
                                                                                           gaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTACGCCCGCACCCGCAGCCCTTTCAGGGTATTGAACGGCGCCTcca                                                       >  4:255224/1‑149 (MQ=255)
                                                                                                 tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCAAGACACCACAAAACCTCAGTTTTTCTCAGCCCGCTAACGCAGCCCTTTAAGGGTATTGAACGGCGCATCCACCACg                                                  >  7:249603/1‑148 (MQ=255)
                                                                                                 tGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCACTTTTTATCCCCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg                                                  >  1:198702/1‑148 (MQ=255)
                                                                                                   cccTCGGTCAGGCTTTTTTATTTGCCCTGTGCCGGCCTCTTCGCGGCCTTGCCAGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                               <  8:241129/149‑1 (MQ=255)
                                                                                                     cTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCCGCCATGCCCGCGAAGAGGCCGGCACAAAAACCACAAAACCACACATTTTCTCCGCCCGCACCCGCAGCCCATTAAGGGTATTGAACCGCGCATCCACCACGAACta                                             >  2:285735/1‑148 (MQ=255)
                                                                                                      tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttg                                            >  8:149691/1‑149 (MQ=255)
                                                                                                      tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCTTTCCCGCGAAGAGGCCGGCACAGAAACCACAAAACAACAGTTTTTATACCCCCCCTCACGCACCCCGTTCAGTGTATTTAACGGCGCATc                                                         >  8:221383/1‑136 (MQ=255)
                                                                                                      tCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTATGCGCGGCCTTGCCCGCTAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                                                         <  7:221383/136‑1 (MQ=255)
                                                                                                          gcctgcattttttATTTGCCCTGTGCGGGCCTATTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg                                         <  7:301340/142‑1 (MQ=255)
                                                                                                          gCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTCGCGGACATGCCCCCGAAGCGGACGGCACACCCACCACAAAACCACACTTTTTCTCCCCCCGCTAACGCAGCCCTTTCAGGGTAATTAACGGCGCATCCACCACGAACTTGATgg                                        >  2:100556/1‑149 (MQ=255)
                                                                                                               ccttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                                                          <  8:267592/110‑1 (MQ=255)
                                                                                                               ccttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAggg                                                                          >  7:267592/3‑112 (MQ=255)
                                                                                                                 ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                                                                   <  6:177011/85‑1 (MQ=38)
                                                                                                                 ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTcc                                                                                                   >  5:177011/1‑85 (MQ=38)
                                                                                                                 ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacca                                  >  5:186011/1‑149 (MQ=255)
                                                                                                                      tATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGc                              <  6:186011/148‑1 (MQ=255)
                                                                                                                              ctgTGCCGGCCTCTTAGCGGCCTTTCCCCCGAAGAGGCCGGCCCAGACACCACAAAAACACAGTTTTTATACGCCCGCTCACGCAGCCCGTTAAGGGTATTTAAAGGCGCATCCACCACGAACTTGTTGGCCACCATAGCCGGCCCGc                      >  4:350583/3‑148 (MQ=255)
                                                                                                                                   ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCCGTTTGTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccac                >  5:254081/1‑149 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCAGCCCGCTCACGCAGCCCTTTCAGGGTATTGACCGGCGCATCCACCACGAACTTGTTTGCCACCATCGCCGGCACGCTACCACCTg             <  7:179529/149‑1 (MQ=39)
                                                                                                                                         tctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt          >  7:8529/1‑149 (MQ=255)
                                                                                                                                         tctTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt          >  6:271917/1‑149 (MQ=255)
                                                                                                                                          ctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTgg            >  7:183251/1‑146 (MQ=255)
                                                                                                                                          ctTCGCGGGCAAGGCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTgg            <  8:183251/146‑1 (MQ=255)
                                                                                                                                           ttCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg        >  4:196786/1‑149 (MQ=255)
                                                                                                                                              gcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtg     <  1:54581/149‑1 (MQ=255)
                                                                                                                                               cgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtg     <  1:285735/149‑1 (MQ=255)
                                                                                                                                                 gCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGc  <  3:196786/149‑1 (MQ=255)
                                                                                                                                                  ccTTGCCCGCGAAGAGGCCGGCACAGACACGCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGtgtg   >  1:31862/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGTGC  >  NC_002947/4740693‑4740986

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: