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,4360GA100.0% 90.9 / ‑6.7 33H71Y (CAC→TAC) furDNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/17);  minor base C (1/0);  total (16/17)
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.

TTTTGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGA  >  NC_000913/710296‑710572
                                                                                                                                            |                                                                                                                                         
ttttGCGGCAATTTCTCGCTGACGCGCTTCGATGGAATCATCACTAAATGCGATAACCTTGCCGCAGTAGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGc                                                                                                                                         <  2:119325/143‑1 (MQ=255)
      ggCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAACTGTCGCGGGTGa                                                                                                                                   >  1:645873/1‑143 (MQ=255)
       gCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATCCGATACCCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGAc                                                                                                                                  <  1:327089/143‑1 (MQ=255)
           tttCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCCGTCGAGGCAGATCACGTGATCGTGGTGATGTTGCTGTGTCAGTTCACATACGGATTTCCCCCCTTCAAAATTGTAGCCGGTGACGATa                                                                                                                              >  1:249969/1‑143 (MQ=255)
             tCACGCTGACGCGCTTCGATGGGATCATCACTAAATGCGATGACCTTGCGGCAGTCGGGGCAGATCAGGTGATCGTGGTGATGTGGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGGGACGATAcc                                                                                                                            >  1:496929/1‑143 (MQ=255)
               aCGCTGACGCGCTTCGATGGAATCATCACTAAATTCGCTAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAg                                                                                                                          <  2:44506/143‑1 (MQ=255)
               aCGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCGGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTAAGTTCAACTAGGGATTTACCGTCGTCAAAAGTGTAGCGGGTGGCGATACCAg                                                                                                                          >  2:606527/1‑143 (MQ=255)
                       gcgcTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCa                                                                                                                  >  1:198172/1‑143 (MQ=255)
                           ttCGCTGTAATCATCACTAAATTCGATACCCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACt                                                                                                              <  1:145588/143‑1 (MQ=255)
                             cGATGGAATCATCACTAAATTCGATAACCTTGGCGCAGTCGAGGCAGATCTGGCCCTTCGCGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTgg                                                                                                           <  1:490210/143‑1 (MQ=255)
                             cGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTgg                                                                                                           <  1:574322/143‑1 (MQ=255)
                               aTGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGtt                                                                                                         >  1:459937/1‑143 (MQ=255)
                                 ggATTCATCATTAAATTCGATAACCTTGCCGTAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTCCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCa                                                                                                       <  2:338101/143‑1 (MQ=255)
                                 ggAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCa                                                                                                       <  1:256763/143‑1 (MQ=255)
                                                    aTAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGAAGATACCAGCGTCGTCAAACTGGTTCACGACGCGATATACCGTAGc                                                                                    >  2:247102/1‑143 (MQ=255)
                                                      agctttGCCGCAGTCGCGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCa                                                                                  <  2:604845/139‑1 (MQ=255)
                                                         cTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGAc                                                                               >  1:215614/1‑143 (MQ=255)
                                                              cGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATAACAGCGTCGTAAAACTGGTTCAGTACGCGCTATACCGTAGCCAGACCaaaa                                                                          >  2:88714/1‑141 (MQ=255)
                                                                           gATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAACATTGTAGCGGGTGACGATACCAGCGTCGTCAAACAGGTTCAGTACGCGATATACCGTAGCCCGACCGATTTCTTCACCCatat                                                             >  2:127051/1‑143 (MQ=255)
                                                                                 gTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGCTCa                                                       >  1:574500/1‑143 (MQ=255)
                                                                                      cggggtgATGTTGTTGTGTCAGTTCACATACGGGTTTACCGCCTTCAAACTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGt                                                  <  2:361709/140‑1 (MQ=255)
                                                                                            gATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGtat                                            >  1:325881/1‑143 (MQ=255)
                                                                                                  gCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTAGTTCAGTACGCGATATACCGGAGGCAGACCAATTTCTTCACCGATAGCGCGCAGACGCTAGCATaaagat                                      >  2:63951/1‑140 (MQ=255)
                                                                                                     gtgtAAGTTCAAATAGGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATAGATCAGACGTTTGTATAAATCTTcc                                   <  1:156100/143‑1 (MQ=255)
                                                                                                               aaaGACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGt                         <  1:106277/143‑1 (MQ=255)
                                                                                                                 aTACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGa                       <  2:198172/143‑1 (MQ=255)
                                                                                                                    cGGATTTAGCGCCTTCAATATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATgg                    <  2:392868/143‑1 (MQ=255)
                                                                                                                     ggATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGt                   >  2:496196/1‑143 (MQ=255)
                                                                                                                      taTTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGtt                  <  2:215614/142‑1 (MQ=255)
                                                                                                                             cGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGc           <  2:459937/143‑1 (MQ=255)
                                                                                                                              gCCTTCCAAATGGTAGCGGGTGACGATACCCGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCt          <  1:198815/143‑1 (MQ=255)
                                                                                                                                    aaaaTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgac    >  1:469905/1‑142 (MQ=255)
                                                                                                                                      aaTTGTAGCGGGTGACGATACCCGCCTATCAAAATACGTACAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGAGGGTTGTCCGGCTCCTgagga  >  2:205894/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                         
TTTTGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGA  >  NC_000913/710296‑710572

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: