Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,655,155 T→A 10.9% S131C (AGC→TGC)  epd ← D‑erythrose 4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,655,1550TA10.9% 148.4 / 8.5 55S131C (AGC→TGC) epdD‑erythrose 4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/24);  new base A (2/4);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

GATGCCAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTGGCACTGC  >  NC_002947/5655012‑5655301
                                                                                                                                               |                                                                                                                                                  
gATGCCAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCCGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACACCCGTGGCGTCAACGTCGGCCTCGCTGGCCa                                                                                                                                               >  3:6665/1‑149 (MQ=255)
      aaaCGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCAGCCCATggc                                                                                                                                           >  7:207093/1‑146 (MQ=255)
        aCGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCAAGGGCTGc                                                                                                                                       <  1:180674/149‑1 (MQ=255)
           ccTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATg                                                                                                                                             >  4:249132/1‑140 (MQ=255)
           ccTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATg                                                                                                                                             <  3:249132/140‑1 (MQ=255)
           ccTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCAGCCCATGGCCAGCgag                                                                                                                                    <  2:111123/149‑1 (MQ=255)
            cTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCgaga                                                                                                                                   <  5:172918/149‑1 (MQ=255)
            cTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCAGCCCATGGCCAGCgagg                                                                                                                                   <  2:300075/149‑2 (MQ=255)
             tGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGaa                                                                                                                                  >  7:283657/1‑149 (MQ=255)
               gTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCTCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACa                                                                                                                                <  6:240588/149‑1 (MQ=255)
               gTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCAGCCCATGGCCAGCgaggc                                                                                                                                  <  8:207093/147‑3 (MQ=255)
               gTCCAGCACAAGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACa                                                                                                                                <  8:288560/149‑1 (MQ=255)
                tCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAg                                                                                                                               <  5:304249/149‑1 (MQ=255)
                 ccAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGc                                                                                                                              <  4:6665/149‑1 (MQ=255)
                  cAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCa                                                                                                                             <  7:215150/149‑1 (MQ=255)
                      acacGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGTTGGCCATGGGCTGCGAGAACAGCACCCt                                                                                                                         >  3:368252/1‑149 (MQ=255)
                          gcagcaGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCGTGGGCTGCGAGAACAGCACCCTTGgc                                                                                                                     <  7:234734/149‑1 (MQ=255)
                          gcagcaGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGgc                                                                                                                     <  3:313681/149‑1 (MQ=255)
                          gcagcaGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGgc                                                                                                                     >  5:83854/1‑149 (MQ=255)
                          gcagcaGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGgc                                                                                                                     >  4:313681/1‑149 (MQ=255)
                           cagcagCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGgcg                                                                                                                    <  2:164109/149‑1 (MQ=255)
                           cagcagCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGgcg                                                                                                                    >  3:310455/1‑149 (MQ=255)
                            agcagcGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGgcgc                                                                                                                   <  4:310455/149‑1 (MQ=255)
                                   gCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTgg                                                                                                                      <  5:51427/139‑1 (MQ=255)
                                   gCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTgg                                                                                                                      >  6:51427/1‑139 (MQ=255)
                                    cACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCtt                                                                                                           >  5:24531/1‑149 (MQ=255)
                                       gccgcAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTc                                                                                                          >  4:193026/1‑147 (MQ=255)
                                         cgcAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAg                                                                                                       >  5:310516/1‑148 (MQ=255)
                                                           aGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCa                                                                                    <  6:83854/149‑1 (MQ=255)
                                                             gCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATTGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACg                                                                                  >  6:213930/1‑149 (MQ=255)
                                                                ttGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGgt                                                                               >  2:59287/1‑149 (MQ=255)
                                                                ttGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGgt                                                                               >  8:230023/1‑149 (MQ=255)
                                                                      acCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAATGGCTGGCCATCGGCACGGGTGTtata                                                                         <  2:327315/149‑1 (MQ=255)
                                                                      acCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTtata                                                                         <  3:193026/149‑1 (MQ=255)
                                                                      acCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTtata                                                                         <  6:310516/149‑1 (MQ=255)
                                                                       ccAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGTTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAg                                                                        <  4:368252/149‑1 (MQ=255)
                                                                        cAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCCTCGGCACGGGTGTTATAgg                                                                       >  1:273765/1‑149 (MQ=255)
                                                                        cAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAgg                                                                       >  5:49882/1‑149 (MQ=255)
                                                                                     ccGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACAGGTGTTATAGGCACCGGAACACTc                                                          <  8:183849/149‑1 (MQ=255)
                                                                                             gcaggcCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCa                                                        <  6:366342/143‑1 (MQ=255)
                                                                                             gcaggcCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCa                                                        >  5:366342/1‑143 (MQ=255)
                                                                                                      gTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGAt                                           <  8:283657/147‑1 (MQ=255)
                                                                                                          aTGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAg                                     >  4:169996/1‑149 (MQ=255)
                                                                                                              cGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAgcg                                   >  5:343934/1‑147 (MQ=255)
                                                                                                                tAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGc                               <  1:59287/149‑1 (MQ=255)
                                                                                                                          tGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATg                                          >  5:93073/1‑128 (MQ=255)
                                                                                                                          tGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATg                                          <  6:93073/128‑1 (MQ=255)
                                                                                                                                    gTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCggggg           >  4:79100/1‑149 (MQ=255)
                                                                                                                                        ttcTCGCAGCCCATGGCCAGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTgg        <  6:171245/146‑1 (MQ=255)
                                                                                                                                        ttcTCGCAGCCCATGGCCAGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTgg        >  5:171245/3‑148 (MQ=255)
                                                                                                                                         ccTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATg                                          <  5:27611/113‑1 (MQ=255)
                                                                                                                                         ccTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATg                                          >  6:27611/1‑113 (MQ=255)
                                                                                                                                         ccTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTGGCa      >  6:56398/1‑149 (MQ=255)
                                                                                                                                           tCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTGGCACt    >  8:86081/1‑149 (MQ=255)
                                                                                                                                             gCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTGGCACTgc  >  6:80325/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GATGCCAAACGCCTGGTCCAGCACACGCAGCAGCGGCACGCCGCAGTTGGTGGTGCACGAGGCGTTGGACACCAGGCGCTCGCTGCCGGTAAGGCAGGCCTGGTTGATGCCGTAGACAACCGTGGCGTCAACGTCGGCCTCGCTGGCCATGGGCTGCGAGAACAGCACCCTTGGCGCACCGGCTTCGAGGAAACGCTGGCCATCGGCACGGGTGTTATAGGCACCGGAACACTCCAGCACCAGGTCGATGCCCAGCGCAGCCCAGTCGATGCCTTCGGGGGTGGCACTGC  >  NC_002947/5655012‑5655301

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