Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,513,841 A→G R156R (CGT→CGC fecB ← ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,513,8410AG93.5% 79.9 / ‑2.4 32R156R (CGT→CGCfecBferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (14/15);  total (17/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-02

TGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGG  >  NC_000913/4513602‑4514078
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
ggggCACCCAGCAGCAGGGCAGGATTGACCGCCAGCATTTGCTCCAGGCCGAGGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCAt                                                                                                                                                                                                                                    <  1:174423/250‑1 (MQ=255)
                        ttGACCGCCAGCAGTTGCCCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCAttt                                                                                                                                                                                                             <  1:167334/250‑1 (MQ=255)
                           aCCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGcc                                                                                                                                                                                                         <  1:203798/251‑1 (MQ=255)
                            ccGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCg                                                                                                                                                                                                        >  1:116426/1‑251 (MQ=255)
                               ccAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCCGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatg                                                                                                                                                                                                     >  1:91052/1‑251 (MQ=255)
                               ccAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatg                                                                                                                                                                                                     >  1:13380/1‑251 (MQ=255)
                                cagcagTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatg                                                                                                                                                                                                     >  1:16377/1‑250 (MQ=255)
                                  gcatttGCCCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATa                                                                                                                                                                                                  <  1:282213/247‑1 (MQ=255)
                                  gcagTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCgatgat                                                                                                                                                                                                   >  1:63176/1‑250 (MQ=255)
                                             aGGCCGATGGACGGCCTGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCCCGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCACTTTGTGAAGCTGGCTGGCCCCCTGCTCCATCCTCCCTTTATGTTGTTCCTGTCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGTTAGCAGCAGATTg                                                                                                                                                                                       >  2:113472/1‑251 (MQ=255)
                                             aGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTg                                                                                                                                                                                       >  1:31664/1‑251 (MQ=255)
                                                ccaggggccggcagggtcgaccccGCCATCGCAGCGGGAACGGTCAGCCCCAGAGAGGCCTGCACGCTGCCGGTCCAGGCCTCCTGAGTATGCAGGTTGAATTGCGGTTCGCTTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCGGGCTGGCCCCCTGCGCCTTCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCaa                                                                                                                                                                                    <  2:284801/230‑1 (MQ=255)
                                                ccGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCa                                                                                                                                                                                     >  1:62335/1‑250 (MQ=255)
                                                     ggatggCAGGGACGCCCCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCGTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAAtttt                                                                                                                                                                               <  1:204491/247‑1 (MQ=255)
                                                                                  cGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGGCTCCTGAGTATTCAGGTTGAATTGCTGGTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCGGGTCGGCCCACTGCGCCAACCTCTTTTTATGTCGTTCCAGGCTGGCCTGCATCTCTCGCTTTTTCCCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAagc                                                                                                                                                  <  1:350804/251‑1 (MQ=255)
                                                                                                    gagagGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGAtt                                                                                                                                <  1:79122/251‑1 (MQ=255)
                                                                                                         gCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTtgctg                                                                                                                           >  1:151840/1‑251 (MQ=255)
                                                                                                                 gCTGACGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTCGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCg                                                                                                                    >  1:7261/1‑250 (MQ=255)
                                                                                                                    gCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTTCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATg                                                                                                                 <  1:331468/250‑1 (MQ=255)
                                                                                                                     ccGGTTCAGGTCTCCTGGGTATGTAGGTGGAATTGCTGTTCGGGTGATGTGCCAAAGGCCACGCTTGTCCCTTTTGGAAGCTGGCTGGCCCACTGCGCCTTCCGCGCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTGTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTa                                                                                                               <  2:110937/251‑1 (MQ=255)
                                                                                                                                       gTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGAct                                                                                             <  1:203828/251‑1 (MQ=255)
                                                                                                                                               gTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCgg                                                                                     >  1:79995/1‑251 (MQ=255)
                                                                                                                                                     ttGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATg                                                                                >  1:145078/1‑250 (MQ=255)
                                                                                                                                                             cgcgTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTc                                                                        <  1:272828/250‑1 (MQ=255)
                                                                                                                                                              gcttgaattgccaatggctccacgcgtCCTTTTGGGAAGCTGGCCTGCCCACGTCTCCATCCTCTCTTTATTTTTTTCCAGGCGTGCCTGAATCTATCGCGTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTTCAAATTTTCAGCGTATTTTTCGTTGCGGGACTTAATCAGCAGTAGCGGCGCGATTGGCTGCAAGGCGATGTAACCCCCCGCATTGCTACTGCTGTCGGCAATGATCAGGTCTg                                                                      <  2:357866/226‑1 (MQ=255)
                                                                                                                                                                 tGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCTTCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTCCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGtt                                                                   >  1:88355/1‑251 (MQ=255)
                                                                                                                                                                           aaGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGc                                                         >  1:67206/1‑251 (MQ=255)
                                                                                                                                                                            aGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCa                                                        <  1:113472/251‑1 (MQ=255)
                                                                                                                                                                               ccACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTCCCCACCATTTCGCCGATGATAGCCGCCGATTGCCAATTTTCAGCGTAGGTTTCGTTGCAGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGGTGTAAACCCCCGCATGGCGACTGCTGTCGGCACTGATTCGGTCTGGTTTAGGAGCGGCCATg                                                     >  2:321968/1‑251 (MQ=255)
                                                                                                                                                                                                      tggctggcCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTgcgcgcg                              >  1:260429/1‑251 (MQ=255)
                                                                                                                                                                                                                          actctctTTATTTTGTTCCATGCGTGCCTGCAGCTCTCGCTTTTTTCCCGCCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGTGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCTTCAATGTCTGCCAGGCTCGTCTGCGCGCGTGTTACGACCGACTGCCAc           <  2:361210/249‑1 (MQ=255)
                                                                                                                                                                                                                                  tATGTTGTTCCAGGCGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGCTCCAACGGACTGCCACGGTTTCAgg  >  1:210848/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                             
TGGGCAACCAGCAGCCAGGCAGGATTGACCGCCAGCAGTTGCTCCAGGCCGATGGACGGCATGGACGCACCCGCCATCGCAGCGGGAACGTTCAGCCCCAGAGAGGCCAGCACGCTGCCGGTCCAGGTCTCCTGAGTATGCAGGTTGAATTGCTGTTCGCGTGATGTGCCAAAGGCCACGCGTGTCCCTTTGGGAAGCTGGCTGGCCCACTGCGCCATCCTCTCTTTATGTTGTTCCAGACGTGCCTGCATCTCTCGCTTTTTACCCACCATTTCGCCGATGATAGCCGCAGATTGCAAATTTTCAGCGTAGGTTTCGTTGCGGGACTTAAGCAGCAGTACCGGCGCGATTTGCTGCAAGGCGATGTAAACCCCCGCATGGCGACTGCTGTCGGCAATGATCAGGTCTGGTTTCAGAGCGGCAATGGCTTCCAGGCTCGGCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGG  >  NC_000913/4513602‑4514078

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