Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,144,248 T→A 7.4% I93I (ATT→ATA)  PP_3648 → carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,144,2480TA 7.4% 155.3 / 3.3 54I93I (ATT→ATA) PP_3648carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (27/23);  new base A (2/2);  total (29/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 = 9.54e-01

CTGGGGTCACGTCTGGGGCCGCGAAGGCTTGTCGCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGTTCGAACCACTGCATGGAAAGCCTCG  >  NC_002947/4144101‑4144395
                                                                                                                                                   |                                                                                                                                                   
cTGGGGTCACGTCTGGGGCCGCGAAGGCTTGTCGCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTc                                                                                                                                                    >  3:361025/1‑149 (MQ=255)
               gggCCGCGAAGGCTTGTCGCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTcc                                                                                                                                     <  4:361025/149‑1 (MQ=255)
               gggCCGCGAAGGCTTGTCGCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTc                                                                                                                                      >  4:309501/1‑148 (MQ=255)
                    gcgAAGGCTTGTCGCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTc                                                                                                                                <  6:19425/149‑1 (MQ=255)
                              gTCGCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCAt                                                                                                                           <  5:284928/144‑1 (MQ=255)
                              gTCGCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCAt                                                                                                                           >  6:284928/1‑144 (MQ=255)
                                 gCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTcgc                                                                                                                                                  >  4:263826/1‑118 (MQ=255)
                                 gCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTcgc                                                                                                                                                  <  3:263826/118‑1 (MQ=255)
                                    gAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTgcggcg                                                                                                                >  7:64581/1‑149 (MQ=255)
                                       agagCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATACGCGaaa                                                                                                                                              >  3:141366/1‑116 (MQ=255)
                                       agagCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATACGCGaaa                                                                                                                                              <  4:141366/116‑1 (MQ=255)
                                       agagCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATACGCGaaa                                                                                                                                              <  4:135704/116‑1 (MQ=255)
                                       agagCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATACGCGaaa                                                                                                                                              >  3:135704/1‑116 (MQ=255)
                                                caTGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTcg                                                                                                                                                   <  7:20825/102‑1 (MQ=255)
                                                caTGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTcg                                                                                                                                                   >  8:20825/1‑102 (MQ=255)
                                                  tGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTg                                                                                                  <  5:60074/149‑1 (MQ=255)
                                                  tGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTg                                                                                                  >  1:251801/1‑149 (MQ=255)
                                                     tAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGa                                                                                                >  5:152937/1‑148 (MQ=255)
                                                         cTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCcgc                                                                                                         <  8:163593/135‑1 (MQ=255)
                                                         cTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCcgc                                                                                                         >  7:163593/1‑135 (MQ=255)
                                                              ccATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGc                                                                                      <  6:152937/149‑1 (MQ=255)
                                                               cATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGcc                                                                                     >  1:89946/1‑149 (MQ=255)
                                                                    tttCCGCGCTCAACAGGCCCCATGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGc                                                                                 <  5:44349/148‑1 (MQ=255)
                                                                       ccGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGt                                                                             >  2:36752/1‑149 (MQ=255)
                                                                        cgcgCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCTTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGt                                                                             >  8:115935/1‑148 (MQ=255)
                                                                        cgcgCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTg                                                                            <  2:251801/149‑1 (MQ=255)
                                                                             tCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTa                                                                                                                        >  6:161420/1‑100 (MQ=255)
                                                                               aaCAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGt                                                                                                   <  1:334003/119‑1 (MQ=255)
                                                                               aaCAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGt                                                                                                   >  2:334003/1‑119 (MQ=255)
                                                                               aaCAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGaa                                                                          <  7:230708/144‑1 (MQ=255)
                                                                               aaCAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGaa                                                                          >  8:230708/1‑144 (MQ=255)
                                                                                aCAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGttt                                                                    <  2:89946/149‑1 (MQ=255)
                                                                                 cAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTg                                                                   >  1:322510/1‑149 (MQ=255)
                                                                                  aGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGc                                                                  <  8:64581/149‑1 (MQ=255)
                                                                                  aGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGc                                                                  >  1:187267/1‑149 (MQ=255)
                                                                                    gCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTgcggc                                                                                                                 >  8:164966/1‑100 (MQ=255)
                                                                                    gCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGcc                                                                 >  8:17284/1‑148 (MQ=255)
                                                                                        cACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGc                                                            >  7:144521/1‑149 (MQ=255)
                                                                                         aCGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGcc                                                           >  1:44242/1‑149 (MQ=255)
                                                                                             gCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATg                                                       <  2:187267/149‑1 (MQ=255)
                                                                                               tCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCg                                                     >  7:258985/1‑149 (MQ=255)
                                                                                               tCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCg                                                     >  7:348194/1‑149 (MQ=255)
                                                                                                 aaGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCg                                                     <  4:107151/147‑1 (MQ=255)
                                                                                                           ttCGCGGTGCCTTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCAATGCGGAGCCAACCCGt                                         <  7:115935/149‑1 (MQ=255)
                                                                                                             cgcgGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGtt                                        <  7:17284/148‑1 (MQ=255)
                                                                                                                      cTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGc                              >  8:224737/1‑149 (MQ=255)
                                                                                                                          gcAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGtt                          >  2:220690/1‑149 (MQ=255)
                                                                                                                           cAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGTTc                         >  8:228841/1‑149 (MQ=255)
                                                                                                                              tGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGTTCGaa                      <  2:44242/149‑1 (MQ=255)
                                                                                                                                ggCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGTTCGAAc                     <  3:309501/148‑1 (MQ=255)
                                                                                                                                  cTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGTTCGAACCAc                  <  8:258985/149‑1 (MQ=255)
                                                                                                                                        cGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGTTCTAACCACTGCATg            <  1:36752/149‑1 (MQ=255)
                                                                                                                                             aCAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGtt                          <  1:127891/130‑1 (MQ=255)
                                                                                                                                             aCAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGtt                          >  2:127891/1‑130 (MQ=255)
                                                                                                                                                  ttCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGTTCGAACCACTGCATGGAAAGCCTCg  <  8:144521/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGGGGTCACGTCTGGGGCCGCGAAGGCTTGTCGCTGAAAGAGCGCAGCATGATAAACTTGGCCATGATTTCCGCGCTCAACAGGCCCCACGAGCTCAAGTTGCACATTCGCGGTGCCCTGCGCAATGGGCTGAGCCGTGAACAGATTCGCGAAATCCTGCTCCAGGTCGGCATTTATTGCGGCGTACCCGCGGCGGTGGACAGTTTCCGCCTGGCCCGTGAAGCGTTTGCCGAAGCCGATGCGGAGCCAACCCGTTAACAACGGGCTGTTCGAACCACTGCATGGAAAGCCTCG  >  NC_002947/4144101‑4144395

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