Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,715,260:1 +T 12.0% coding (61/882 nt) yfiE ← putative LysR‑type DNA‑binding transcriptional regulator YfiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,715,2601.T12.0% 158.5 / 20.9 50coding (61/882 nt)yfiEputative LysR‑type DNA‑binding transcriptional regulator YfiE
Reads supporting (aligned to +/- strand):  ref base . (32/12);  new base T (4/2);  total (36/14)
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.82e-01

CAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCT  >  NC_000913/2715126‑2715400
                                                                                                                                       |                                                                                                                                            
cagcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaa                                                                                                                                      <  2:238303/143‑1 (MQ=255)
 agcaACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaa                                                                                                                                     <  2:217854/143‑1 (MQ=255)
    aaCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAACCTCGCAAAATGGa                                                                                                                                  >  1:767592/1‑143 (MQ=255)
     aCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaattgt                                                                                                                                 <  1:425000/143‑5 (MQ=255)
            tCCTTCGCGGGTAAGGCACACCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCtt                                                                                                                          <  1:915953/143‑1 (MQ=255)
               ttCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAGCAATCCTGCTCAAGCTGGTGGATATGACAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAGAGGAGCCCTCGTcca                                                                                                                       >  1:53845/1‑143 (MQ=255)
                    gggTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACg                                                                                                                  <  1:526960/143‑1 (MQ=255)
                    gggTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACg                                                                                                                  <  1:526923/143‑1 (MQ=255)
                             acaTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCACATAGGAACCCTCTTCCTCCTCGGTTTTTAGc                                                                                                         >  1:346448/1‑143 (MQ=255)
                               aTCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAGAGGAACCCTCTTCCACCACGGCTTTAAGCGt                                                                                                       >  2:97245/1‑143 (MQ=255)
                                tCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGATTACAAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGta                                                                                                      <  1:908975/143‑1 (MQ=255)
                                         cAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCCCACAAAGGAACCCTCTTCCACCACGGTTTTAAGACTAATAAAACga                                                                                             >  1:788685/1‑142 (MQ=255)
                                                   aaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGGAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATccg                                                                                   >  1:803111/1‑142 (MQ=255)
                                                   aaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCa                                                                                   >  1:588068/1‑143 (MQ=255)
                                                    aaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGACGCTCGCAACAAGGAACCCTCCTCCACCACGGTTTTAAGCGTAATAACACGGCGCAGATCCac                                                                                  >  1:258972/1‑142 (MQ=255)
                                                     aTAACTGCACTGAAACCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATa                                                                                 <  1:679551/143‑1 (MQ=255)
                                                         cTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAccc                                                                             >  2:78444/1‑143 (MQ=255)
                                                          tGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAACAGGAACCCTCTTCCACCACGGTTTTACGCGTAATAAAACGGCGCAGAGCCACAAcccc                                                                            >  2:192700/1‑143 (MQ=255)
                                                               tGAAAACTCCTGCTCAAGCTGCTGAATATGAACAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGAAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGACGCAGATCAATAACCCCAAAcc                                                                       >  1:320944/1‑143 (MQ=255)
                                                               tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAAcc                                                                       >  1:464074/1‑143 (MQ=255)
                                                                gAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCt                                                                      <  1:633595/143‑1 (MQ=255)
                                                                 aaaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCAACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTa                                                                     >  1:786918/1‑143 (MQ=255)
                                                                      tCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACATCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAAATCCATAACCCCAAACCTATCGaa                                                                >  1:715994/1‑143 (MQ=255)
                                                                      tCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGTAAAAAGGAACCCTCTTCCACCACGGTTTTCAGCGTACTAACACGGCGCAGATCCATACCGCCAAAACCATCaaa                                                                >  1:750919/1‑143 (MQ=255)
                                                                         tGCTCAAGCGGCTGAATATGAACAGTCACCGTCGATTGCGTACAGCCCAATTTTTGCGTAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAat                                                             <  2:438010/143‑1 (MQ=255)
                                                                           cTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAatat                                                           <  1:721623/143‑1 (MQ=255)
                                                                           cTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGATTTAAGCGTAATAAAACGGCGCAGATCCATAACCCGAAACCTCTCGAAAatat                                                           >  2:628845/1‑143 (MQ=255)
                                                                            tCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATc                                                          >  1:675588/1‑143 (MQ=255)
                                                                             cAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCg                                                         >  1:221973/1‑143 (MQ=255)
                                                                                  tgctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATCAACCACAAGATATCCATAACCCCAAACCTATCGAAAACATCGAATca                                                    >  2:886481/1‑142 (MQ=255)
                                                                                    ctgAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCGACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAAAATCGAATCGAg                                                  >  1:545816/1‑143 (MQ=255)
                                                                                       aaTATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCATCGTCCGCAACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTAGCGAAGATATCGAATCTAGAat                                               >  2:493397/1‑143 (MQ=255)
                                                                                           tGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaa                                           >  2:265784/1‑143 (MQ=255)
                                                                                            gAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAACTATCGAATCTAGAATATaaa                                          >  1:5886/1‑143 (MQ=255)
                                                                                              aaaGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGGTTTAAGCGTAATAAAACGGCGCCGATCCATAACCCCAAACCTCTCGAAAAAATCGAATCTACAATATaaaaa                                        >  2:313484/1‑143 (MQ=255)
                                                                                               aaGTCACCGTCGATTGTGTACAGCACAATTTTAGCGAAGCTCGCAAGAGGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAc                                       <  1:75042/143‑1 (MQ=255)
                                                                                                  tCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGTACATATCGAATCTAGAACACAAAAAcatt                                    >  1:809819/1‑143 (MQ=255)
                                                                                                        tCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAtttt                              >  1:563748/1‑143 (MQ=255)
                                                                                                            ttGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTGATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATCTTTTTa                          >  1:521068/1‑143 (MQ=255)
                                                                                                              gtgtACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACTGCTCATATCCATAACCCCAAACCTATCGAAAATAGCGAATCTAGAATATAAAAACATTCATTTTTTTaaa                        <  1:105615/143‑1 (MQ=255)
                                                                                                                      cacaATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTACACTACAGTTTTAAGCGTAATAAAACGGCGCAGATACATAACCCCAAAACTATCGAAAATATCGAAGTGAGAATATAAACACAAACATTTTATTAAATCTTCCgt                >  2:239368/1‑143 (MQ=255)
                                                                                                                         aaTTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTc             >  1:244176/1‑143 (MQ=255)
                                                                                                                         aaTTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTc             >  1:843069/1‑143 (MQ=255)
                                                                                                                         aaTTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCCATCCAAAACATCGAATCAAGAATATAAACACATTCCTTTTCTAAAATGTTCCatggc             >  2:668028/1‑141 (MQ=255)
                                                                                                                           tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg           >  1:683239/1‑143 (MQ=255)
                                                                                                                            ttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACTCCAAACCTATCGAAAATATCGAATCTAGAATATAACAACATTCCTTTTTTTAAATGTTCCGTCTCggg          >  2:291295/1‑143 (MQ=255)
                                                                                                                            ttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATCAAAACACTCATTTTTTTAAATGTTCCGGGTCggg          >  2:769731/1‑143 (MQ=255)
                                                                                                                            ttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCggg          >  1:713597/1‑143 (MQ=255)
                                                                                                                                   aaGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATAGAAAATAGAATTTAAAAACATACATATTTTTAAATGTTCCGTGTCGGGTACTGTc   >  1:291062/1‑143 (MQ=255)
                                                                                                                                    aGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCt  <  1:599819/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                            
CAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCT  >  NC_000913/2715126‑2715400

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