Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 624,659 A→T 6.7% intergenic (‑149/‑226) fepB ← / → entC ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913624,6590AT 6.7% 208.1 / 5.2 60intergenic (‑149/‑226)fepB/entCferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base A (31/25);  new base T (2/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

TTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTG  >  NC_000913/624519‑624805
                                                                                                                                            |                                                                                                                                                  
ttCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                                                                                            <  2:370431/149‑1 (MQ=255)
  ccTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAgg                                                                                                                                          <  4:549111/149‑1 (MQ=255)
      ttattaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                                                                                                      >  2:610716/1‑149 (MQ=255)
          taataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAgggg                                                                                                                                  <  2:548107/149‑1 (MQ=255)
            ataaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGc                                                                                                                                >  4:137988/1‑149 (MQ=255)
             taaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCTCCCTTCCAGGGTGAgggctgg                                                                                                                               <  2:551629/149‑5 (MQ=25)
                 gTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAAc                                                                                                                           <  1:432348/149‑1 (MQ=255)
                  ttAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGc                                                                                                                                >  1:429272/1‑143 (MQ=255)
                  ttAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGc                                                                                                                                <  2:429272/143‑1 (MQ=255)
                   tAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACgg                                                                                                                         <  4:591305/149‑1 (MQ=255)
                    aaGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgc                                                                                                                        <  1:256865/149‑1 (MQ=255)
                        gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                                                                                    >  2:69272/1‑149 (MQ=255)
                        gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAg                                                                                                                    <  3:331569/149‑1 (MQ=255)
                           tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACgg                                                                                                                         <  1:143278/141‑1 (MQ=255)
                           tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACgg                                                                                                                         >  2:143278/1‑141 (MQ=255)
                           tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGAc                                                                                                                 >  2:193034/1‑149 (MQ=255)
                              tGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATc                                                                                                              <  3:63232/149‑1 (MQ=255)
                               gACAAATTCGACAAAGCGCACAATCCGTCCCCTCGACCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATca                                                                                                             >  3:339332/1‑149 (MQ=255)
                                     ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgc                                                                                                       >  1:499367/1‑149 (MQ=255)
                                     ttCCACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAGTGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAAGACATCACATTgc                                                                                                       >  4:392916/1‑149 (MQ=255)
                                       cGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcgc                                                                                                     <  3:392388/149‑1 (MQ=255)
                                                    aaTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGTCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCa                                                                                        >  1:63697/1‑149 (MQ=255)
                                                            cccTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaataat                                                                                >  3:409746/1‑149 (MQ=255)
                                                                   cccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgc                                                                                                                        <  3:417570/102‑1 (MQ=255)
                                                                   cccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgc                                                                                                                        >  4:417570/1‑102 (MQ=255)
                                                                   cccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                                                                                               <  3:425012/111‑1 (MQ=255)
                                                                   cccTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACAt                                                                                                               >  4:425012/1‑111 (MQ=255)
                                                                      tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                                                                                       <  4:300406/132‑1 (MQ=255)
                                                                      tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                                                                                       >  3:300406/1‑132 (MQ=255)
                                                                                gggCTAGGGTGAGGGCAACAGCCAGCACTGGAGCGAACATTAACCCTCACACCAGACCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCAATCAATAATGCTTCTCATTTTCATTGTaa                                                            <  2:191357/149‑1 (MQ=255)
                                                                                 ggTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAAc                                                           <  2:589752/149‑1 (MQ=255)
                                                                                  gTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAAcc                                                          >  2:221203/1‑149 (MQ=255)
                                                                                   ttAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTAt                                                                                                 >  4:529892/1‑109 (MQ=255)
                                                                                   ttAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTAt                                                                                                 <  3:529892/109‑1 (MQ=255)
                                                                                     aGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTg                                                                                                        <  4:564262/100‑1 (MQ=255)
                                                                                     aGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTg                                                                                                        >  3:564262/1‑100 (MQ=255)
                                                                                         tGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCa                                                   >  4:180440/1‑149 (MQ=255)
                                                                                           aGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGa                                                 >  2:572023/1‑149 (MQ=255)
                                                                                           aGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGa                                                 >  3:604563/1‑149 (MQ=255)
                                                                                              ggAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGc                                              >  1:421273/1‑149 (MQ=255)
                                                                                               gAACAGCCAGCACTGATGCGAACATTAACCATCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCACAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCACAGTTTCATTGTAACCACAACCAGATGCa                                             >  4:519903/1‑149 (MQ=255)
                                                                                                 aCAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAAc                                           <  3:153531/149‑1 (MQ=255)
                                                                                                     ccagcACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCg                                       >  1:31334/1‑149 (MQ=255)
                                                                                                     ccagcACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCg                                       >  4:311103/1‑149 (MQ=255)
                                                                                                                  gAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGc                          <  1:588660/149‑1 (MQ=255)
                                                                                                                    aCATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGt                        >  2:351479/1‑149 (MQ=255)
                                                                                                                    aCATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGt                        <  1:610716/149‑1 (MQ=255)
                                                                                                                      gttAACCCTCACCACAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTAAGCTAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTCCGTTa                      >  4:401585/2‑149 (MQ=255)
                                                                                                                       ttAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTAc                     >  4:535488/1‑149 (MQ=255)
                                                                                                                        tAACCCTCTCCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCTAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTAcc                    <  3:519903/149‑1 (MQ=255)
                                                                                                                            ccTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCaa                >  2:73893/1‑149 (MQ=255)
                                                                                                                             cTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAag               >  2:216394/1‑149 (MQ=255)
                                                                                                                                   ccAGCCCTCACCCTGGCAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACTCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGa         <  3:392916/149‑1 (MQ=255)
                                                                                                                                   ccAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGa         <  1:402026/149‑1 (MQ=255)
                                                                                                                                       cccTCTCCCTTCCAGGGTGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCTAATcc                                                                                         <  3:323518/65‑1 (MQ=25)
                                                                                                                                       cccTCTCCCTTCCAGGGTGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                                                                                         >  4:323518/1‑65 (MQ=25)
                                                                                                                                       cccTCTCCCTTCCAGGGTGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATa     >  2:152071/1‑149 (MQ=255)
                                                                                                                                          tCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTg  >  2:321465/1‑149 (MQ=255)
                                                                                                                                           cACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAAc                                                     >  2:222639/1‑97 (MQ=255)
                                                                                                                                           cACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAAc                                                     <  1:222639/97‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
TTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTG  >  NC_000913/624519‑624805

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