Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 710,436 G→A 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,4360GA96.4% 77.5 / NA 28H71Y (CAC→TAC) furDNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/13);  minor base C (1/0);  total (15/13)
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

TTGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAA  >  NC_000913/710298‑710573
                                                                                                                                          |                                                                                                                                         
ttGCGGCTACTTCCCGCTGACGCGCCTCGATGCAATCATCACTAAAGTCGATAATCTTGCCGCAGTCGAGCCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCgg                                                                                                                                       <  1:135866/143‑1 (MQ=255)
          ttCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGCTCGTGGTGATGTTGCTGTGTCAATTCAAATACGGATTTACCGCCTCCTACATTGTAAGGGCTGCCGATAc                                                                                                                             >  2:416876/1‑143 (MQ=255)
                    cgcgCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCgtcgtc                                                                                                                   <  2:600018/143‑1 (MQ=255)
                     gcgcTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTCACGGGTGACGATACCAACGTCGTCa                                                                                                                  >  2:123471/1‑143 (MQ=255)
                             aTGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGGTACCAGCGTCGTCAAACTGGtt                                                                                                          >  2:314800/1‑143 (MQ=255)
                             aTGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGtt                                                                                                          >  2:477116/1‑143 (MQ=255)
                             aTGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGATGACGATACCACCGGCGTCAAACTgact                                                                                                          >  2:514318/1‑140 (MQ=255)
                              tGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAACATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGGTc                                                                                                         >  2:28974/1‑143 (MQ=255)
                               ggAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCa                                                                                                        >  1:284942/1‑143 (MQ=255)
                                       caCTAAATGCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCGGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGa                                                                                                <  2:414795/143‑1 (MQ=255)
                                             aTTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATAcc                                                                                          <  1:264699/143‑1 (MQ=255)
                                                cGATAACCTTGCCGCAGTCGAGGCAGATCAGCTGACCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTa                                                                                       <  1:660730/143‑1 (MQ=255)
                                                       cTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCGAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGAc                                                                                >  1:343434/1‑143 (MQ=255)
                                                       cTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATGTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGAc                                                                                <  2:484141/143‑1 (MQ=255)
                                                           ccGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAAt                                                                            <  1:130438/143‑1 (MQ=255)
                                                             gCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGCTACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAAttt                                                                          <  1:399003/143‑1 (MQ=255)
                                                              cAGTCCGGGCAGATCAGGTGAGCGTGGTGATGGTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATttc                                                                         <  1:234270/143‑1 (MQ=255)
                                                                                   tCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACg                                                    >  2:429718/1‑143 (MQ=255)
                                                                                      tggtgATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGGAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGttt                                                 <  1:308503/143‑1 (MQ=255)
                                                                                          gATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGCTCAGCCGTTTGtat                                             >  2:600175/1‑143 (MQ=255)
                                                                                                          ttCAAATACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATATCTTCACCTATATCGATCAGACGGTGGTATAAATCTTCCGCACTGa                             >  1:371965/1‑143 (MQ=255)
                                                                                                                tACGGATTTACCGCCTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGAt                       >  1:328980/1‑143 (MQ=255)
                                                                                                                        tACCGCCTTCAAACATGTAGCGGATGACGAGACCAGCGTCGCCAAACTCGCTCAGTACGCGATATACCCTACCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTcc               <  1:14037/143‑1 (MQ=255)
                                                                                                                             ccTTCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTGCCGCACTGACGTGGTGGTTGTTCGGCTc          >  1:54483/1‑143 (MQ=255)
                                                                                                                               ttCAAAATTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCt        >  1:315383/1‑143 (MQ=255)
                                                                                                                                  aaaaTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaa     >  1:162246/1‑143 (MQ=255)
                                                                                                                                     aTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATGTACCGTAGCCAGCCTATTTGCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa  <  1:337508/143‑1 (MQ=255)
                                                                                                                                     aTTGTAGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTgaagaa  <  1:314800/143‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
TTGCGGCAATTTCACGCTGACGCGCTTCGATGGAATCATCACTAAATTCGATAACCTTGCCGCAGTCGAGGCAGATCAGGTGATCGTGGTGATGTTGCTGTGTCAGTTCAAATACGGATTTACCGCCTTCAAAATTGTGGCGGGTGACGATACCAGCGTCGTCAAACTGGTTCAGTACGCGATATACCGTAGCCAGACCAATTTCTTCACCCATATCGATCAGACGTTTGTATAAATCTTCCGCACTGACGTGATGGTTGTCCGGCTCCTGAAGAA  >  NC_000913/710298‑710573

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