Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 629,546 A→C T83P (ACA→CCA)  entH → enterobactin synthesis proofreading thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913629,5460AC98.1% 193.2 / ‑4.0 53T83P (ACA→CCA) entHenterobactin synthesis proofreading thioesterase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (33/19);  total (33/20)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTG  >  NC_000913/629402‑629692
                                                                                                                                                |                                                                                                                                                  
gtgCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAAc                                                                                                                                              >  1:8466/1‑151 (MQ=255)
        aGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCa                                                                                                                                      <  2:729704/151‑1 (MQ=255)
         gCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCa                                                                                                                                      >  1:342751/1‑150 (MQ=255)
           cGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAaca                                                                                                                                   <  2:675151/151‑1 (MQ=255)
           cGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAac                                                                                                                                    >  1:591381/1‑150 (MQ=255)
            taaaTGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAacac                                                                                                                                  >  2:276018/2‑151 (MQ=255)
            gAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAacac                                                                                                                                  >  2:789628/1‑151 (MQ=255)
            gAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAacac                                                                                                                                  >  1:573848/1‑151 (MQ=255)
               aTGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACAcc                                                                                                                               >  2:703478/1‑151 (MQ=255)
                   cGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATc                                                                                                                            >  2:748223/1‑150 (MQ=255)
                   cGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCg                                                                                                                           >  2:529657/1‑151 (MQ=255)
                    ggTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGc                                                                                                                          <  1:815787/151‑1 (MQ=255)
                       tGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCg                                                                                                                       >  1:478473/1‑151 (MQ=255)
                         aTACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGgt                                                                                                                     >  1:898104/1‑151 (MQ=255)
                            cccGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGgtg                                                                                                                    <  2:665610/149‑1 (MQ=255)
                             ccGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCt                                                                                                                 <  2:501466/151‑1 (MQ=255)
                             ccGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCt                                                                                                                 <  2:390746/151‑1 (MQ=255)
                              cGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTg                                                                                                                <  2:607741/151‑1 (MQ=255)
                                 aCTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAgg                                                                                                             >  2:521372/1‑151 (MQ=255)
                                       caGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAgg                                                                                                       >  1:914505/1‑151 (MQ=255)
                                       caGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAg                                                                                                        <  2:418999/150‑1 (MQ=255)
                                        aGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGt                                                                                                      >  1:9655/1‑151 (MQ=255)
                                           cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTAcg                                                                                                   <  1:387458/151‑1 (MQ=255)
                                           cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTAcg                                                                                                   >  1:566513/1‑151 (MQ=255)
                                            gTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTAcgc                                                                                                  <  2:75406/151‑1 (MQ=255)
                                            gTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTAcgc                                                                                                  <  1:230647/151‑1 (MQ=255)
                                            gTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTAcgc                                                                                                  >  2:337402/1‑151 (MQ=255)
                                            gTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTAcgc                                                                                                  >  2:123438/1‑151 (MQ=255)
                                               cGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACgcgg                                                                                                <  1:776944/150‑1 (MQ=255)
                                                   ttACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCt                                                                                           >  1:646068/1‑151 (MQ=255)
                                                   ttACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCt                                                                                           >  2:312936/1‑151 (MQ=255)
                                                           tGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCg                                                                                   >  1:208550/1‑151 (MQ=255)
                                                                gcgcGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGc                                                                               >  1:886612/1‑150 (MQ=255)
                                                                 cgcgTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCAt                                                                             <  1:480722/151‑1 (MQ=255)
                                                                  gcgTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATc                                                                            >  1:479471/1‑151 (MQ=255)
                                                                   cgTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCt                                                                           >  1:813386/1‑151 (MQ=255)
                                                                    gTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCtt                                                                          <  1:170414/151‑1 (MQ=255)
                                                                    gTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCtt                                                                          >  1:283449/1‑151 (MQ=255)
                                                                           gcgcTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGc                                                                   >  2:163038/1‑151 (MQ=255)
                                                                             gcTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCa                                                                  >  1:921774/1‑150 (MQ=255)
                                                                                    gAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                          <  2:97330/151‑1 (MQ=255)
                                                                                            gggATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATTGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaa                                                  <  1:769012/151‑1 (MQ=255)
                                                                                                      gCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                                        >  2:394246/1‑151 (MQ=255)
                                                                                                          gATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCg                                     >  1:264718/1‑150 (MQ=255)
                                                                                                           aTTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGAt                                   >  1:673651/1‑151 (MQ=255)
                                                                                                           aTTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGAt                                   <  2:721993/151‑1 (MQ=255)
                                                                                                                    gaCCCGCGACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                          >  2:268114/1‑151 (MQ=255)
                                                                                                                         gcgACGGACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg                     >  1:773726/1‑151 (MQ=255)
                                                                                                                               gACAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgctg               >  1:362754/1‑151 (MQ=255)
                                                                                                                                aCAGTGTGTGGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgctgc              <  2:226987/151‑1 (MQ=255)
                                                                                                                                     gtgtgGTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTg         <  1:353149/151‑1 (MQ=255)
                                                                                                                                         ggTAGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCg      <  1:789646/150‑1 (MQ=255)
                                                                                                                                            aGGCCCAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTg  >  2:671174/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCACTTGTCGGCTG  >  NC_000913/629402‑629692

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: