Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 231,008 T→G 8.4% intergenic (‑127/‑114) 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,0080TG 8.4% 124.4 / 2.8 49intergenic (‑127/‑114)yafC/yafDputative LysR family transcriptional regulator YafC/endonuclease/exonuclease/phosphatase domain‑containing protein YafD
Reads supporting (aligned to +/- strand):  ref base T (21/23);  new base G (2/2);  total (23/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 = 6.32e-01

CGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTA  >  NC_000913/230867‑231150
                                                                                                                                              |                                                                                                                                              
cGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGTTTTAATGTTCTCCATAATGAGCAAACTTCTGACCGTTGTATGCACtt                                                                                                                                               >  1:30004/1‑143 (MQ=255)
         tttCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATa                                                                                                                                      >  1:386756/1‑143 (MQ=255)
         tttCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGCAGGATGTCATGTGCTCTATGAGGTGCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATa                                                                                                                                      >  2:784880/1‑143 (MQ=255)
          ttCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACaccg                                                                                                                                     <  2:331925/143‑4 (MQ=255)
          ttCATTTTTGCTCCTGTGGGAAAATAGAAGAGATATTTTGAAGTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACaccg                                                                                                                                     <  1:3562/143‑4 (MQ=255)
                  tGCTCCCGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCCCGCGATTTAATGTTCTCCATAATGAGCGAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATAc                                                                                                                             <  2:595525/143‑1 (MQ=255)
                        tGTTGCAAAATGGAAGAGATCTTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGGCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGc                                                                                                                       <  2:649481/143‑1 (MQ=255)
                         gTTGCAAAATAGAAGAGATATTTTGAATTTATTGGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATCAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCt                                                                                                                      <  2:498261/143‑1 (MQ=255)
                          ttGCAAAATAGAAGCGATATTTTGAATTTATTTGTCATTAAACCTTCAGGATGCGTGTTTTGTCTTGCGTTTTAATAATATCAATAATAAGAAAAATTCTGAACGGTGTAAGCACTTGCTTACATAACAATATACAATTGATc                                                                                                                     >  1:242477/1‑143 (MQ=255)
                          ttGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAACATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAAATGCTc                                                                                                                     >  1:45954/1‑143 (MQ=255)
                          ttGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCGCTTGCTTACATAACCATATACAATTGCTc                                                                                                                     >  1:145803/1‑143 (MQ=255)
                             cagaatagaatagATCTCTTTAATTTTTTTGTCATTAATCCATCAGGAGGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGtt                                                                                                                  <  2:298294/140‑1 (MQ=255)
                                  tgtcagagATATTTTTAATTTATTTGTCAGTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAag                                                                                                             <  2:283202/139‑1 (MQ=255)
                                    gaagaGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAagag                                                                                                           <  2:62763/143‑1 (MQ=255)
                                     aagaGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGt                                                                                                          >  1:656457/1‑143 (MQ=255)
                                     aagaGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGt                                                                                                          <  1:142132/143‑1 (MQ=255)
                                       gagaTATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCCAAATGAGCAAAATTCTGCACGGTGAAATCGCTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGa                                                                                                        >  1:338856/1‑143 (MQ=255)
                                        agaTATTTTGATTTTATTTGTCTTTAATCCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAGTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAg                                                                                                       <  2:19857/143‑1 (MQ=255)
                                        agaTATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAg                                                                                                       >  1:158128/1‑143 (MQ=255)
                                          atatTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGGGATGTCATGCGATTTAATGTTCTCAATAATCAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAAAAATATATAATTGCTCGTTGCAAGAGTGAGCt                                                                                                     >  2:129597/1‑143 (MQ=255)
                                           tatTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTa                                                                                                    <  1:195367/143‑1 (MQ=255)
                                           tatGTTGAATTTATTTGTCATTCAACCCTTTGGAGCACTGTTATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTa                                                                                                    <  1:335910/143‑1 (MQ=255)
                                                 gAATTTATTTGTCACTAAACCATCCGGATGCGTGATATGTCCTGCGATTTAATGTTCTCCATGCTGAGCAAAGTTATGACCGGTGTAAGCACTTGCTTAAATAAGACGCGTTATATTCTCGTTGAAAGAGTGAGCTAAAATcc                                                                                              <  2:612268/143‑1 (MQ=255)
                                                  aaTTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATccc                                                                                             <  1:807533/143‑1 (MQ=255)
                                                          tGTCATTAATCCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATACTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAg                                                                                     <  2:783472/143‑1 (MQ=255)
                                                            tCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATGGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTa                                                                                   <  2:226086/143‑1 (MQ=255)
                                                                                  gATATGTCATGCGATCTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTACTCGTTGAAAGAAGGAGCTAAAAACACCATAACAGTAGAACCCACCCGGGTgggacagg                                                             >  2:680728/1‑136 (MQ=255)
                                                                                     atGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGGAGGGtt                                                          >  2:409308/1‑143 (MQ=255)
                                                                                      ttgcATGCGATTTAATGGTCTCCATAATGAGCAAAATACTGACCGGGGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTg                                                         <  2:167548/140‑1 (MQ=255)
                                                                                      tGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAGATGCACACGTAACAAGATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTg                                                         <  1:801325/143‑1 (MQ=255)
                                                                                      tGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTg                                                         >  1:171362/1‑143 (MQ=255)
                                                                                      tGTCATGCGATTGAATGTGCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTg                                                         <  1:607379/143‑1 (MQ=255)
                                                                                       ggcttGCGATTTCATGTTCTCCAGAATGAGCAAGATTCTGACCGGTGTAAGCACTTGCTTACATAACAATACACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATACCAGTAGAACCCTCCCGAGTGCGGAAGGGTTGa                                                        <  2:388650/139‑1 (MQ=255)
                                                                                              gATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGCAGGGTTGCCGTAAtc                                                 >  1:127803/1‑142 (MQ=255)
                                                                                               aTTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAACTAAAATCCATATAACAGTAGAACCCTCCCGAGGGCGGAAGAGTTGACGTGATag                                                >  2:71150/1‑143 (MQ=255)
                                                                                                  tAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAgg                                             <  1:248726/143‑1 (MQ=255)
                                                                                                   aaTGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGt                                            >  2:587255/1‑143 (MQ=255)
                                                                                                    aTGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAAGATACAATTGCCCGTTGACAGAGTGAGCTAAAATCTCTATAACAGTAGAACCCTCCCGCGTGCGGCAGGGCTGACGAAACAGAggct                                           >  2:523638/1‑141 (MQ=255)
                                                                                                       ttCTCCACAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCa                                        <  1:268021/143‑1 (MQ=255)
                                                                                                         ctcCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCaaa                                      <  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)
                                                                                                                                              ttGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAACTCCCTATAACAGTAGAACGCTCCCGAGTGCGGAAGGGTTGAAGTAATAGAGGTTTCAAAGTCAAAGGAGCGAAAAGAGACCTATGGCATggtcaa  >  2:555624/1‑138 (MQ=255)
                                                                                                                                              |                                                                                                                                              
CGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTA  >  NC_000913/230867‑231150

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: