Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 839,746 C→T 9.8% E596K (GAA→AAA)  fiu ← catecholate siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913839,7460CT 9.8% 103.9 / 2.2 42E596K (GAA→AAA) fiucatecholate siderophore receptor
Reads supporting (aligned to +/- strand):  ref base C (15/20);  new base T (2/2);  total (20/22)
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.77e-01

TTTTTGATGGTTGCTTTTTGCTGGGTATAGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTGCGC  >  NC_000913/839604‑839890
                                                                                                                                              |                                                                                                                                                
tttttGATGGTTGCTTTTTGCTGGGTATAGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGGATTGCGAGTAAGTTCCGTCATCATTTTACTCAGCGTAATTTTc                                                                                                                                            >  1:106731/1‑149 (MQ=255)
tttttGATGGTTGCTTTTTACTGGGTATAGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTc                                                                                                                                            <  1:152731/149‑1 (MQ=255)
 ttttGATGGTTGCTTTTTGCTGGGTATAGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCCTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCg                                                                                                                                           <  1:802240/149‑1 (MQ=255)
     gATGGTTGCTTTTTGCTGGGTATAGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGGGTAAGTTCAGTCACCATCTTGCTCAACTTCCTTTTCTattt                                                                                                                                       >  1:299061/1‑147 (MQ=255)
      aTGGTTGCTTTTTGCTGGGTATATCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTTATTTTCGATATc                                                                                                                                      >  2:115161/1‑149 (MQ=255)
      aTGGTTGCTTTTTGCTGGGTATATCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATATCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTTATTTTCGATATc                                                                                                                                      >  2:115490/1‑149 (MQ=255)
                    cTGGCTATGGCCGCCAATCACCGGCCCCGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAgcg                                                                                                                        <  2:746370/149‑1 (MQ=255)
                       ggTATAGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTAATTTTCGATATCAGTGGGGAACGGCccgg                                                                                                                     >  1:393651/1‑149 (MQ=255)
                        gtagcgccgccAATCCCCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGTGGAACAGCgcggc                                                                                                                    <  2:634256/144‑1 (MQ=255)
                           taGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGt                                                                                                                 <  1:274521/149‑1 (MQ=255)
                           taGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGt                                                                                                                 >  2:384947/1‑149 (MQ=255)
                              ccgccAATCACCTGCCCCGCGGGAGTGATATTCCCGGCCACGGACATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAg                                                                                                              <  1:738636/149‑1 (MQ=255)
                              ccgccAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAg                                                                                                              <  1:738658/149‑1 (MQ=255)
                                  cAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAAc                                                                                                          <  1:828699/149‑1 (MQ=255)
                                   aaTCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGGGAGCAACa                                                                                                         >  1:441143/1‑149 (MQ=255)
                                     tCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCAATTTCGATATCAGAGCGGAACAGCGCGGCGGTGAGCAACAGa                                                                                                       >  2:187286/1‑149 (MQ=255)
                                           gCCACCCGGGAGGGATATGCNCNCCCNCGGATCTCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTa                                                                                                 <  1:457217/149‑1 (MQ=255)
                                            ccACGCGGGAGCACCATTCCCGGCCACGTATAGCCGATACCCTTCGACGCCTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTAt                                                                                                <  1:392978/149‑1 (MQ=255)
                                             cACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATc                                                                                               <  1:338598/149‑1 (MQ=255)
                                                  gggCGTGATATGCCCGGCCACGGATATCTCATAGCCCTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGaa                                                                                          <  2:449836/149‑1 (MQ=255)
                                                      gTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTg                                                                                      >  1:340340/1‑149 (MQ=255)
                                                                   ccACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTTATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCa                                                                         <  1:115161/149‑1 (MQ=255)
                                                                   ccACGGATACCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTTATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCa                                                                         <  1:115490/149‑1 (MQ=255)
                                                                      cGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGTGGCCAAtg                                                                      >  2:584230/1‑148 (MQ=255)
                                                                      cGGATATCTCATAGCCTTCGACGCGCTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTGGGTGCCAAtc                                                                      >  1:457923/1‑149 (MQ=255)
                                                                       gtgtatCTCATAGCCTTCGCCGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCCTTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAAtct                                                                     <  2:227209/146‑1 (MQ=255)
                                                                       ggATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCGTTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTACCagcct                                                                     >  2:438752/1‑145 (MQ=255)
                                                                              tcATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGCCGTTTATCGAGGAACTGCCCTTTGGTGCCAATCTCGCTGgt                                                              >  1:261720/1‑149 (MQ=255)
                                                                               cATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGgtg                                                             <  1:80976/149‑1 (MQ=255)
                                                                                         gACGCCTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCttttt                                                   <  2:72317/149‑1 (MQ=255)
                                                                                          aCGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAACCTCGCTGGTGTTGGCtttttc                                                  >  2:508455/1‑148 (MQ=255)
                                                                                          aCGCGTTTCCTACCGTAGTGCGAGTAAGTTCCTTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTg                                                  <  2:270942/149‑1 (MQ=255)
                                                                                            gcgTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTAATTTTCGATATCGGTGCGGAACATCGCGGCGGTGAGCAACCGACGTTTATCCAGAACCTGCCATTTGGTGCCCATCTCGCTGGTGTTGGCTTTTTggg                                                >  2:161265/1‑147 (MQ=255)
                                                                                            gcgTTTCTTACCGTATTGCGAGGAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCg                                                <  2:567825/149‑1 (MQ=255)
                                                                                              gTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGt                                              >  2:59904/1‑149 (MQ=255)
                                                                                                      ccctATTGCGAGTAATTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTACCCAGACCCTCACATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATc                                      <  2:379870/146‑1 (MQ=255)
                                                                                                                     gTTCCGTCACCATTTTGCTCAACCTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTGGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACt                       <  2:298340/149‑1 (MQ=255)
                                                                                                                              tcatTTTGCTCAACTTCATTTCCGATATCAGTGCGGAACAGCGCGGCGGCGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAACCTCGCCGGTGTTGGCTTTTTGCGGCTTAAACTCGGTGCGGTCGGCACTCTCAccccc              >  2:91410/1‑146 (MQ=255)
                                                                                                                                       tCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTg     >  1:29968/1‑149 (MQ=255)
                                                                                                                                       ccAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTg     <  1:149131/148‑1 (MQ=255)
                                                                                                                                          aCTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTgcgc  >  2:206271/1‑149 (MQ=255)
                                                                                                                                          aCTTCATTTTCGACATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACATGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCCGACtcggc  >  2:559656/1‑145 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TTTTTGATGGTTGCTTTTTGCTGGGTATAGCCGCCAATCACCTGCCACGCGGGAGTGATATTCCCGGCCACGGATATCTCATAGCCTTCGACGCGTTTCTTACCGTATTGCGAGTAAGTTCCGTCATCATTTTGCTCAACTTCATTTTCGATATCAGTGCGGAACAGCGCGGCGGTGAGCAACAGACGTTTATCCAGAACCTGCCATTTGGTGCCAATCTCGCTGGTGTTGGCTTTTTGCGGTTTAAAATCGGTGCGGTTGGCACTGTTACCGCTGCCAGACTGCGC  >  NC_000913/839604‑839890

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑