Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,512,350:1 +G 19.6% coding (65/957 nt) fecD ← ferric citrate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,3501.G19.6% 128.8 / 40.3 46coding (65/957 nt)fecDferric citrate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base . (16/21);  new base G (7/2);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAG  >  NC_000913/4512202‑4512495
                                                                                                                                                     |                                                                                                                                                 
cAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCatat                                                                                                                                                    <  1:352386/149‑1 (MQ=255)
 aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑g                                                                                                                                                  <  3:388458/149‑1 (MQ=255)
  gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCatat                                                                                                                                                    <  3:221796/147‑1 (MQ=255)
  gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCatat                                                                                                                                                    >  4:221796/1‑147 (MQ=255)
  gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑gg                                                                                                                                                 >  1:140453/1‑149 (MQ=255)
  gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑gg                                                                                                                                                 >  3:163603/1‑149 (MQ=255)
  gCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑gg                                                                                                                                                 >  4:36663/1‑149 (MQ=255)
    aCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAg                                                                                                                                               <  2:465205/149‑1 (MQ=255)
        ccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGAt                                                                                                                                           <  2:213791/149‑1 (MQ=255)
                       cTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAAAACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGccc                                                                                                                            <  3:113396/149‑1 (MQ=255)
                         gCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGc                                                                                                                          >  1:298362/1‑149 (MQ=255)
                            ccGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAg                                                                                                                       <  2:298362/149‑1 (MQ=255)
                                   aCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAgg                                                                                                                >  3:50011/1‑149 (MQ=255)
                                   aCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAgg                                                                                                                >  3:538878/1‑149 (MQ=255)
                                   aCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAgg                                                                                                                <  4:50011/149‑1 (MQ=255)
                                          cagcagCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCAACGGCACGGGGAACACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatga                                                                                                         >  2:122397/1‑149 (MQ=255)
                                           agcagcAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaa                                                                                                        >  2:167548/1‑149 (MQ=255)
                                            gcagcaAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaa                                                                                                       <  2:128418/149‑1 (MQ=255)
                                             cagcaAGCGCGGCAGTCGGTACTCCATCAATACATCATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaaaa                                                                                                      <  3:400095/149‑1 (MQ=255)
                                                    cgcgGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAg                                                                                               <  3:33128/149‑1 (MQ=255)
                                                          aGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAt                                                                                         <  1:523474/149‑1 (MQ=255)
                                                              ggTACTCCATCAATACATAATAATGCTCGTTTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTc                                                                                     <  3:161677/149‑1 (MQ=255)
                                                                     cATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCg                                                                              >  2:70126/1‑149 (MQ=255)
                                                                     cATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGAGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCg                                                                              >  3:293324/1‑149 (MQ=255)
                                                                      aTCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGc                                                                             >  2:287675/1‑149 (MQ=255)
                                                                        cAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcct                                                                           <  1:527231/149‑1 (MQ=255)
                                                                         aaTACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctc                                                                          <  4:383264/149‑1 (MQ=255)
                                                                              ataataATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGTGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcac                                                                     >  3:363479/1‑149 (MQ=255)
                                                                                aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACaa                                                                   <  3:152745/149‑1 (MQ=255)
                                                                                aataatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACaa                                                                   <  4:538878/149‑1 (MQ=255)
                                                                                   aatGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGTGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGcc                                                                <  4:363479/149‑1 (MQ=255)
                                                                                      gCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGa                                                             >  2:413953/1‑149 (MQ=255)
                                                                                        tCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                           >  1:374246/1‑149 (MQ=255)
                                                                                               cGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAg                                                    <  3:400117/149‑1 (MQ=255)
                                                                                                ggCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAACAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAgg                                                   >  2:359978/1‑149 (MQ=255)
                                                                                                             gTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAgc                                      >  2:52794/1‑149 (MQ=255)
                                                                                                                    gcgcgcgcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCa                               >  3:186192/1‑149 (MQ=255)
                                                                                                                        gcgcCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgc                           >  3:32378/1‑149 (MQ=255)
                                                                                                                            cacggcacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCt                       <  2:537519/149‑1 (MQ=255)
                                                                                                                                gcacggGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACATGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCgg                   >  3:482916/1‑149 (MQ=255)
                                                                                                                                     gggATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGa              <  3:47647/149‑1 (MQ=255)
                                                                                                                                      ttaTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctc                                                                          >  2:61435/3‑88 (MQ=38)
                                                                                                                                      ttaTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGcctc                                                                          <  1:61435/86‑1 (MQ=38)
                                                                                                                                      ttaTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAt             >  1:484364/3‑149 (MQ=255)
                                                                                                                                      ggATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAt             <  1:287675/149‑1 (MQ=255)
                                                                                                                                       taTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCCGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtc            >  3:126966/2‑149 (MQ=255)
                                                                                                                                       taTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtc            >  1:517448/2‑149 (MQ=255)
                                                                                                                                        aTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtct           >  4:316328/1‑149 (MQ=255)
                                                                                                                                        aTCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATAAGCGCCAGCACTGCGCCTGCGGGCAGAtct           >  1:333913/1‑149 (MQ=255)
                                                                                                                                         tCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCa                                  >  3:454054/1‑125 (MQ=255)
                                                                                                                                         tCACT‑CCATATGGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCa                                  <  4:454054/125‑1 (MQ=255)
                                                                                                                                            cTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCgg       <  1:70126/149‑1 (MQ=255)
                                                                                                                                             tCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCggg      <  3:150846/149‑1 (MQ=255)
                                                                                                                                                 atat‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAg  <  1:399941/149‑1 (MQ=255)
                                                                                                                                                 atat‑GGAGTGATAACAGCGCACAGCCCGCCAGAGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAg  <  4:293324/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAT‑GGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAG  >  NC_000913/4512202‑4512495

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