Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,595 G→T 100% P18T (CCT→ACT)  fur ← DNA‑binding transcriptional dual regulator Fur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,5950GT92.9% 106.2 / ‑0.5 45P18T (CCT→ACT) furDNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (25/15);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01
Rejected as polymorphism: E-value score below prediction cutoff.

TCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAT  >  NC_000913/710459‑710733
                                                                                                                                        |                                                                                                                                          
tcAAACTGGTTCAGTATGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGACGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTCTTCAATCAGGAAGCGt                                                                                                                                      >  2:565227/1‑143 (MQ=255)
   accTGGTGCAGTCCGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTCACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTAc                                                                                                                                   <  2:10529/141‑1 (MQ=255)
      tGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGATTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACttt                                                                                                                                >  1:164232/1‑143 (MQ=255)
        ttgcAGTACTCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCa                                                                                                                              <  2:460361/140‑1 (MQ=255)
         ttCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAg                                                                                                                             >  1:547863/1‑143 (MQ=255)
         ttCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAg                                                                                                                             <  1:334269/143‑1 (MQ=255)
             gTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCCTATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCa                                                                                                                         <  2:190607/143‑1 (MQ=255)
                        aCCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAg                                                                                                              >  2:562327/1‑143 (MQ=255)
                                 aGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTGtt                                                                                                     >  1:351564/1‑143 (MQ=255)
                                  gACCAGTTTCTTCACCCCGAGCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTg                                                                                                    <  2:250081/143‑1 (MQ=255)
                                     cAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTa                                                                                                 >  1:9591/1‑143 (MQ=255)
                                          tcttcACCCGTATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTGGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGt                                                                                            <  2:282786/143‑1 (MQ=255)
                                          ccttcGCCCTCATCGATCTGACGTTTGTATAAATCTTCCGCACTGACGTGAGGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGt                                                                                            <  2:560226/142‑1 (MQ=255)
                                            ttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGCCGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCa                                                                                          >  1:277449/1‑143 (MQ=255)
                                            ttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGGCAGCTTTCTTTAGGGCGGTGTTGTTATCAGTCa                                                                                          >  1:290158/1‑143 (MQ=255)
                                            ttcACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTCGGGCGGTAGTGTTATCAGTCa                                                                                          >  2:119706/1‑143 (MQ=255)
                                               aCCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGc                                                                                       >  1:203484/1‑143 (MQ=255)
                                                 ccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTATGGCGGTATTGTTATCAGTCATACgg                                                                                     >  1:523062/1‑143 (MQ=255)
                                                 ccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCgg                                                                                     >  2:444868/1‑143 (MQ=255)
                                                 ccATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCgg                                                                                     >  1:143968/1‑143 (MQ=255)
                                                      tCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTCTTATCAGTAATGCGGAATCt                                                                                >  2:448280/1‑143 (MQ=255)
                                                             gACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGt                                                                         >  1:141177/1‑143 (MQ=255)
                                                             gACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGt                                                                         >  1:517272/1‑143 (MQ=255)
                                                              aCGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGtt                                                                        <  2:77157/143‑1 (MQ=255)
                                                                         aaTCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGtt                                                             <  2:371870/143‑1 (MQ=255)
                                                                           tCTTCCGCACTGCGGTGATGGTGGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTTCGTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCa                                                           <  2:357990/143‑1 (MQ=255)
                                                                           tCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATGTTTAAACGAGTAAGCGTCACTTTCAGGCAAGCTTACTTTAGGGAGGTAATGTTATCAGTCATCCGGACTCTGTCCCATTACAAcgcgcatca                                                           >  2:434473/1‑138 (MQ=255)
                                                                              tCCGCTCTGACGTGATGGTTGTCCGGCTCCTGAAGAACTGCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGACCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACtt                                                        <  2:630057/143‑1 (MQ=255)
                                                                               ccGCACGGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTAGTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTc                                                       <  1:27142/143‑1 (MQ=255)
                                                                                ctcccTGACGTGATGGGTGTCCGGCTCCTGAAGTACTTCCAGGCTTCTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCa                                                      <  2:386085/139‑1 (MQ=255)
                                                                                cGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCa                                                      >  2:4564/1‑143 (MQ=255)
                                                                                  cACTGACGTGTTGGTTGTCCGGCTCCTAAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTGTTGTTATCAGTCGTGCGGAGTCTGTCATGTTACTAAGCGGATCAcgcggct                                                    >  2:350599/1‑137 (MQ=255)
                                                                                  cACTGACGTGATTGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGTTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCAtt                                                    <  1:103786/143‑1 (MQ=255)
                                                                                  cACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTACAGGGCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAACCTGTCCTGTTACTAAGCGGTTCACATCAtt                                                    >  2:205184/1‑143 (MQ=255)
                                                                                            aTGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGGCATGCGGCATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                          >  2:95515/1‑143 (MQ=255)
                                                                                            aTGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGACAGCTTTCTGGAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTg                                          >  1:247383/1‑143 (MQ=255)
                                                                                              ggTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTGATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGAc                                        >  2:395178/1‑143 (MQ=255)
                                                                                              ggTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGAc                                        >  2:68135/1‑143 (MQ=255)
                                                                                                   tCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAAGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTGGAAGAAGTGAccggat                                   >  2:307861/1‑138 (MQ=255)
                                                                                                              aagaaCTTCCAGGATTTTTAAACGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTgag                        >  1:624655/1‑143 (MQ=255)
                                                                                                                 accTTCCAGGATTTTTAGACGAGTAAGCGTGACTTTCAGGCCGGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGGGACAGAATTTGCTCTTGAGATa                     <  2:283027/141‑1 (MQ=255)
                                                                                                                 aaCTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCAATTTAAGAATAGAAGTGACAGAATTTGCTCTTGAGATa                     <  2:203484/143‑1 (MQ=255)
                                                                                                                 aaCTTACAGGTTTGTTAAAGGAGTAAGCGTTACTTTCAGGCCAGCGTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATa                     <  2:280823/143‑1 (MQ=255)
                                                                                                                                tAAACGAGGAAGCGTTACTTGCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATtata      <  2:351564/143‑1 (MQ=255)
                                                                                                                                    cGAGTAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAt  >  1:454760/1‑143 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAACTTCCAGGATTTTTAAACGAGGAAGCGTTACTTTCAGGCCAGCTTTCTTTAGGGCGGTATTGTTATCAGTCATGCGGAATCTGTCCTGTTACTAAGCGGTTCACTTCATTAGAAGAAGTGACAGAATTTGCTCTTGAGATAATGCGTATCATTATAGAAT  >  NC_000913/710459‑710733

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