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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,0070CA 7.0% 194.5 / 5.2 57intergenic (‑126/‑115)yafC/yafDLysR family putative transcriptional regulator/endo/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (22/31);  new base A (2/2);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

TTCTTCCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTA  >  NC_000913/230861‑231150
                                                                                                                                                  |                                                                                                                                               
ttcttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTTTGACCGGTGTAAGCACtt                                                                                                                                               >  2:504065/1‑149 (MQ=255)
  cttcCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGc                                                                                                                                             <  4:184322/149‑1 (MQ=255)
      cGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTAc                                                                                                                                         <  1:232880/149‑1 (MQ=255)
         cGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACacc                                                                                                                                      <  2:281236/149‑3 (MQ=255)
          gTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAGCCATCAGGATGCGTGATATGGCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGGTTACATaa                                                                                                                                     <  3:285522/149‑1 (MQ=255)
          gTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATaa                                                                                                                                     <  2:240755/149‑1 (MQ=255)
          gTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATaa                                                                                                                                     <  2:65911/149‑1 (MQ=255)
           tGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAAc                                                                                                                                    <  4:232227/149‑1 (MQ=255)
           tGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCAAGTGCTTACaccgg                                                                                                                                    <  4:227757/149‑5 (MQ=255)
              cTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAat                                                                                                                                 <  2:71346/149‑1 (MQ=255)
              cTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAat                                                                                                                                 <  2:431735/149‑1 (MQ=255)
               tttCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAata                                                                                                                                <  4:322538/149‑1 (MQ=255)
               tttCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAata                                                                                                                                >  3:564459/1‑149 (MQ=255)
                 tCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAAtata                                                                                                                              <  3:512400/149‑1 (MQ=255)
                 tCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAAtata                                                                                                                              <  1:192102/149‑1 (MQ=255)
                    tttttGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACaa                                                                                                                           <  4:483988/149‑1 (MQ=255)
                          cTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTc                                                                                                                     <  4:213674/149‑1 (MQ=255)
                           tCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCg                                                                                                                    <  4:516183/149‑1 (MQ=255)
                             cTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGtt                                                                                                                  <  3:556815/149‑1 (MQ=255)
                             cTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGtt                                                                                                                  <  1:186950/149‑1 (MQ=255)
                                ttGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGaa                                                                                                               >  2:282352/1‑149 (MQ=255)
                                 tGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGaaa                                                                                                              <  4:225171/149‑1 (MQ=255)
                                      aaTAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCCCGTTGAAAGAGTg                                                                                                         <  2:273649/149‑1 (MQ=255)
                                       aTAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGa                                                                                                        <  2:292279/149‑1 (MQ=255)
                                       aTAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGa                                                                                                        <  2:66123/149‑1 (MQ=255)
                                                atatTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATc                                                                                               >  2:226457/1‑149 (MQ=255)
                                                 tatTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATcc                                                                                              <  4:578363/149‑1 (MQ=255)
                                                  atTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATccc                                                                                             <  2:154/149‑1 (MQ=255)
                                                   ttttGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCt                                                                                            >  1:196003/1‑149 (MQ=255)
                                                          tttatttGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAg                                                                                     >  1:471958/1‑149 (MQ=255)
                                                            tatttGTCATTAAACCATCAGGGTGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCGAAATTCTGACCGGTGTAAACACTTGCTTACATAACAATATACAATTGCTCGTTGAACGAGTGAGCTAAAATCCCTATAACAGTa                                                                                   >  2:272332/1‑149 (MQ=255)
                                                                   cATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCt                                                                            <  1:355141/149‑1 (MQ=255)
                                                                     ttAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTcc                                                                          >  2:373297/1‑149 (MQ=255)
                                                                            aTCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGc                                                                   >  2:554881/1‑149 (MQ=255)
                                                                                                 tGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATagag                                              >  1:136175/1‑149 (MQ=255)
                                                                                                         aaTGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCaaa                                      >  1:66392/1‑149 (MQ=255)
                                                                                                           tGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGt                                    >  3:202343/1‑149 (MQ=255)
                                                                                                            gTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTc                                   >  2:211489/1‑149 (MQ=255)
                                                                                                             ttCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCa                                  <  1:226457/149‑1 (MQ=255)
                                                                                                                  cATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGt                             >  2:582123/1‑149 (MQ=255)
                                                                                                                     aaTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACTTAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCg                          >  2:241093/1‑149 (MQ=255)
                                                                                                                     aaTGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCg                          >  2:485867/1‑149 (MQ=255)
                                                                                                                      aTGAGCAAAATACTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGa                         >  3:514070/1‑149 (MQ=255)
                                                                                                                       tGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGaa                        <  2:196003/149‑1 (MQ=255)
                                                                                                                            aaaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAAtata                                                                                                                              <  2:308825/42‑1 (MQ=255)
                                                                                                                            aaaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAAtata                                                                                                                              >  1:308825/1‑42 (MQ=255)
                                                                                                                            aaaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCaaa                                      <  4:529111/130‑1 (MQ=255)
                                                                                                                            aaaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCaaa                                      >  3:529111/1‑130 (MQ=255)
                                                                                                                            aaaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAac                   >  1:174895/1‑149 (MQ=255)
                                                                                                                             aaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAaca                  <  3:380475/149‑1 (MQ=255)
                                                                                                                             aaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAaca                  <  1:584496/149‑1 (MQ=255)
                                                                                                                             aaaTTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAaca                  <  1:504065/149‑1 (MQ=255)
                                                                                                                                        tatgtAAGCAAGTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATg       >  1:51738/3‑149 (MQ=255)
                                                                                                                                         atgtAAGCAAGTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATgc      >  1:504813/2‑149 (MQ=255)
                                                                                                                                             aaGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTa  >  2:305318/1‑149 (MQ=255)
                                                                                                                                                 aCTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAac                   >  2:241177/1‑128 (MQ=255)
                                                                                                                                                 aCTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAac                   <  1:241177/128‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TTCTTCCGACGTGGCTTTCATTTTTGCTCCTGTTGCAAAATAGAAGAGATATTTTGAATTTATTTGTCATTAAACCATCAGGATGCGTGATATGTCATGCGATTTAATGTTCTCCATAATGAGCAAAATTCTGACCGGTGTAAGCACTTGCTTACATAACAATATACAATTGCTCGTTGAAAGAGTGAGCTAAAATCCCTATAACAGTAGAACCCTCCCGAGTGCGGAAGGGTTGACGTAATAGAGGTTTCAAAGTCAAAAGTGCGAAAAAACACCTATGCCATGCGCTA  >  NC_000913/230861‑231150

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑