Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,525,908 C→G 20.1% intergenic (+17/+22) yfeN → / ← yfeR conserved outer membrane protein YfeN/putative LysR‑type DNA‑binding transcriptional regulator YfeR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,525,9080CG20.1% 41.1 / 3.5 23intergenic (+17/+22)yfeN/yfeRconserved outer membrane protein YfeN/putative LysR‑type DNA‑binding transcriptional regulator YfeR
Reads supporting (aligned to +/- strand):  ref base C (9/9);  new base G (3/2);  total (12/11)
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

TGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTC  >  NC_000913/2525768‑2526038
                                                                                                                                             |                                                                                                                                  
tGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTc                                                                                                                                  >  1:314528/1‑143 (MQ=255)
      gtggtgGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCa                                                                                                                            <  2:114717/143‑1 (MQ=255)
       gtgtggTATGCAACATGACACTTTACTCCGGGGATGCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTCCTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAg                                                                                                                           <  2:390848/141‑1 (MQ=255)
        tgtggAATCCATCATCACAGATTACTACGGGGATTCAGTATCGAGATGCGGATGACAAACTGGGCGAAGCTTTTTCCCATGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGt                                                                                                                          <  2:400517/142‑1 (MQ=255)
         gtggAATGGAACATCACACATTAGTAGGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGAGCCCCGGCAGTa                                                                                                                         <  2:489191/143‑1 (MQ=255)
          tggAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGAGGACAAACTGGGCGAAGATTCTTACCAGGATGCGATAACCTACTACATCACGTTTAACGTCTGATTGATTACTGCCGGCGCTCCCGGCAGTAt                                                                                                                        >  2:336758/1‑143 (MQ=255)
                  aaCATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCAc                                                                                                                <  1:319180/143‑1 (MQ=255)
                             ttACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGCATGACGTCACGCTATATCTGa                                                                                                     >  2:561149/1‑143 (MQ=255)
                               actacGGGGATTCAGTATCGAGATGCGGATTGCAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATa                                                                                                   <  2:310177/143‑1 (MQ=255)
                                 tacGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATCTCTGATACa                                                                                                 >  2:402468/1‑143 (MQ=255)
                                        aTTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATc                                                                                          <  2:570048/143‑1 (MQ=255)
                                                tCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGGTCGCGGCAGTTTGAATTCAGTCTATATCTGTGACAACGGATCCTCTTtcc                                                                                  >  1:124301/1‑143 (MQ=255)
                                                     aggCGGATTACAAACTGGGCGAAGATTGTTACCAGGATGCTATAATCTACTCCATCAAATTTAACTTCTTATTGATTACTGCCGGTGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCa                                                                             <  2:103343/141‑1 (MQ=255)
                                                        ctgATGACAAACTGGGCGAAGAGTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGcc                                                                          <  2:127113/141‑1 (MQ=255)
                                                                         gAAGATTTTTACCGGGATGGGATAATCTTCGCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTGTATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTg                                                         <  2:86615/143‑1 (MQ=255)
                                                                                            tgctgctCTACTCCATCAAATGTTACTTCTGATTGGTTACTGCCGGGGGTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCAt                                      <  2:243577/137‑1 (MQ=255)
                                                                                                       tCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCCCGTCCCATACAGCTCCAgc                           >  2:463400/1‑143 (MQ=255)
                                                                                                              aaTTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGCCTCGTGCCACGCCCCATAACGCTTCGGCGGCGGTg                    >  1:15345/1‑143 (MQ=255)
                                                                                                                            attcatACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAATTGAACGGtt      >  1:99737/6‑143 (MQ=255)
                                                                                                                             ttcatACTGCCGGGAGCCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCGCGTCCCCTAAAGCTTCAGCGGCGGTGGAAAGTGGACGGgtc     >  1:261220/5‑143 (MQ=255)
                                                                                                                               catACTGCCGGGAGCCCCGGCAGTATGAATGCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACATCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGttctt   >  2:505795/3‑143 (MQ=255)
                                                                                                                                ttACTGCCGGGGCTCCCGGCAGTATTAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCCCACCGCCCTGGTCTCGTACCCCCCCCCACAAAGCTTCAGCGGCGGTGGAACGTGAACGGttcttc  >  2:174367/1‑143 (MQ=255)
                                                                                                                                ttACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGttcttc  >  1:180898/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                  
TGCACTGTGGTGGAATGCAACATCACACATTACTACGGGGATTCAGTATCGATATGCGGATGACAAACTGGGCGAAGATTTTTACCAGGATGCGATAATCTACTCCATCAAATTTAACTTCTGATTGATTACTGCCGGGGCTCCCGGCAGTATGAATTCACTCTATATCTGATACAACGGATCCCCTTCCCGCCCAGCCATTAGCGCATTGCCCTGGTCTCGTACCACGTCCCATAAAGCTTCAGCGGCGGTGGAAAGTGAACGGTTCTTC  >  NC_000913/2525768‑2526038

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: