Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,436 G→A 100% H71Y (CAC→TAC)  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,4360GA98.6% 92.7 / ‑6.6 35H71Y (CAC→TAC) furDNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (13/19);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAA  >  NC_000913/710297‑710573
                                                                                                                                           |                                                                                                                                         
tttGCGGCAACTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGTAGATCAGGTGATCGTGGTGATGTTGCTGTGCCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCg                                                                                                                                        <  1:668645/143‑1 (MQ=255)
     ggCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCCGCTGCGGCGATCCTGGGGATCTTGCTGCGTCAGTTCAAATACGGATTTACCCCCTACACAATTGTGGCGGGGGa                                                                                                                                   >  2:714304/1‑143 (MQ=255)
      gCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGAc                                                                                                                                  >  1:452846/1‑143 (MQ=255)
         aTTTCACGCTGACGCGCTTCGATGGAATCATCACTAGAGTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATATTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGAt                                                                                                                               <  2:771997/143‑1 (MQ=255)
          tttCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGAGAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATa                                                                                                                              <  1:37555/143‑1 (MQ=255)
           ttCACGCGGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATAc                                                                                                                             <  1:778878/143‑1 (MQ=255)
                      gcgcTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGAGAGGTTGCTGGGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCa                                                                                                                  >  2:715831/1‑143 (MQ=255)
                              aTGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGtt                                                                                                          <  2:507200/143‑1 (MQ=255)
                              aTGGAATCATCACGAAATTCGACAACCTCGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGtt                                                                                                          <  1:304183/143‑1 (MQ=255)
                                ggAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATGTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCa                                                                                                        <  1:665355/143‑1 (MQ=255)
                                          cTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGata                                                                                              <  1:198756/143‑1 (MQ=255)
                                                       ccgtGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAATTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGa                                                                                 <  2:232192/140‑1 (MQ=255)
                                                        cTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCACCTTCAAAATTGTAGAGGGTGACGATACCAGCGTCGTAAAACTGGTTAAAACCGCGATATGACGTaagcaaac                                                                                >  2:3350/1‑136 (MQ=255)
                                                            ccGCAGTCGAGGCAGATCAGGTGCTCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAAt                                                                            <  1:749620/143‑1 (MQ=255)
                                                               ctgTCGAGCTAGATCATGTGATCGTGGTGATGTTGCTGTGTCTGTTCAACTCAGGCTTTCCCGTCCAGCAAATTGTAGCGGGTGACGATACCAGCGTCGTCACACTCGTACAGTACGCGATATACCGTAGCCAGACCAATttc                                                                         <  2:668717/141‑1 (MQ=255)
                                                                 gTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATttctt                                                                       >  1:424563/1‑143 (MQ=255)
                                                                  tCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCGCATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGGTATACCGTCGCCAGGCCAATttcttc                                                                      >  1:798157/1‑143 (MQ=255)
                                                                  tCGAGCCAGATCAGGTGCTCGTGGTGATGTTGCTGTGTGCGTTCAAATACGGATTTACCGCCTTTAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATttcttc                                                                      <  1:495454/143‑1 (MQ=255)
                                                                      ggCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCAccc                                                                  >  2:110766/1‑143 (MQ=255)
                                                                                    tCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATGGATCAGACg                                                    >  1:389706/1‑143 (MQ=255)
                                                                                       tggtgATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGttt                                                 <  1:866500/143‑1 (MQ=255)
                                                                                                  cTGTGTCAGTTCAAATACGGATTTACCGCCTTGAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCtt                                      <  2:25513/143‑1 (MQ=255)
                                                                                                   tgtgTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTc                                     >  1:375298/1‑143 (MQ=255)
                                                                                                       tCAGTTCAAAGACGGATTTACCGCCTTCAAAATTGTGGCGGGAGACCATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCa                                 <  2:387299/143‑1 (MQ=255)
                                                                                                           ttCAAAGACGGATTTGCCGCCTTCAAAATTGTAGCGGGTGACGTTCCCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGa                             <  1:254006/143‑1 (MQ=255)
                                                                                                            tCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGAc                            >  1:558738/1‑143 (MQ=255)
                                                                                                                  aCGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATg                      <  1:715831/143‑1 (MQ=255)
                                                                                                                  aCGGAGTCAGCGCCTTCAAAAGTGTGGGGGGTGACGATACCAGCGTCGTCAACCTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATg                      <  2:258738/143‑1 (MQ=255)
                                                                                                                         tACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATAGACAGTAGACAGACCAATTTCTTCACCCATATCGATCAGCCGTTTGTATAAATCATACGGACAGACGTGAAAGTTGTcg               >  2:808427/1‑142 (MQ=255)
                                                                                                                         tACCGCCTTCAAAAGTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGGAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTcc               <  1:821690/143‑1 (MQ=255)
                                                                                                                         tACCGCCTACAAAAGTGTAGCGGGTGACGATACCAGCGGCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTcc               <  1:491627/143‑1 (MQ=255)
                                                                                                                           ccGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGAACGCGATATACCGTAGACAGACCACGTTCTGCACCCATATCGATCAGACGTTTGTATCAATCTTCCGCACTTACTTGAGTGTTGTCCgc             >  1:504625/1‑142 (MQ=255)
                                                                                                                            cGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGc            >  1:319350/1‑143 (MQ=255)
                                                                                                                                     cactgtgGCGGGGTACCACACCAGCGGCGTCAAACTCGTTCAGTACGCGATATACCGTACCCCGACCAATTTCTGCACCCATATCGATCAGACGTTTGTATAAATCTCCCGCACTGACGTGATGGTTGTCCGTCTCCTgaaga   <  2:432004/140‑1 (MQ=255)
                                                                                                                                     aaTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCCGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgacga   >  1:223877/1‑143 (MQ=255)
                                                                                                                                      aTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa  <  1:476125/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                         
TTTGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAA  >  NC_000913/710297‑710573

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