Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,007 C→A 8.2% intergenic (‑126/‑115) yafC ← / → yafD putative LysR family transcriptional regulator YafC/endonuclease/exonuclease/phosphatase domain‑containing protein YafD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,0070CA 8.2% 131.1 / 3.9 50intergenic (‑126/‑115)yafC/yafDputative LysR family transcriptional regulator YafC/endonuclease/exonuclease/phosphatase domain‑containing protein YafD
Reads supporting (aligned to +/- strand):  ref base C (22/23);  new base A (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

CGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATG  >  NC_000913/230867‑231140
                                                                                                                                             |                                                                                                                                     
cGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGTTTTAATGTTCT‑CCATAATGAGCAAACTTCTGACCGTTGTATGCACtt                                                                                                                                     >  1:30004/1‑143 (MQ=255)
         tttCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATa                                                                                                                            >  1:386756/1‑143 (MQ=255)
         tttCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGCAGGATGTCATGTGCTCT‑ATGAGGTGCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATa                                                                                                                            >  2:784880/1‑143 (MQ=255)
          ttCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACaccg                                                                                                                           <  2:331925/143‑4 (MQ=255)
          ttCATTTTTGCTCCTGTGGGAAAATAGAAGAGATATTTTGAAGTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACaccg                                                                                                                           <  1:3562/143‑4 (MQ=255)
                  tGCTCCCGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCCCGCGATTTAATGTTCT‑CCATAATGAGCGAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATAc                                                                                                                   <  2:595525/143‑1 (MQ=255)
                        tGTTGCAAAATGGAAGAGATCTTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGGCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGc                                                                                                             <  2:649481/143‑1 (MQ=255)
                         gTTGCAAAATAGAAGAGATATTTTGAATTTATTGGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATCAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCt                                                                                                            <  2:498261/143‑1 (MQ=255)
                          ttGCAAAATAGAAGCGATATTTTGAATTTATTTGTCATTAAACCTTCAGGATGCGTGTTTTGTCTTGCGTTTTAATAATAT‑CAATAATAAGAAAAATTCTGAACGGTGTAAGCACTTGCTTACATAACAATATACAATTGATc                                                                                                           >  1:242477/1‑143 (MQ=255)
                          ttGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAACATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAAATGCTc                                                                                                           >  1:45954/1‑143 (MQ=255)
                          ttGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCGCTTGCTTACATAACCATATACAATTGCTc                                                                                                           >  1:145803/1‑143 (MQ=255)
                             cagaatagaatagATCTCTTTAATTTTTTTGTCATTAATCCATCAGGAGGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGtt                                                                                                        <  2:298294/140‑1 (MQ=255)
                                  tgtcagagATATTTTTAATTTATTTGTCAGTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAag                                                                                                   <  2:283202/139‑1 (MQ=255)
                                    gaagaGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAagag                                                                                                 <  2:62763/143‑1 (MQ=255)
                                     aagaGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGt                                                                                                >  1:656457/1‑143 (MQ=255)
                                     aagaGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGt                                                                                                <  1:142132/143‑1 (MQ=255)
                                       gagaTATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCCAAATGAGCAAAATTCTGCACGGTGAAATCGCTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGa                                                                                              >  1:338856/1‑143 (MQ=255)
                                        agaTATTTTGATTTTATTTGTCTTTAATCCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAGTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAg                                                                                             <  2:19857/143‑1 (MQ=255)
                                        agaTATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAg                                                                                             >  1:158128/1‑143 (MQ=255)
                                          atatTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGGGATGTCATGCGATTTAATGTTCT‑CAATAATCAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAAAAATATATAATTGCTCGTTGCAAGAGTGAGCt                                                                                           >  2:129597/1‑143 (MQ=255)
                                           tatTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTa                                                                                          <  1:195367/143‑1 (MQ=255)
                                           tatGTTGAATTTATTTGTCATTCAACCCTTTGGAGCACTGTTATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTa                                                                                          <  1:335910/143‑1 (MQ=255)
                                                 gAATTTATTTGTCACTAAACCATCCGGATGCGTGATATGTCCTGCGATTTAATGTTCT‑CCATGCTGAGCAAAGTTATGACCGGTGTAAGCACTTGCTTAAATAAGACGCGTTATATTCTCGTTGAAAGAGTGAGCTAAAATcc                                                                                    <  2:612268/143‑1 (MQ=255)
                                                  aaTTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATccc                                                                                   <  1:807533/143‑1 (MQ=255)
                                                          tGTCATTAATCCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATACTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAg                                                                           <  2:783472/143‑1 (MQ=255)
                                                            tCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATGGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTa                                                                         <  2:226086/143‑1 (MQ=255)
                                                                                  gATATGTCATGCGATCTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTACTCGTTGAAAGAAGGAGCTAAAAACACCATAACAGTAGAACCCACCCGGGTgggacagg                                                   >  2:680728/1‑136 (MQ=255)
                                                                                     atGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGGAGGGtt                                                >  2:409308/1‑143 (MQ=255)
                                                                                      ttgcATGCGATTTAATGGTCT‑CCATAATGAGCAAAATACTGACCGGGGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTg                                               <  2:167548/140‑1 (MQ=255)
                                                                                      tGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCAGATGCACACGTAACAAGATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTg                                               <  1:801325/143‑1 (MQ=255)
                                                                                      tGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTg                                               >  1:171362/1‑143 (MQ=255)
                                                                                      tGTCATGCGATTGAATGTGCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTg                                               <  1:607379/143‑1 (MQ=255)
                                                                                       ggcttGCGATTTCATGTTCT‑CCAGAATGAGCAAGATTCTGACCGGTGTAAGCACTTGCTTACATAACAATACACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATACCAGTAGAACCCTCCCGAGTGCGGAAGGGTTGa                                              <  2:388650/139‑1 (MQ=255)
                                                                                              gATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGCAGGGTTGCCGTAAtc                                       >  1:127803/1‑142 (MQ=255)
                                                                                               aTTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAACTAAAATCCATATAACAGTAGAACCCTCCCGAGGGCGGAAGAGTTGACGTGATag                                      >  2:71150/1‑143 (MQ=255)
                                                                                                  tAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAgg                                   <  1:248726/143‑1 (MQ=255)
                                                                                                   aaTGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGt                                  >  2:587255/1‑143 (MQ=255)
                                                                                                    aTGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAAGATACAATTGCCCGTTGACAGAGTGAGCTAAAATCTCTATAACAGTAGAACCCTCCCGCGTGCGGCAGGGCTGACGAAACAGAggct                                 >  2:523638/1‑141 (MQ=255)
                                                                                                       ttCT‑CCACAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCa                              <  1:268021/143‑1 (MQ=255)
                                                                                                         ct‑cCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCaaa                            <  2:343846/143‑1 (MQ=255)
                                                                                                               tAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCaa                       <  1:784880/143‑1 (MQ=255)
                                                                                                                aaTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGCGTCAAAGCCaca                      >  2:504794/1‑141 (MQ=255)
                                                                                                                    aGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTg                  <  1:600659/143‑1 (MQ=255)
                                                                                                                    aGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTg                  >  2:707615/1‑143 (MQ=255)
                                                                                                                      cAAAATTCTGACCGGTGTAAGTACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACACTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCg                >  1:358803/1‑143 (MQ=255)
                                                                                                                       aaaaTTCTGACCGGTGTAAGCACTTGCTTACATCACAATCTACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAGGTGCGa               >  1:44301/1‑143 (MQ=255)
                                                                                                                         aaTTCTGACGGGTGTAAGCACTGGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaaa             <  2:759527/143‑1 (MQ=255)
                                                                                                                           ttCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGATGACGCAATAGAGGTTTCAAAGTCCGAGGTGCGacaca           >  1:105715/1‑139 (MQ=255)
                                                                                                                                   tatgtAAGCAAGTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTAt   >  1:93221/3‑143 (MQ=255)
                                                                                                                                    atgtAAGCAAGTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTGAGAGAGGTTTCAAAGTCGAAAGTGCGAAAAGACACCTATg  >  1:206351/2‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                     
CGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCT‑CCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATG  >  NC_000913/230867‑231140

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: ‑