Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 626,267 C→T 9.3% R66R (CGC→CGT)  entE → 2,3‑dihydroxybenzoate‑AMP ligase component of enterobactin synthase multienzyme complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913626,2670CT 9.3% 124.5 / 3.3 43R66R (CGC→CGT) entE2,3‑dihydroxybenzoate‑AMP ligase component of enterobactin synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base C (20/19);  new base T (2/2);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

GGGAAAAAGGCTACTGGCAGGATTTGCCGCTGACCGACATTCTGACGCGACATGCTGCGAGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCATCAGCGTAGTGAACTGAAC  >  NC_000913/626119‑626408
                                                                                                                                                    |                                                                                                                                             
gggAAAAAGGCTACTGGCAGGATTTGCCGCTGACCGACATTCTGACGCGACATGCTGCGAGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGc                                                                                                                                               <  5:200285/149‑1 (MQ=255)
           tACTGGCAGGATTTGCCGCTGACCGACATTCTGACGCGACATGCTGCGAGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCAt                                                                                                                                    >  1:221346/1‑149 (MQ=255)
                     aTTTGCCGCTGACCGACATTCTGACGCGACATGCTGCGAGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGt                                                                                                                          <  8:126569/149‑1 (MQ=255)
                                  cGACATTCTGACGCGACATGCTGCGAGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTgg                                                                                                             >  2:168519/1‑149 (MQ=255)
                                                       tgcGAGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCt                                                                                         >  8:105643/1‑148 (MQ=255)
                                                       tgcGAGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCt                                                                                         >  4:159369/1‑148 (MQ=255)
                                                          gAGTGACAGCATCGCCGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCGGGTACAACTGGGTAACGTCGCTGAAt                                                                                     >  4:135192/1‑149 (MQ=255)
                                                           aGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAAtt                                                                                    <  4:107927/149‑1 (MQ=255)
                                                            gTGACAGCATCGCGGTTATTGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTATGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTg                                                                                   >  2:128132/1‑149 (MQ=255)
                                                                cAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGtata                                                                               <  3:148938/149‑1 (MQ=255)
                                                                   cATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTa                                                                            >  5:177381/1‑149 (MQ=255)
                                                                   cATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTa                                                                            >  1:75797/1‑149 (MQ=255)
                                                                         ggTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCTTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCtttt                                                                      <  3:159369/149‑1 (MQ=255)
                                                                                   ggCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTa                                                                                                           >  2:193441/1‑102 (MQ=255)
                                                                                   ggCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTa                                                                                                           <  1:193441/102‑1 (MQ=255)
                                                                                        gcgaCAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATAtt                                                                             >  5:103825/1‑127 (MQ=255)
                                                                                        gcgaCAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATAtt                                                                             <  6:103825/127‑1 (MQ=255)
                                                                                          gaCAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTAcc                                                                          <  1:84602/128‑1 (MQ=255)
                                                                                          gaCAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTAcc                                                                          >  2:84602/1‑128 (MQ=255)
                                                                                             aGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGc                                                  <  1:168519/149‑1 (MQ=255)
                                                                                              gTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCg                                                 >  1:74761/1‑149 (MQ=255)
                                                                                               ttGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGt                                                >  1:136361/1‑149 (MQ=255)
                                                                                                    ttATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTgcgc                                           <  3:91297/149‑1 (MQ=255)
                                                                                                            aGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTg                                   <  2:221346/149‑1 (MQ=255)
                                                                                                               tGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGta                                                                                 <  2:83861/100‑1 (MQ=255)
                                                                                                                  aTCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCgttg                             <  2:75797/149‑1 (MQ=255)
                                                                                                                     aGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGTGTTGCGCCGGTGCTGGCgttgtt                           >  8:78333/1‑148 (MQ=255)
                                                                                                                      ggcgtcGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCa                         <  5:105715/149‑1 (MQ=255)
                                                                                                                      ggcggcGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCGGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCa                         <  1:49338/149‑1 (MQ=255)
                                                                                                                            ggATAACCTCGCGTGTAGTTTATGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCATc                   <  1:128132/149‑1 (MQ=255)
                                                                                                                            ggATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCAt                    >  6:95914/1‑148 (MQ=255)
                                                                                                                                 aCCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACC‑TTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCATCAGCGTa             <  6:177381/149‑1 (MQ=255)
                                                                                                                                   cTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCATCAGCGTAg            >  7:31020/1‑149 (MQ=255)
                                                                                                                                      gcgTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCg                                         >  2:62871/1‑117 (MQ=255)
                                                                                                                                      gcgTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCg                                         <  1:62871/117‑1 (MQ=255)
                                                                                                                                        gtgtAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCATCATCGTAGtgaac       >  1:87988/1‑149 (MQ=255)
                                                                                                                                            aGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCATCAGCGTAGtgaactgaa   >  4:203868/1‑149 (MQ=255)
                                                                                                                                             gTTTACGTCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTgg                                                    <  2:87659/99‑1 (MQ=255)
                                                                                                                                             gTTTACGTCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTgg                                                    <  8:157291/99‑1 (MQ=255)
                                                                                                                                             gTTTACGTCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTgg                                                    >  7:157291/1‑99 (MQ=255)
                                                                                                                                             gTTTACGTCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTgg                                                    >  1:87659/1‑99 (MQ=255)
                                                                                                                                             gTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCATCAGCGTAGtgaactgaac  <  2:136361/149‑1 (MQ=255)
                                                                                                                                              tttACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCg                                <  7:249658/118‑1 (MQ=255)
                                                                                                                                              tttACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCg                                >  8:249658/1‑118 (MQ=255)
                                                                                                                                                    ccGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGtat                                                                                >  4:35763/1‑64 (MQ=255)
                                                                                                                                                    ccGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGtat                                                                                <  3:35763/64‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGGAAAAAGGCTACTGGCAGGATTTGCCGCTGACCGACATTCTGACGCGACATGCTGCGAGTGACAGCATCGCGGTTATCGACGGCGAGCGACAGTTGAGTTATCGGGAGCTGAATCAGGCGGCGGATAACCTCGCGTGTAGTTTACGCCGTCAGGGCATTAAACCTGGTGAAACCGCGCTGGTACAACTGGGTAACGTCGCTGAATTGTATATTACCTTTTTCGCGCTGCTGAAACTGGGCGTTGCGCCGGTGCTGGCGTTGTTCAGCCATCAGCGTAGTGAACTGAAC  >  NC_000913/626119‑626408

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